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Abstract 

This paper presents a brief review of the newly developed Extended Electrodynamics. The 
relativistic and non-relativistic approaches to the extension of Maxwell equations are considered 
briefly, and the further study is carried out in relativistic terms in Minkowski space-time. The non- 
linear vacuum solutions are considered and fully described. It is specially pointed out that solitary 
waves with various, in fact arbitrary, spatial structure and photon- like propagation properties exist. 
The null character of all non-linear vacuum solutions is established and extensively used further. 
Coordinate-free definitions are given to the important quantities amplitude and phase. The new 
quantity, named scale factor, is introduced and used as a criterion for availability of rotational 
component of propagation of some of the nonlinear, i.e. nonmaxwellian, vacuum solutions. The 
group structure properties of the nonlinear vacuum solutions are analyzed in some detail, showing 
explicitly the connection of the vacuum solutions with some complex valued functions. Connection- 
curvature interpretations are given and a special attention is paid to the curvature interpretation of 
the intrinsic rotational (spin) properties of some of the nonlinear solutions. Several approaches to 
coordinate-free local description and computation of the integral spin momentum are considered. 
Finally, a larage family of nonvacuum spatial soliton-like solutions is explicitly written down, and 
a procedure to get (3+1) versions of the known (1+1) soliton solutions is obtained. 
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1 Introduction 



1.1 The concept of Physical Object 

When we speak about free physical objects, e.g. classical particles, solid bodies, elementary particles, 
etc., we always keep in mind that, although we consider them as free, they can not in principle be 
absolutely free. What is really understood under "free object" (or more precisely "relatively free 
object") is, that some definite properties (e.g. mass, velocity) of the object under consideration do not 
change in time during its evolution under the influence of the existing environment. The availability 
of such time-stable features of any physical object guarantees its identification during its existence 
in time. Without such an availability of constant in time properties (features), which are due to the 
object's resistance abilities, we could not speak about objects and knowledge at all. So, for example, 
a classical mass particle in external gravitational field is free with respect to its mass, and it is not free 
with respect to its behavior as a whole (its velocity changes), because in classical mechanics formalism 
its mass does not change during the influence of the external field on its accelerated way of motion. 

The above view implies that three kinds of quantities will be necessary to describe as fully as 
possible the existence and the evolution of a given physical object: 

1. Proper (identifying) characteristics, i.e. quantities which do NOT change during the 
entire existence of the object. The availability of such quantities allows to distinguish a physical 
object among the other ones. 

2. Kinematical characteristics, i.e. quantities, which describe the allowed space-time evolution, 
where "allowed" means consistent with the constancy of the identifying characteristics. 

3. Dynamical characteristics, i.e. quantities which are functions (explicit or implicit) of the 
proper and of the kinematical characteristics. 

Some of the dynamical characteristics have the following two important properties: they are uni- 
versal, i.e. every physical object carries nonzero value of them (e.g. energy-momentum), and they 
are conservative, i.e. they may just be transferred from one physical object to another (in various 
forms) but no loss is allowed. 

Hence, the evolution of a physical object subject to bearable/acceptable exterior influence (per- 
turbation), coming from the existing environment, has three aspects: 

1. constancy of the proper (identifying) characteristics, 

2. allowed kinematical evolution, and 

3. exchange of dynamical quantities with the physical environment. 

If the physical object under consideration is space-extended (continuous) and is described by 
a many-component mathematical object we may talk about internal exchange of some dynamical 
characteristic among the various components of the object. For example, the electromagnetic field has 
two vector components (E, B) and internal energy-momentum exchange between E and B should be 
considered as possible. 

The third feature suggests that the dynamical equations, describing locally the evolution of the 
object, may come from giving an explicit form of the quantities controlling the local internal and 
external exchange processes i.e. writing down corresponding local balance equations. Hence, denoting 
the local quantities that describe the external exchange processes by Qi, i = 1, 2, . . . , the object should 
be considered to be Qi-free, i = 1, 2, . . . , if the corresponding integral values are time constant, which 
can be achieved only if Qi obey differential equations presenting appropriately (implicitly or explicitly) 
corresponding local versions of the conservation laws (continuity equations). In the case of absence 
of external exchange these equations should describe corresponding internal exchange processes. The 
corresponding evolution in this latter case may be called proper evolution. 

In trying to formalize these views we have to give some initial explicit formulations of some most 
basic features (properties) of what we call physical object, which features would lead us to a, more or 
less, adequate theoretical notion of our intuitive notion of a physical object. Anyway, the following 
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properties of the theoretical concept "physical object" we consider as necessary: 

1. It can be created. 

2. It can be destroyed. 

3. It occupies finite 3-volumes at any moment of its existence, so it has structure. 

4. It has a definite stability to withstand some external disturbances. 

5. It has definite conservation properties. 

6. It necessarily carries energy-momentum, and, possibly, other measurable (conservative or non- 
conservative) physical quantities. 

7. It exists in an appropriate environment (called usually vacuum), which provides all necessary 
existence needs. 

8. It can be detected by the rest of the world through allowed energy-momentum exchanges. 

9. It may combine with other appropriate objects to form new objects of higher level structure. 

10. Its death gives necessarily birth to new objects following definite rules of conservation. 

Remark. The property to be finite we consider as a very essential one. So, the above features do NOT 
allow the classical material points and the infinite classical fields (e.g. plane waves) to be considered 
as physical objects since the former have no structure and cannot be destroyed, and the latter carry 
infinite energy, so they cannot be created. Hence, the Born-Infeld "principle of finiteness" stating that 
a satisfactory theory should avoid letting physical quantities become infinite, declared in [9], may be 
strengthened as follows: all real physical objects are spatially finite entities and NO infinite values of 
the physical quantities carried by them are allowed. 

Clearly, together with the purely qualitative features physical objects carry important quantita- 
tively described physical properties, and any external interaction may be considered as an exchange 
of such quantities provided both the object and the environment carry them. Hence, the more uni- 
versal is a physical quantity the more useful for us it is, and this moment determines the exclusively 
important role of energy- momentum, which modern physics considers as the most universal one, i.e. 
every physical object necessarily carries energy-momentum. 

From pure formal point of view a description of the evolution of a given continuous physical object 
must include obligatory two mathematical objects: 

1. First, the mathematical object which is meant to represent as fully as possible the integrity 
of the object under consideration, and when subject to appropriate operators, must reproduce 
explicitly all important information about the structure and admissible dynamical evolution of the 
physical object; 

2. Second, the mathematical object D*$> which represents the admissible changes, where D 
is appropriately chosen differential operator acting mainly on the kinematical and dynamical char- 
acteristics. If D does not depend on \I> and its derivatives, and, so, on the corresponding proper 
characteristics of the relation Dty = then would mean that those kinematical and dynamical 
properties of which feel the action of D, are constant with respect to D, so the evolution prescribed 
by = 0, would have "constant" character and would say nothing about possible changes of those 
characteristics, which do not feel D. 

From principle point of view, it does not seem so important if the changes are zero, or not 
zero, the significantly important point is that the changes D^> are admissible. The changes 
are admissible in the following two cases: first, when related to the very object through some 
"projection" P upon the "projections" P(D^) vanish, and then the object is called free (with 
respect to those characteristics which feel D); second, when the projections P(D^) do not vanish, 
but the object still survives, then the object is called not free (with respect to the same characteristics). 
In the first case the corresponding admissible changes D^> have to be considered as having an intrinsic 
for the object nature, and they should be generated by some necessary for the very existence of the 
object internal energy-momentum redistribution during evolution (recall point 6 of the above stated 
10 properties of a physical object). In the second case the admissible changes D^, in addition to 
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the intrinsic factors, depend also on external factors, so, the corresponding projections P(D^) should 
describe explicitly or implicitly some energy-momentum exchange with the environment. 

Following this line of considerations we come to a conclusion that every description of a free 
physical object must include some mathematical expression of the kind ¥(^,D^;S) = 0, specifying 
(through the additional quantities S) what and how changes, and specifying also what is projected 
and how it is projected. If the object is not free but it survives when subject to the external influ- 
ence, then it is very important the quantity F(\l/, D^; S) / to have, as much as possible, universal 
character and to present a change of a conservative quantity, so that this same quantity to be ex- 
pressible through the characteristics T of the external object(s). Hence, specifying differentially some 
conservation/balance properties of the object under consideration, and specifying at every space- 
time point the corresponding admissible exchange process with the environment through the equation 
F(^, D^;Q) = G(F, dT; D^, ...) we obtain corresponding equations of motion being consistent 
with the corresponding integral conservation properties. 

In accordance with these considerations it seems reasonable to note also that this two-sided 
change- conservation nature of a physical object could be mathematically interpreted in terms of the 
integrability-nonintegrability properties of a system of PDE. In fact, according to the Frobenius inte- 
grability theorems (see Sec. 9), integrability means that every Lie bracket of two vector fields of the 
corresponding differential system is linearly expressible through the vector fields generating the very 
system, so, the corresponding curvature is zero and nothing flows out of the system. From physi- 
cal point of view this would mean that such physical objects are difficult to study, because, strictly 
speaking, they seem to be strictly isolated, and, therefore, unobservable. When subject to appropriate 
external perturbation they lose isolation, so they are no more integrable, and have to transform to 
other physical objects, which demonstrate new integrability properties. 

On the contrary, nonintegrability means nonzero curvature, so there are outside directed Lie brack- 
ats [X, Y] (called usually vertical projections) of the generators of the differential system, and the 
correspondig flows of these Lie brackets carry the points of the expected integral manifold out of it, 
so, something is possible to flow out of the system. These sticking out Lie brackets look like system's 
tools, or "tentakles", which carry out the necessary for its very existence contacts with the rest of the 
world. We can say that namely these "nonintegrability" features of the physical system makes it ac- 
cessible to be studied. Figuratively speaking, "the nonintegrability properties protect and support the 
integrability properties" of the system, so these both kind of properties are important for the system's 
existence, and their coavailability/coexistence demonstrates the dual nature of physical systems. 

We note that this viewpoint assumes somehow, that for an accessible to study time-stable physical 
system it is in principle expected both, appropriate integrable and nonintegrable differential/Pfaff 
systems to be associated. And this fits quite well with the structure of a fiber bundle: as we know, 
the vertical subbundle of a fiber bundle is always integrable with integral manifolds the corresponding 
fibers; a horizontal subbundle chosen may carry nonzero curvature, and this curvature characterizes 
the system's properties that make possible some local studies to be performed. 

We recall that this idea of change-conservation has been used firstly by Newton in his momentum 
balance equation p = F which is the restriction of the nonlinear partial differential system V p p = mF, 
or pWipi = mF- 7 , on some trajectory. This Newton's system of equations just says that there are 
physical objects in Nature which admit the "point-like" approximation, and which can exchange 
energy-momentum with "the rest of the world" but keep unchanged some other of their properties, 
and this allows these objects to be identified in space-time and studied as a whole, i.e. as point-like 
ones. The corresponding "tentakles" of these objects seem to carry out contacts with the rest of the 
world by means of energy-momentum exchange. 

In macrophysics, as a rule, the external influences are such that they do NOT destroy the sys- 
tem, and in microphysics a full restructuring is allowed: the old ingredients of the system may fully 
transform to new ones, (e.g. the electron-positron annihilation) provided the energy-momentum con- 
servation holds. The essential point is that whatever the interaction is, it always results in appearing of 
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relatively stable objects, carrying energy-momentum and some other particular physically measurable 
quantities. This conclusion emphasizes once again the importance of having an adequate notion of 
what is called a physical object, and of its appropriate mathematical representation. 

Modern science requires a good adequacy between the real objects and the corresponding math- 
ematical model objects. So, the mathematical model objects * must necessarily be spatially finite, 
and even temporally finite if the physical object considered has by its intrinsic nature finite life-time. 
This most probably means that ^ must satisfy nonlinear partial differential equation(s), which should 
define in a consistent way the admissible changes and the conservation properties of the object un- 
der consideration. Hence, talking about physical objects we mean time-stable spatially finite entities 
which have a well established balance between change and conservation, and this balance is kept by a 
permanent and strictly fixed interaction with the environment. 

This notion of a (finite continuous) physical object sets the problem to try to consider and represent 
its integral characteristics through its local dynamical properties, i.e. to try to understand its nature 
and integral appearance as determined and caused by its dynamical structure. The integral appearance 
of the local features may take various forms, in particular, it might influence the spatial structure 
of the object. In view of this, the propagational behavior of the object as a whole, considered in 
terms of its local translational and rotational components of propagation, which, in turn, should be 
related to the internal energy-momentum redistribution during propagation, could be stably consistent 
only with some distinguished spatial structures. We note that the local rotational components of 
propagation NOT always produce integral rotation of the object, but if they are available, time stable, 
and consistent with the translational components of propagation, some specific conserved quantity 
should exist. And if the rotational component of propagation shows some consistent with the object's 
spatial structure periodicity, clearly, the corresponding frequency may be used to introduce such a 
quantity. 

Finally we note that the rotational component of propagation of a (continuous finite) physical 
object may be of two different origins: relative and intrinsic. In the "relative" case the corresponding 
physical quantity, called angular momentum, depends on the choice of some external to the system 
factors, usually these are relative axis and relative point. In the "intrinsic" case the rotational com- 
ponent (if it is not zero) is meant to carry intrinsic information about the system considered, so, the 
corresponding physical quantity, usually called spin, should NOT depend on any external factors. 

As an idealized (mathematical) example of an object as outlined above, let's consider a spatial 
region D of one-step piece of a helical cylinder with some proper (or internal) diameter r , and let 
D be winded around some straightline axis Z. Let at some (initial) moment t Q our mathematical 
model-object \I/ be different from zero only inside D. Let now at t > t Q the object \E', i.e. the region D, 
begin moving as a whole along the helical cylinder with some constant along Z (translational) velocity 
c in such a way that every point of D follows its own (helical) trajectory around Z and never crosses 
the (helical) trajectory of any other point of D. Obviously, the rotational component of propagation 
is available, but the object does NOT rotate as a whole. Moreover, since the translational velocity c 
along Z is constant, the spatial periodicity A, i.e. the height of D along Z, should be proportional to 
the time periodicity T, and for the corresponding frequency v = 1/T we obtain v = c/A. Clearly, this 
is an idealized example of an object with a space-time consistent dynamical structure, but any physical 
interpretation would require to have explicitly defined D, corresponding dynamical equations, local 
and integral conserved quantities, A, c and, probably, some other parameters. We shall show that in 
the frame of EED such finite solutions of photon-like nature do exist and can be explicitly written 
(see the pictures on pp. 62-63). 
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1.2 Notes on Faraday-Maxwell Electrodynamics 

The 19th century physics, due mainly to Faraday and Maxwell, created the theoretical concept of 
classical field as a model of a class of spatially continuous (extended) physical objects having dynamical 
structure. As a rule, the corresponding free fields are considered to satisfy linear dynamical equations, 
and the corresponding time dependent solutions are in most cases infinite. The most important 
examples seem to be the free electromagnetic fields, which satisfy the charge-free Maxwell equations. 
The concepts of flux of a vector field through a 2- dimensional surface and circulation of a vector field 
along a closed curve were coined and used extensively. Maxwell's equations in their integral form 
establish where the time-changes of the fluxes of the electric and magnetic fields go to, or come from, 
in both cases of a closed 2-surface and a 2-surface with a boundary. We note that these fluxes are 
specific to the continuous character of the physical object under consideration and it is important 
also to note that Maxwell's field equations have not the sense of direct energy-momentum 
balance relations as the Newton's law p = F has. Nevertheless, they are consistent with 
energy-momentum conservation in finite regions, as it is well known, but they impose much stronger 
requirements on the field components, so, the class of admissible solutions is much more limited 
than the local conservation equations would admit. Moreover, from qualitative point of view, the 
electromagnetic induction experiments show that a definite quantity of the field energy-momentum 
is transformed to mechanical energy- momentum carried away by the charged particles, but, on the 
contrary, the corresponding differential equation does not allow such energy-momentum transfer (see 
further). Hence, the introduced by Newton basic approach to derive dynamics from local energy- 
momentum balance equations had been left off, and this step is still respected today. 

Talking about description of free fields by means of differential equations we mean that these 
equations must describe locally the intrinsic dynamics of the field, so NO boundary conditions or any 
other external influence should be present. The pure field Maxwell equations (although very useful 
for considerations in finite regions with boundary conditions) have time-dependent vacuum solutions 
(i.e free-field solutions, so, solutions in the whole space) that give an inadequate description of the 
real free fields. As a rule, if these solutions are time-stable, they occupy the whole 3-space or an 
infinite subregion of it, and they do not go to zero at infinity, hence, they carry infinite energy and 
momentum. As an example we recall the (transverse with zero invariants) plane wave solution, given 
by the electric and magnetic fields of the form 

E = ju(ci + ez),p(ct + ez), j; B = jep(rf + ez), -eu(ct + ez), oj , e = ±1, 

where u and p are arbitrary differentiable functions. Even if u and p are soliton-like with respect to 
the coordinate z, they do not depend on the other two spatial coordinates (x,y). Hence, the solution 
occupies the whole M 3 , or its infinite subregion, and clearly it carries infinite integral energy 

E= (u 2 + p 2 )dxdydz = oo. 

In particular, the popular harmonic plane wave 

u = U cos(Lvt ± k z .z), p = P sin(ujt ± k z .z), k 2 = uj 2 , U D = const, P D = const, 
clearly occupies the whole 3-space, carries infinite energy 

E = j (U 2 + P 2 )dxdydz = oo 

and, therefore, could hardly be a model of a really created field, i.e. of really existing field objects. 

We recall also that, according to Poisson's theorem for the D'Alembert wave equation (which 
equation is necessarily satisfied by every component of E and B in the pure field case), a spatially 
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finite and smooth enough initial field configuration is strongly time- unstable [1]: the initial condition 
blows up radially and goes to infinity. Hence, Maxwell's equations cannot describe spatially finite 
and time-stable time- dependent free field configurations with soliton-like behavior. The contradictions 
between theory and experiment that became clear at the end of the 19th century were a challenge to 
theoretical physics. Planck and Einstein created the notion of elementary field quanta, named later 
by Lewis [2] the photon. The concept of photon proved to be adequate enough and very seminal, 
and has been widely used in the 20th century physics. However, even now, after a century, we still 
do not have a complete and satisfactory self-consistent theory of single (or individual) photons. It is 
worth recalling at this point the Einstein's remarks of dissatisfaction concerning the linear character 
of Maxwell theory which makes it not able to describe the microstructure of radiation [3] . Along this 
line we may also note here some other results, opinions and nonlinearizations [4]- [14]. 

According to the non-relativistic formulation of Classical Electrodynamics (CED) the electromag- 
netic field has two aspects: electric and magnetic. These two aspects of the field are described by two 
vector fields, or by the corresponding through the Euclidean metric 1-forms, on M 3 : the electric field 
E and the magnetic field B, and a parametric dependence of E and B on time is admitted. Maxwell's 
equations read (with j - the electric current and in Gauss units) 

„ 19E 4f. , ^ 
rotB--— = — j, divE = 47rp, (1) 



c 



The pure field equations are 



1 <9B 

rotE + -— = 0, divB = 0. (2) 

c at 

1 <9E 

rotB ^- = 0, divE = 0, (3) 

c at 

1 dB 

rotE + -— = 0, divB = 0. (4) 

c at 

The first of equations (2) is meant to represent mathematically the results of electromagnetic induction 
experiments. But, it does NOT seem to have done it quite adequately. In fact, as we mentioned earlier, 
these experiments and studies show that field energy-momentum is transformed to mechanical energy- 
momentum, so that the corresponding material bodies (magnet needles, conductors, etc.) change 
their state of motion. The vector equation (2), which is assumed to hold inside all media, including 
the vacuum, manifests different feature, namely, it gives NO information about any field energy- 
momentum transformation to mechanical one in this way: the time change of B is compensated by 
the spatial nonhomogenity of E, in other words, this is an intrinsic dynamic property of the field 
and it has nothing to do with any direct energy-momentum exchange with other physical systems. 
This is quite clearly seen in the relativistic formulation of equations (2) and the divergence of the 
stress-energy-momentum tensor (see further eqns. (8)-(9)). 

This interpretation of the electromagnetic induction experiments has had a very serious impact on 
the further development of field theory. The dual vector equation (3), introduced by Maxwell through 
the so called "displacement current", completes the time-dependent differential description of the field 
inside the frame of CED. The two additional equations divB = 0, divE = impose restrictions on 
the admissible spatial configurations (initial conditions, spatial shape) in the pure field case. The 
usual motivation for these last two equations is that there are no magnetic charges, and NO electric 
charges are present inside the region under consideration. This interpretation is also much stronger 
than needed: a vector field may have non-zero divergence even if it has no singularities, i.e. points 
where it is not defined, e.g. the vector field (sech(x), 0, 0) has NO singularities and is not divergence- 
free. The approximation for continuously distributed electric charge, i.e. div E = Airp, does not justify 
the divergence-free case p = 0, on the contrary, it is consistent with the latter, since NO non-zero 
divergences are allowed outside of the continuous charge distribution. 
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All these clearly unjustified by the experiment and additionally tacitly imposed requirements lead, 
in our opinion, to an inadequate mathematical formulation of the observational facts: as we mentioned 
above, a mathematical consequence of (3)- (4) is that every component of the field necessarily satisfies 
the wave D'Alembert equation = 0, which is highly undesirable in view of the required by this 
equation " blow up" of the physically sensible initial conditions (see the comments in the next section) . 

Another feature of the vacuum Maxwell system (3)- (4), which brings some feeling of dissatisfaction 
(especially to the mathematically inclined physicists) is the following. Let's consider E and B as 1- 
forms on the manifold M 4 with coordinates (x, y, z,£ = ct). It seems natural to expect that every 
solution to (3)-(4) would define completely integrable Pfaff system (E, B). For example, the plane 
wave solution given above defines such a completely integrable Pfaff system, which is easy to check 
through verifying the relations 

dE A E A B = 0, dB A E A B = 0, 

with the plane wave E and B. Now, if we change slightly this plane wave E to E' = E + const.dz, 
const / 0, and keep the same plane wave B, the new couple (E',B) will give again vacuum solution, 
but will no more define completely integrable Pfaff system, since we obtain 

dE' AE'AB = -^const.(u 2 + p 2 ) (; dx A dy A dz A d£, 

dB A E' A B = const. (pu^ — up^)dx A dy A dz A e££. 

These expressions are different from zero in the general case, i.e. for arbitrary functions u and p. 

As for the charge-present case, i.e. equations (l)-(2), some dissatisfaction of a quite different nature 
arises. In fact, every mathematical equation A = B means that A and B are different notations for the 
same element, in other words, it means that on the two sides of "=" stays the same mathematical 
object, but expressed in two different forms. So, every physically meaningful equation necessarily 
implies that the physical nature of A and B must be the same. For example, the Newton's law p = F 
equalizes the same change of momentum expressed in two ways: through the characteristics of the 
particle, and through the characteristics of the external field, and this energy-momentum exchange 
does not destroy neither the particle nor the field. Now we ask: which physical quantity, that is well 
defined as a local characteristic of the field F^u, as well as a local characteristic of the charged particles, 
is represented simultaneously through the field's derivatives as SF^, and through the characteristics of 
the charged particles as j M ? All known experiments show that the field F^ u , considered as a physical 
object, can NOT carry electric charge, only mass objects may be carriers of electric charge. So, from 
this (physical) point of view SF^ should NEVER be made equal to ju, these two quantities have 
different physical nature and can NOT be equivalent. Equations (1) may be considered as acceptable 
only if they formally, i.e. without paying attention to their physical nature, present some sufficient 
condition for another properly from physical point of view justified relation to hold. For example, the 
relativistic version of (1) is given by 8F^ = 4-7rj M , and it is sufficient for the relation F fJil ,5F u = AnF^ v j u 
to hold, but both sides of this last equation have the same physical nature (namely, energy momentum 
change), therefore, it seems better justified and more reliable. Note, however, that an admission of 
this last equation means that the very concept of electromagnetic field, considered as independent 
physical object, has been DRASTICALLY changed (extended), e.g. it propagates NO MORE with 
the velocity of light which is a Lorentz-invariant property. In our view, in such a case, any solution F 
represents just one aspect, the electromagnetic one, of a very complicated physical system and should 
not be considered as electromagnetic field in the proper sense, i.e. as a physical object. 

The important observation at such a situation is that the adequacy of Maxwell description is 
established experimentally NOT directly through the field equations, BUT indirectly, mainly through 
the corresponding Lorentz force, i.e. through the local energy-momentum consevation law. Therefore, 
it seems naturally to pass to a direct energy-momentum description, i.e to extend Maxwell equations 
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to energy-momentum exchange equations, and in this way to overcome the above mentioned physically 
qualitative difference between SFn and j^. 

In this paper we present a concise review of the newly developed Extended Electrodynamics (EED) 
[15]- [20], which was built in trying to find such a reasonable nonlinear extension of Maxwell equations 
together with paying the corresponding respect to the soliton-like view on the structure of the 
electromagnetic field. In the free field case EED follows the rules: 

1. NO new objects should be introduced in the new equations. 

2. The Maxwell energy-momentum quantities should be kept the same, and the local conservation 
relations together with their symmetry (see further relation (7)) should hold. 

So, the first main goal of EED is, of course, to give a consistent field description of single time-stable 
and spatially finite free-field configurations that show a consistent translational-rotational propagation 
as their intrinsic property. As for the electromagnetic fields interacting continuously with some clas- 
sical medium, i.e. a medium carrying continuous bounded electric charges and magnetic moments, 
and having abilities to polarize and magnetize when subject to external electromagnetic field, as we 
mentioned above, our opinion is that the very concept of electromagnetic field in such cases is quite 
different. The observed macro-phenomena concerning properties and propagation may have very little 
to do with those of the free fields. We are not even sure that such fields may be considered as physical 
objects (in the sense considered in Sec. 1.1 to this paper), and to expect appearance of some kind of 
proper dynamics. The very complicated nature of such systems seems hardly to be fully described at 
elementary level, where at every space-time micro-region so many "annihilation-creation" elementary 
acts are performed, the exchanged energy-momentum acquires mechanical degrees of freedom, the 
entropy change becomes significant, some reorganization of the medium may take place and new, not 
electromagnetic, processes may start off, etc. So, the nature of the field-object that survives during 
these continuous "annihilation-creation" events is, in our view, very much different from the one we 
know as electromagnetic free field; the corresponding solution F may be considered just as a seriously 
transformed image of the latter, and it is just one component of a many-component system, and, in 
some cases, even not the most important one. Shortly speaking, the free field solutions are meant to 
represent physical objects, while the other solutions represent just some aspects of a much more com- 
plicated physical system. In our view, only some macro-picture of what really happens at micro-level 
we should try to obtain through solving the corresponding equations. And this picture could be the 
closer to reality, the more adequate to reality are our free field notions and the energy-momentum 
balance equations. That's why the emphasis in our approach is on the internal energy-momentum 
redistribution in the free field case. 

As for the 4-potential approach to Maxwell equations we recall that in the vacuum case the two 
equations dF = 0, d*F = 0, considered on topologically non-trival regions, admit solutions having NO 
4-potentials. For example, the field of a point source is given by (in spherical coordinates originating 
at the source) F = (q/r 2 )dr A d(ct),*F = qsinOdO A dtp and is defined on S 2 x R 2 . Now the new 
couple T = aF + b* F, *J r = —bF + a* F with a, b = const, is a solution admitting NO 4-potential 
on S 2 x I 2 . This remark suggests that the 4-potential does not seem to be a fundamental quantity. 

Turning back to the free field case we note that, in accordance with the soliton-like view, which 
we more or less assume and shall try to follow, the components of the field at every moment t have 
to be represented by smooth, i.e. nonsingular, functions, different from zero, or concentrated, only 
inside a finite 3-dimensional region C R 3 , and this must be achieved by the nonlinearization as- 
sumed. Moreover, we have to take care for the various spatial configurations to be admissible, so, it 
seems desirable the nonlinearization to admit time-stable solutions with arbitrary initial spatial shape. 
Finally, the observed space-time periodical nature of the electromagnetic radiation (the Planck's for- 
mula W = hu in case of photons) requires the new nonlinear equations to be able to describe such a 
periodical nature together with the corresponding frequency dependence during propagation. 

The above remarks go along with Einstein's view that "the whole theory must be based on partial 
differential equations and their singularity-free solutions" [21]. 
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2 From Classical Electrodynamics to Extended Electrodynamics 



The above mentioned idea, together with the assumption for a different interpretation of the concept 
of field: first, considered as independent physical object (field in vacuum), and second, considered 
in presence of external charged matter, is strongly supported by the fact that Maxwell's equations 
together with the Lorentz force give very good and widely accepted expressions for the energy density 
w, the momentum density (the Poynting vector) s, and for the intra-field energy-momentum exchange 
during propagation - the Poynting equation: 

W= i (E2+B2) ' ~ S= ^ BXB ^ ^ = - diVS ~ < 5 > 

Let's consider now the D'Alembert wave equarion □<& = 0. A solution of this equation, meant to 
describe a finite field configuration with arbitrary spatial structure and moving as a whole along some 
spatial direction, which we choose for z-coodinate, would look like as $>(x,y, z + ect),s = ±1, where 
$ is a bounded finite with respect to (x, y, z) function. Substituting this in the equation, clearly, 
the second derivatives of $ with respect to z and t will cancel each other, so, the equation requires $ 
to be a harmonic function with respect to (x, y) : & xx + & yy = 0. According to the Liouville theorem 
in the theory of harmonic functions on M 2 if a harmonic on the whole ]R 2 function is bounded it must 
be constant. Therefore, $ shall not depend on (x,y), i.e. <E> = &(z + ect). 

Considering now the Poynting equation we'll find that it really admits spatially finite solutions 
with soliton-like propagation along the direction chosen. In fact, it is immediately verified that if in 
the above given plane wave expressions for E and B we allow arbitrary dependence of u and p on the 
three spatial coordinates (x, y, z), i.e. u = u(x, y, £ + ez), p = p(x, y, £ + ez), then the corresponding 
E and B will satisfy the Poynting equation and will NOT satisfy Maxwell's vacuum equations. So, 
the Poynting equation is able to describe collections of co-moving photons, which may be considered 
classically as spatially finite radiation pulses with photon-like behavior. Such finite pulses do not "blow 
up" radially, so, this suggests to look for new dynamical equations which also must be consistent with 
the Maxwell local conservation quantities and laws. Hence, if the equations we look for should be 
consistent with the photon structure of electromagnetic radiation, i.e. with the view that any EM- 
pulse consists of a collection of real spatially finite objects, and each one propayates as a whole in 
accordance with the eneryy-momentum that it has been endowed with in the process of its creation, we 
have to find new field equations. 

We begin now the process of extension of Maxwell equations. First, we specially note the well 
known invariance, usually called electric-maynetic duality, of the pure field equations (3)- (4) with 
respect to the transformation 

(E,B)-(-B,E): (6) 

the first couple (3) is transformed into the second one (4), and vice versa. Hence, instead of assuming 
equations (4) we could require invariance of (3) with respect to the transformation (6). Moreover, the 
energy density w, the Poynting vector s, as well as the Poynting relation are also invariant with respect 
to this transformation. This symmetry property of Maxwell equations we consider as a structure one 
and we'll make a serious use of it. 

The above symmetry transformation (6) is extended to the more general duality transformation 



(E,B) - (E',B') = (E,B).a(a,6) = (E,B) 



a b 
-b a 



(aE — 6B, 6E + aB), a = const, b = const. 

(7) 

Note tnat we assume a "right action" of the matrix a on the couple (E, B). The couple (E', B') gives 
again a solution to Maxwell vacuum equations, the new energy density and Poynting vector are equal 
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to the old ones multiplied by (a 2 + b 2 ), and the Poynting equation is satisfied of course. Hence, the 
space of all vacuum solutions factors over the action of the group of matrices of the kind 



a(a, b) 



a b 
—b a 



(a 2 + b 2 ) + 0. 



All such matrices with nonzero determinant form a group with respect to the usual matrix product 
and further this group will be denoted by G. 

This invariance is very important from the point of view of passing to equations having direct 
energy- momentum exchange sense. In fact, in the vacuum case the Lorentz force 



F = pE + -(jxB) 



1 

47T 



E divE + rot B 



1 5E\ n 



should be equal to zero for every solution, so, the transformed (E', B') must give also zero Lorentz 
force for every couple (a, b). This observation will bring us later to the new vacuum equations. 

In the relativistic formulation of CED the difference between the electric and magnetic components 
of the field is already quite conditional, and from the invariant-theoretical point of view there is 
no any difference. However, the 2-aspect (or the 2-component) character of the field is kept in a 
new sense and manifests itself at a different level. To show this we go to Minkowski space-time M, 
spanned by the standard coordinates x M = (x, y, z, £ = ct) with the pseudometric r\ having components 
i]^ = (—1, —1, —1, 1), and r/^ = for /j, ^ v in these coordinates. We shall use also the Hodge star 
operator, defined by 

a A *j3 = —r](a, (3)^J\det{q^\dx Ady Adz A d£, 

where a and (3 are differential forms on M. With this definition of * we obtain for a 2-form F in 
canonical coordinates {*F)^ V = — \e iilUa fiF a '^ . The exterior derivative d and the Hodge * combine 
to give the coderivative (or divergence) operator 5 = (— l) p * _1 d*, which in our case reduces to 
$ = *d*, and on 2-forms we have (SF)^ = — V U F U ^. Note that, reduced on 2-forms, the Hodge 
♦-operator induces a complex structure: ** 2 = —id\2^^, and the same property has transformation 
(6). Explicitly, F and *F look as follows: 



(*F) 
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Fu = B3, F13 = — B2, F23 = Bi, Fu = Ei, F24 = E2, -F34 = E3. 
E 3 , (*F) 13 = -E 2 , (*F) 23 = Ei, (*F) U = -Bi, (*F) 24 = -B 2 , (*F) 34 



B, 



So, if the non-relativistic vector couple (E, B) corresponds to the relativistic 2-form F, then the 
transformed by transformation (6) couple (— B,E) obviously corresponds to the relativistic 2-form 
*F. Hence, the Hodge * is the relativistic image of transformation (6). 

Recall now the relativistic energy-momentum tensor, which can be represented in the following 
two equivalent forms (see the end of this section): 



Qu 



1 

47T 



1 

8^ 



-F^F™ -(*F)^(*F) L 



(8) 



It is quite clearly seen, that F and *F participate in the same way in Q v , so Q 1 ^ is invariant with 
respect to F — ► *F and the full energy-momentum densities of the field are obtained through summing 
up the energy-momentum densities, carried by F and *F. That is how the above mentioned invariance 
of the energy density and Poynting vector look like in the relativistic formalism. Now, the Poynting 
relation corresponds to the zero-divergence of Q v ^. 



1 



{SF) V + {*F) IW {S*F) V 



0. 



(9) 
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The obvious invariance of the Poynting relation with respect to (6) corresponds here to the invariance 
of (9) with respect to the Hodge *. Moreover, the above relation (9) obviously suggests that the field 
is potentially able to exchange energy-momentum through F, as well as through *F independently, 
provided appropriate external field-object is present. 

Now, the idea of the generalization CED^EED could be formulated as follows. The above men- 
tioned momentum balance interpretation of the second Newton's law p = F says: the momentum 
gained by the particle is lost by the external field, this is the true sense of this relation. If there is no 
external field, then F = 0, and giving the explicit form of the particle momentum p = mv = mr we 
obtain the equations of motion p = mv = of a free particle, implying rh = 0. The same argument 
should work with respect to the field: if there is no particles then m = 0, p = and p = 0, so, we 
must obtain the equations of motion (at least one of them) for the free field: F = 0, after we give 
the explicit form of the dependence of F on the field functions and their derivatives. In other words, 
we have to determine in terms of the field functions and their derivatives how much 
energy-momentum the field is potentially able to transfer to another physical system in 
case of the presence of the latter and its ability to interact with the field, i.e. to accept 
the corresponding energy-momentum. 

Following this idea in non-relativistic terms we have to compute the expression 

d 

making use of Maxwell equations (l)-(2), and to put the result equal to zero. In the free field case 
jfcPmech = 0. The corresponding relation may be found in the textbooks (for example see [22], ch.6, 
Paragraph 9) and it reads 

d E\ / B \ 

rotB-— J x B + E div E + (rotE + — J xE + BdivB = 0. (10) 

In the frame of Maxwell theory the first two terms on the left of (10) give the Lorentz force in field 
terms, so, in the free field case we must have 

rotB - — J x B + EdivE = 0. (11) 

Hence, it follows now from (10) and (11) that we must also have 

<9B \ 

rot E + — J x E + B div B = 0. (12) 

If we forget now about the sourceless Maxwell equations (3)- (4) we may say that these two nonlinear 
vector equations (11)-(12) describe locally some intrinsic energy-momentum redistribution during the 
propagation. 

In order to have the complete local picture of the intrinsic energy-momentum redistribution we 
must take also in view terms of the following form: 

1 (8E dB\ 

- — xE + Bx— , Ex rotB, B x rotE, EdivB, BdivE. 

c \ot at J 

A reliable approach to obtain the right expression seems to be the same invariance with respect to the 
transformation (7) to hold for (11)-(12). This gives the additional vector equation: 

<9B \ / <9E \ 

n .1 E +- — J x B - E div B + (rot B - — J x E - B div E = 0. (13) 
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The above consideration may be put in the following way. According to our approach, the more 
important aspects of the field from the point of view to write its proper dynamical equations are its 
abilities to exchange energy-momentum with the rest of the world: after finding out the corresponding 
quantities which describe quantitatively these abilities we put them equal to zero. 

The important moment in such an approach is to have an adequate notion of the mathematical 
nature of the field. The dual symmetry of Maxwell equations suggests that in the nonrelativistic 
formalism the mathematical nature of the field is 1-form uj = E (g> e\ + B (g> e 2 on M 3 with values in 
an appropriate 2-dimensional vector space Q with the basis (ei, e 2 ), where the group G acts as linear 
transformations. In fact, if a(a,b) G G, then consider the new basis {e' l} e' 2 ) given by 



a 2 + b 2 



(aei-be 2 ), e 2 



1 



a 2 + b 2 



(bei + ae 2 ). 



Accordingly, a € G transforms the basis through right action by means of (a x )* 
the "new" solution u/(E',B') is, in fact, the "old" solution w(E,B): 



a/det(a). Then 



J = E' ® ei + B' ® e 2 = (aE - 6B) °^ - ^ 2 + (6E + aB) ® ^ + - E < , +- B , 2 - 



a 2 + 6 2 



i.e., the "new" solution a/(E',B'), represented in the new basis (e[,e 2 ) coincides with the "old" 
solution w(E, B), represented in the old basis (ei,e2). 

In view of this we may consider transformations (7) as nonessential, i.e. we may consider (E, B) 
and (E',B') as two different representations in corresponding bases of Q of the same solution w, and 
the two bases are connected with appropriate element of G. 

Such an interpretation is approporiate and useful if the field shows some invariant properties with 
respect to this class of transformations. We shall show further that the nonlinear solutions of (11)-(13) 
do have such invariant properties. For example, if the field has zero invariants (the so called "null 
field"): 

h = (B 2 - E 2 ) = 0, I 2 = 2E.B = 0, 

then all transformations (7) keep unchanged these zero- values of I\ and I 2 . In fact, under the trans- 
formation (7) the two invariants transform in the following way: 

l[ = (a 2 - b 2 ) h + 2a6/ 2 , I 2 = -2abh + (a 2 - b 2 ) I 2 , 

and the determinant of this transformation is (a 2 + 6 2 ) 2 / 0. So, a null field stays a null field under 
the dual transformations (7). Moreover, NO non-null field can be transformed to a null field by means 
of transformations (7), and, conversely, NO null field can be transformed to a non-null field in this 
way. Further we shall show that the "null fields" admit also other G- invariant properties being closely 
connected with availability of rotational component of propagation. 

From this G-covariant point of view, the well known Lorentz force appears to be just one component 
of the corresponding covariant with respect to transformations (7), or to the group G, G-covariant 
Lorentz force. So, in order to obtain this G-covariant Lorentz force we have to replace E and B in the 
expression 

rotB-— J xB + EdivE, 

by (E',B') given by (7). In the vacuum case, i.e. when the field keeps its integral energy-momentum 
unchanged, the G-covariant Lorentz force F obtained must be zero for any values (a, b) of the trans- 
formation parameters. 

After the corresponding computation we obtain 



<9E , 

rotB - — ) x B + E div E 



+ b 2 



dB , 

rot E + — ) x E + B div B 



+ 
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+ab 



dB 



rot E + — x B - E div B + rot B - — x E - B div E 



dB 



0. 



Note that the G — covariant Lorentz force F has three vector components given by the quantities in 
the brackets. Since the constants (a, b) are arbitrary the equations (11)-(13) follow. 

Another formal way to come to equations (11)-(13) is the following. Consider a 2-dimensional real 
vector space Q with basis (ei,e2), and the objects 



dB 



dE 



rot E + — — <g ei + rot B - — - <g e 2 , wi = E®ei+B®e 2 



uj 2 = -B (g ei + E (g e 2 , 6* = -(divB) (g> ei + (divE) e 2 . 

Now define a new object by means of the bilinear maps "vector product x", "symmetrized tensor 
product V" and usual product of vector fields and functions (denoted by ".") as follows: 



(x, V) 



<9B\ 



dE 



rot E + — — <g ei + rot B - — - <g e 2 , E e 1 + B (g) e 2 



9£ 



+ 



(., V) [-B(gei +E <g> e 2 , -(divB) (g e x + (divE) <g e 2 ] 



9 B 

rot E + -— ) x E + Bdiv B 



dB 



ei V ei + 



9 E 

rotB - — ] x B f- EdivE 



rot E + — x B - E div B + rot B - — x E - B div E 



dE 



d£ 



e 2 V e 2 + 
ei v e 2 = F. 



Now, putting this product equal to zero, we obtain the required equations. 

The left hand sides of equations (11)-(13), i.e. the components of the 3-dimensional (nonrelativistic) 
G-covariant Lorentz force F, define three vector quantities of energy-momentum change, which appear 
to be the field's tools controlling the energy-momentum exchange process during interaction. We can 
also say that these three terms define locally the energy-momentum which the field is potentially able 
to give to some other physical system, so, the field demonstrates three different abilities of possible 
energy-momentum exchange. 

Obviously, equations (11)-(13) contain all Maxwell vacuum solutions, and they also give new, 
non-Maxwellean nonlinear solutions satisfying the following inequalities: 



dB " 

™tE + — ] /(). 



dB 

rot E- — I/O, divE/0, divB + 0. 



(14) 



It is obvious that in the nonlinear case we obtain immediately from (11)-(13) the following relations: 



rotB-^|.E 



0, 



dB 

rotE + — | .B = 0, E.B 



0, 



and combining with (13) we obtain additionally (in view of E.B = 0) 



B. rotB 



<9E\ / _ dB\ „ 
- — - E. rotE + — = B. 



d£ 



dB 



d£ 



.rotB - E.rotE = 0. 



(15) 



(16) 



So, in the nonlinear case, from (15) we obtain readily the Poynting relation, and (in hydrodynamic 
terms) we may say that relation (16) requires the local helicities B.rotB and E.rotE of E and B [23,24] 
to be always equal. 

In relativistic terms we first recall Maxwell's equations in presence of charge distribution: 



S * F = 0, SF = 4irj. 



(17) 
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In the pure field case we have 

5 * F = 0, 6F = 0, 

or, in terms of d 

dF = 0, d*F = 0. 

Recall now from relations (8)-(9) the explicit forms of the stress-energy-momentum tensor and its 
divergence. From Maxwell equations (17) it follows that field energy-momentum is transferred to an- 
other physical system (represented through the electric current j) ONLY through the term 4irF llu j u = 
Fnv(8F) v (the relativistic Lorentz force), since the term (*F)^ U (5* F) u is always equal to zero because 
of the first of equations (17) which is meant to represent the Faraday's induction law. 

Following the same line of consideration as in the nonrelativistic case we extend Maxwell's free 
field equations 5F = 0, 5 * F = as follows. In the field expression of the relativistic Lorentz force 
Ffi, v bF v we replace F by aF + b * F, so we obtain the relativistic components of the G — covariant 
Lorentz force F: 

a 2 F^5F» + b 2 {*F)^{5 * Ff + ab[F^{5 * Ff + (*F)^5F U ] . 

In view of the arbitrariness of the parameters (a, b) after putting this expression equal to zero we 
obtain equations (11)-(13) in relativistic notation. The equation 

F^SF" = (18) 

extends the relativistic Maxwell's equation SF = and in standard coordinates (18) gives equations 
(11) and the first of (15). Equations (12) and the second of (15) are presented by the equation 

(*F)^(S*FY = 0. (19) 

Looking back to the divergence relation (9) for the energy- momentum tensor Q 1 ^ we see that (18) 
and (19) require zero values for the two naturally arising distinguished parts of V U Q^, and so the 
whole divergence is zero. This allows Q 1 ^ to be assumed as energy-momentum tensor for the nonlinear 
solutions. 

The relativistic version of (13) and (16) looks as follows: 

F^(5 * F f + (*F)^ U 5F U = 0. (20) 

This relativistic form of the equations (11)-(13) and (15)-(16) goes along with the above inter- 
pretation: equations (18)- (19), together with equation (9), mean that no field energy- momentum is 
transferred to any other physical system through any of the two components F and *F, so, only 
intra-field energy-momentum redistribution should take place; equation (20) characterizes this intra- 
field energy momentum exchange in the following sense: the energy-momentum quantity, transferred 
locally from F to *F and given by F fll/ (S * F) u , is always equal to that transferred locally from *F to 
F which, in turn, is given by (*F) flu (5F) 1/ . 

In terms of the "change-conservation" concept considered in Sec. 1.1 , we could say that the self- 
projections F^ V (8F) U and (*F)^ U (5 * FY of the admissible changes 5F and 5 * F on F and *F, 
respectively, vanish, and the cross-projections (*F) lll/ (5F) u and (F) I1U (5 * F) u have the same absolute 
value. 

So, we can think of the electromagnetic field as a 2-vector component field O, = F <£> e\ + *F <g> e2, 
where (ei, e?) is an appropriately chosen basis of the same real 2-dimensional vector space Q, and later 
the nature of Q will be specially considered. 

In presence of external fields (media), which exchange energy-momentum with our field F, the 
right hand sides of (18)- (20) will not be zero in general. The energy-momentum local quantities that 
flow out of our field through its two components (F, *F) and are given by the left hand sides of (18)- 
(20), in accordance with the local energy-momentum conservation law have to be absorbed by the 
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external field. Hence, these same quantities have to be expressed in terms of the field functions of the 
external field. 

Maxwell equations give in the most general case only one such expression, namely, when the 
external field is represented by continuously distributed free and bounded charged particles. This 
expression reads: 

^TTF^j , j = j f ree + j bound' 

where jf ree = pu v is the standard 4-current for freely moving charged particles, Ji, oun d ls expressed 
through the polarization P and magnetization M vectors of the (dielectric) medium by 

f <9P ] 

Jbound = Obound, Pbound) = < CrotM + — , -divP > . 

So, Maxwell equations allow NO energy- momentum exchanges through *F. The usual justification 
for this is the absence of magnetic charges. In our view this motivation is insufficient and has to be 
left off. 

In the frame of EED, in accordance with the concept of G—covariant Lorentz force, we reject this 
limitation, in general. Energy-momentum exchanges through F, as well as, through *F are allowed. 
Moreover, EED does not forbid some media to influence the energy-momentum transfers between F 
and *F, favoring exchanges with F or *F in correspondence with medium's own structure. Formally, 
this means that the right hand side of (20), in general, may also be different from zero. Hence, in the 
most general case, EED assumes the following equations: 

F^SF" = F^a 1 )", (21) 
(**V0$*F)" = (*f) F (a 4 )", (22) 
F au {5 * F) v + (*F)^5F» = F„ v {a 2 Y + (**V(a 3 )". (23) 

The new objects a l ,i = 1,2,3,4 are four 1-forms, and they represent the abilities of the corre- 
sponding medium for energy-momentum exchange with the field. Clearly, their number corresponds 
to the four different abilities of the field to exchange energy-momentum, which are represented by 
F^dF", (*F) I11/ (5*FY, (*F) UU 5F U and F tlv (8*FY- These four 1-forms shall be expressed by means 
of the external field functions and their derivatives. This is a difficult problem and its appropriate res- 
olution requires a serious knowledge of the external field considered. In most practical cases however 
people are interested only in how much energy-momentum has been transferred, and the transferred 
energy- momentum is given by some bilinear functions on F and a l ,i = 1, 2, 3, 4 as it is in the Maxwell 
case. 

Here is the 3-dimensional form of the above equations (the bold a 1 denotes the spatial part of the 
corresponding a 1 ): 



<9E \ 

rotB - —J x B + EdivE = a 1 x B + E(a x ) 4 , 

E. (rotB - ^) = E.a\ 
rotE + x E + BdivB = a 4 x E - B(a 4 ) 4 , 

B. (rotE + ^)=B.a 4 , 

/ ^ dB\ n ( „ dE 
rotE + — x B + 



^rotE + - x B + |^rotB - - x E - BdivE - EdivB = 

= a 2 x B + E(a 2 ) 4 + a 3 x E - B(a 3 ) 4 , 
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b - ( rotB - W) - e ( rotE + W) = B - a3 - E - a2 - 

As an example we recall the case of availability of magnetic charges with density p m and magnetic 
current j m . Then the corresponding four 1-forms are 



a 



"4tt 




"47T. 


— Je,4vr/) e 




Je, -47T/? e 


c 




C 



2 4 

a = a = 



4-ir. 



-j m , -47r/9„ 



The corresponding three vector equations (we omit the scalar ones) look like 



( 



rotE + 



dB 



x E + BdivB 



47T 



(-j m x E) + 47rp m B, 



(^rotB - ^ x B + EdivE = ^(j e xB)+ 47rp e E, 
/ <9E\ / <9B\ 

(rotB - — J x E + (rotE + — J x B - BdivE - EdivB = 

4vr 

= — (-j m x B + j e x E) - 4^(p m E + Pe B). 

Returning to the general case we note that the important moment is that although the nature 
of the field may significantly change, the interaction, i.e. the energy-momentum exchange, must 
NOT destroy the medium. So, definite integrability properties of the medium MUST be available, 
and these integrability properties should be expressible through the four 1-forms a l ,i = 1,2,3,4. In 
the Maxwell case, where the charged particles represent any medium, this property implicitly presents 
through the implied stability of the charged particles, and it is mathematically represented by the local 
integrability of the electric current vector field j^: the corresponding system of ordinary differential 
equations = has always solution at given initial conditions. 

If the physical system "electromagnetic field+medium" is energy-momentum isolated, i.e. no 
energy-momentum flows out of it, and the system does NOT destroy itself during interaction, EED 
assumes, in addition to equations (21)-(23), the following: 

Every couple (a l ,a?),i / j, defines a completely integrable 2-dimensional Pfaff system. 

This assumption means that the following equations hold: 

da* A a* A a? = 0, i,j = 1,2,3,4. (24) 

Equations (21)-(24) constitute the basic system of equations of EED. Of course, the various special 
cases can be characterized by adding some new consistent with (21)-(24) equations and relations. 

We are going now to express equations (21)-(23) as one relation, making use of the earlier in- 
troduced object fi = F ® e\ + *F ® e2 and combining the 1-forms a l ,i = 1,2,3,4 in two ^-valued 
1-forms: 

$ = a 1 <g> ei + a 2 (g> e 2 , * = a 3 <S> e\ + a 4 ® e 2 . 
The basis (ei,e2) defines two projections tt\ and 1T2: 

TTlQ = F <S> ei, 7T 2 r2 = *F <S> B2- 

Every bilinear map 93 : Q x Q — > W, where W is some linear space, defines corresponding product in 
the (/-valued differential forms by means of the relation 



(p(Q\ ei, Q 3 2 (g> ej) = Q\ Al^® tp(ei,ej). 
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Now, let ip = V, where "V" is the symmetrized tensor product. Recalling * D * 2 = —id we obtain 

V($, *7TiSl) + V(*, *tt 2 ^) = a 1 A*F®eiVei-Q 4 AF®e 2 Ve 2 + (-a 3 A F + a 2 A *F) <g> e x V e 2 . 

Now, equations (21)-(23) are equivalent to 

V(<50, = V($, *7TiO) + V(*, *7T 2 S1). (25) 

The Maxwell case (with zero magnetic charges) corresponds to a 2 = a 4 = 0, a 1 = a 3 = 4n(5S + j), 
where the 2-form S is defined by the polarization 3- vector P and the magnetization 3-vector M in the 
similar way as F is defined by (E,B). Explicitly, 7r 2 <5r2 = 0, ttiSQ = 4ir(5S + j) <8> e\. 
In vacuum (25) reduces to 

F = V (SSI, *Q) = F flu 5F u dx^ O ei V ei + (*F)^(S * Ffdx^ 1 <g> e 2 V e 2 + 



+ 



(5 * FY + (*F)^5F U 



dx^ ® ei V e 2 = 0, 



which says that the relativistic G-covariant Lorentz force, denoted also by F, is equal to zero. 

As for the energy- momentum tensor Q^ v of the vacuum solutions, considered as a symmetric 
2-form on M, it is defined in terms of Q, as follows: 

Q(X,Y) = ±*g[i(X)n,*i(Y)ti\, 

where (X, Y) are two arbitrary vector fields on M, g is the metric in Q defined by g(a, (5) = ^tr(a./3*), 
and (3* is the transposed to (3. 

Further we are going to consider the pure field (vacuum) case a 1 = fully. As for the non- 
vacuum case, we shall show explicitly some (3+l)-soliton solutions with well defined integral conserved 
quantities. We note the following useful relations in Minkowski space-time. Let a be 1-form, F, G be 
two 2-forms and H be a 3-form. Then we have the identities: 

* (a A *F) = -a a F^dx^; (26) 

* (F A *H) = X -F™H CVil dx»- (27) 

\F^G a Hl = F^G™ - (*G)^(*F) UCT - (28) 

Making use of the identities (26)-(28) equations (21)-(23) can be represented equivalently and 
respectively as follows: 

(*F) ffT {d*F) <rriJl = F IJa/ (a 1 ) v , a<r (29) 
F^(dF) aT)l = (*F)^(a 4 y, a<r (30) 
-(*Fy T (dF) aTIM - F° T (d * F) CTTM = F^{a 2 Y + (*F) /1I/ (a 3 )" , a < r. (31) 

Finally we note the following. The identity (28) directly leads to relation (8). In fact, we put in 
(28) G = F, after that we multiply by ^, now we represent the so obtained first term on the right as 
F va F va - \F^F va , and finally we transfer F im F va to the left side. 

Another suggestion comes if we put in (28) G = *F. Then multiplying by ^, taking in view that 
* * F = —F, and transferring everything to the left we obtain the identity 

\F a(3 (*F)^5^ - F, a (*Fr = 0. 

This last relation suggests that NO interaction energy-momentum exists between F and *F, which is in 
correspondence with our equation (20) which states that even if some energy-momentum is transferred 
locally from F to *F, the same quantity of energy-momentum is simultaneously locally transferred 
from *F to F. We note also that these last remarks are of pure algebraic nature and they do NOT 
depend on whether F satisfies or does not satisfy any equations or additional conditions. 
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3 EED: General properties of the vacuum non-linear solutions 

First we note the following coordinate free form of equations (18)-(20) in terms of 5 and d: 

(*F) A *d * F = 5F A *F = 0, (32) 
F A *dF = —F A 5 * F = 0, (33) 
F A *d * F + (*F) A *dF = 5F A F - 5 * F A *F = 0, (34) 

We begin with establishing some elementary symmetry properties of these equations. 

Proposition 1. Equations (32)-(34) are invariant with respect to a conformal change of the metric. 

Proof. Let rj — > g = f 2 rj, where f 2 (m) > 0, m G M. It is seen from the coordinate-free d-form of 
the vacuum equations that only the restrictions *2 and *3 of the Hodge * participate, moreover, * 3 is 
not differentiated. Since * 2 is conformally invariant, the general conformal invariance of the equations 
will depend on *3. But, for this restriction of * we have (* g )3 = / 2 (*r?)3 , so the left-hand sides of the 
equations are multiplied by / / of a given degree, which does not change the equations. 

Proposition 2. Equations (29)-(31) are invariant with respect to the transformation: 
F -► *F, a 1 -► a 4 , a 2 -»• (-a 3 ), a 3 -► a 2 , a 4 -► (-a 1 ). 

Proof. Obvious. 

Proposition 3. The vacuum equations (32)-(34) are invariant with respect to the transformation 

F^F = aF-b*F, *F -> *F = bF + a * F, a, b e M. 

Proof. We substitute and obtain: 

5FA*F = a 2 (5F A *F) - b 2 (5 * F A F) + ab{5F A F - 5 * F A *F) 

5 * T A T = a 2 (5 * F A F) — b 2 (5F A *F) + ab{5F A F - 5 * F A *F) 

5F A T - 5 * T A *F = (a 2 - b 2 ){5F A F - 5 * F A *F) - 2ab(5F A *F + 5 * F A F). 

It is seen that if F defines a solution then T also defines a solution. Conversely, if T defines a solution 
then subtracting the second equation from the first and taking in view that a 2 + b 2 / we obtain 
SF A *F = (5 * F) A F. So, if 5F A *F = 5*FAF = from the first two equations we obtain 
SF A F - S * F A *F = -(SF A F - S * F A *F) = 0. The proposition follows. 

We recall the following property of vectors in Minkowski space-time (further referred to as BP1): 
Basic property 1.: There are NO mutually orthogonal time-like vectors in Minkowski space-time. 
We recall also [27] the following relation between the eigen properties of F^ and Q 1 ^. 

Basic property 2: All eigen vectors of F and *F are eigen vectors of Q too. 
(Further referred to as BP2). 

It is quite clear that the solutions of our non- linear equations are naturally divided into three subclasses: 

1. Linear (Maxwellean) , i.e. those, satisfying 5F = 0, 5 * F = 0. This subclass will not be of 
interest since it is well known. 

2. Semilinear, i.e. those, satisfying S * F = 0, 8F ^ 0. (or 5 * F ^ 0, 8F = 0). This subclass of 
solutions, as it will soon become clear, does NOT admit rotational component of propagation, so it is 
appropriate to describe just running wave kind behavior. 
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3. Nonlinear, i.e. those, satisfying SF ^ 0, 5 * F / 0. 



Further we shall consider only the nonlinear solutions, and the semilinear solutions will be char- 
acterized as having the additional property 5 * F = 0, (or 5F = 0). 

The following Proposition is of crucial importance for all further studies of the nonlinear solutions. 

Proposition 4. All nonlinear solutions have zero invariants: 

h = \F^F^ = ±^Jdet[(F±*F)^] = 0, 

h = \{*F)„ V F^ = ±2^/det(F^) = 0. 

Proof. Recall the field equations in the form: 

F^SFf = 0, {*F)^(5 * F) v = 0, F^(S * FY + (*F)^(SFY = 0. 

It is clearly seen that the first two groups of these equations may be considered as two linear homoge- 
neous systems with respect to 5F^ and 5*F^ respectively. These homogeneous systems have non-zero 
solutions, which is possible only if det(F flu ) = det((*F) MJy ) = 0, i.e. if I2 = 2E.B = 0. Further, 
summing up these three systems of equations, we obtain 



(F + *F)^(5F + 5*FY = 0. 



If now (SF + 5 * FY + 0, then 
= det||(F + *F)^| 



1 2 



-(F + *F) flu (*F-FY u 



1/7 /.; ; " 



1 2 



If 5F V = -(5 * FY / 0, we sum up the first two systems and obtain (*F - F) )IU (5 * FY = 0. 
Consequently, 



= det||(*F-F) 



flU I 



I/? 

r, V 



This completes the proof. 

Corollary. All nonlinear solutions are null fields, so (*F) jJi(T F l " T = and F^ a F ua = (*F)^ a (*FY a - 
Corollary. The vector 5F^ is an eigen vector of the vector (5 * FY is an eigen vector of 

(*FY- 

Corollary. The vectors 5F^ and (S * FY are eigen vectors of the energy tensor Q^. 

Basic Property 3. [27] In the null field case Q v ^ has just one isotropic eigen direction, defined by 
the isotropic vector (, and all of its other eigen directions are space-like. (Further referred to as BP3). 



Proposition 5. All nonlinear solutions satisfy the conditions 

(SF)„(S * FY = 0, \SF\ = \5*F\ (35) 
Proof. We form the inner product i(5 * F)(SF A *F) = and get 

(6 * FY(SF)^(*F) - 5F A (5 * FY(*F) flu dx u = 0. 

Because of the obvious nullification of the second term the first term will be equal to zero (at non-zero 
*F) only if (SF)^ * FY = 0. 

Further we form the inner product i(5 * F)(SF A F — 5 * F A *F) = and obtain 
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(5 * F f{5F)^F - 5F A (6 * F^F^dx"- 
-(6 * F) 2 {*F) + 5 * F f\{5 * Ff(*F)^dx u = 0. 
Clearly, the first and the last terms are equal to zero. So, the inner product by 6F gives 

(5F) 2 (5 * FfF^dx v - [{SFYiS * F) v 'F^\ 5F + (S * F) 2 (8FY(*F)^dx» = 0. 

The second term of this equality is zero. Besides, (5 * F)^F flu dx u = — (5FY(*F)^ u dx u . So, 

[{5Ff - (5 * F) 2 ] (SFr(*F)^dx» = 0. 

Now, if (5F)> J -(*F)i XU dx u / 0, then the relation \SF\ = \8 * F\ follows immediately. If 
(5FY(*F) fMl/ dx L/ = = -(5 * FYF^ v dx v according to the third equation of (22), we shall show 
that (5F) 2 = (5* F) 2 = 0. In fact, forming the inner product i(5F){5F A *F) = , we get 

(5F) 2 * F - 5F A {5FY{*F)^dx v = (5F) 2 * F = 0. 

In a similar way, forming the inner product i(5 * F)(6 * F A F) = we have 

(5 * F) 2 F - 5{*F) A (5 * FYF^ v dx v = {5* F) 2 F = 0. 

This completes the proof. It follows from this Proposition and from BP1 that 

1. SF and 5 * F can NOT be timedike, 

2. SF and 5 * F are simultaneously space-like, i.e. (5F) 2 = (5 * F) 2 < 0, or simultaneously 
isotropic, i.e. \SF\ = \5 * F\ = 0. We note that in this last case the isotropic vectors 5F and 5 * F 
are also eigen vectors of Q^, and since Q u ^ has just one isotropic eigen direction, which we denoted by 
C, we conclude that 5F, 5 * F and £ are collinear. 

In our further study of the nonlinear solutions we shall make use of the following. As it is shown 
in [27] at zero invariants I\ = I2 = the following representation holds: 

F = A A(, *F = A* AC, 

where A and A* are 1-forms, ( is the corresponding to £ through the pseudometric 77 1-form: ^ = 
rj^vCY ■ It follows that 

faC = *fac = o. 

Remark. Further we are going to skip the "hat" over £, and from the context it will be clear the 
meaning of Q: one-form, or vector field. 

We establish now some useful properties of these quantities. 

Proposition 6. The following relations hold: 
l°.A^ = A*^ = 
2°.A°(*F) atl = 0, (A*YF afi = 
3°.A a A* a = 0, A 2 = (A*) 2 < 0. 
Proof. In order to prove 1° we note 

= h = \F,uF^ = \(A,C V - A„^)(A»C ~ 

= hpA^.M" - 2(^o 2 ) = -(A,cn 2 - 

The second of 1° is proved in the same way just replacing F with *F and A with A*. 
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To prove 2° we make use of relation (26). 



= *{A* A A* A C) = *(A* A *F) = -{A*) a F^dx^ 

Similarly 

= *(A A A A C) = *(A A F) = A a (*F) CTI ,dx^. 

Hence, A* is an eigen vector of F and A is an eigen vector of *F. Therefore, A 2 < 0, (A*) 2 < 0, 
i.e. these two vectors (or 1-forms) are space-like. The case \A\ = is not considered since then A is 
collinear to ( and F = A A £ = 0, so *F = too. 
Now, 3° follows from 2° because 

= -(A*)°F a „ = -(A*r(A^ - A,C,) = -(A\A)^ + A^A* .Q = -(A*.A)^. 

Finally we express Q 1 ^ in terms of A, or A*, and (. Since I± = and Q^(F) = Q V J*F) we have 



1 



Hence, since > and C4C 4 > we obtain again that A 2 = (A*) 2 < 0. 

We normalize £, i.e. we divide by £4 / 0, so, we assume further that ^ = (£1, C2, C3> !)■ Now, 
from the local conservation law V V Q^ V = it follows that the isotropic eigen direction defined by £ 
defines geodesic lines: C, U V U C^ = 0. In fact, 



0. 



This relation holds for any ^ = 1,2,3,4. We consider it for \x = 4 and recall that for /j, = 4 we 
have £ 4 = 1. In our coordinates X7 U = d u , and we obtain that the second term becomes zero, so 
V 1/ {A 2 C i v ) = 0. Therefore, ("V^C 1 = 0, which means that all the trajectories of ( are parallel 
straight isotropic lines. Hence, with every nonlinear solution F we are allowed to introduce F- 
adapted coordinate system by the requirement that the trajectories of £ to be parallel to the plane 
(z, £). In such a coordinate system we obtain ^ = (0,0, e, 1), e = ±1. From 3-dimensional point of 
view this means that the field propagates along the coordinate z, and e = — 1 implies propagation 
along z from — 00 to +00, while e = +1 implies propagation from 00 to —00. 

Corollary. The translational direction of propagation of any null-field is determined intrinsically. 

We note that in Maxwell theory, as well as in Born-Infeld theory [9] (the latter reduces to Maxwell 
theory in the null-field case) all time stable null-field solutions are spatially infinite, otherwise they 
have to blow-up radially according to Poisson theorem, i.e. Q shall not be unique. We shall see that 
EED has no problems in this respect, i.e. spatially finite and time stable null field solutions are 
allowed. 

We express now F, *F, A and A* in the corresponding F- adapted coordinate system making use of 
the relations F = A A C, *F = A* A C, where ( = edz + d£. 



F\2 — — 0, F13 — £^14, F23 — sF 24 , 

(*F)ia = (*F) U = 0, (*F)i 3 = e(*F) u = -F 24 , (*F) 23 = e(*F) 24 = F 14 . 
Moreover, from A.( = A*.( = it follows that in an F-adapted coordinate system we obtain 
A = A 4 dx + A 2 dy + f.( = {F u )dx + (F 24 )dy + f.( 

A* = A\dx + A* 2 dy + f*.( = (-F 23 )dx + (F 13 )dy + f*.( = -{eA 2 )dx + (eAJdy + f*.(, 
where / and /* are two arbitrary functions. 



(36) 



(37) 
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Having this in mind we prove the following 
Proposition 6. All nonlinear solutions satisfy the relations: 

^(SF) ti = 0, C»(6*F) li = 0. (38) 

Proof. We form the inner product i(()(SF A *F) = : 

[CiSF)^] * F - SF A (Cn*F)^dx» = 

= [C(SF)^} A* A C - (SF A C)e(A*)» + (SF A A*)C^ = 0. 

Since the second and the third terms are equal to zero and *F ^ 0, then ^(SF)^ = 0. Similarly, from 
the equation (5 * F) A F = we get C^(S * F)„ = 0. The proposition is proved. 

Proposition 7. If F ^ is a nonlinear solution with \SF\ / then SF A F = S * F A *F / 0. 

Proof. From F / we have (5F) 2 / 0. Recall that *(SF A F) = i(SF) * F = (SFy(*F)^ v dx u . 
Assume now the opposite, i.e. that SF A F = 0. We compute 

= i(5F)(5F A *F) = (SF) 2 * F - SF A i(5F) * F = (SF) 2 * F, 

i.e. (SF) 2 = 0, or *F = 0, which contradicts the assumption F / 0. 

Conversely, if 5F A F / then surely 5F / 0. So, we may have \SF\ ^ 0, or |5F| = 0. If |<5F| = 
then SF must be colinear to (, i.e. (AF^O, which is not possible since F = A A (. 

We also note that the 3-form SF A F is isotropic: (SF A F) 2 = 0. In fact 

(SF A F)^ CT (<5F A F)^ = F^ [(SF) 2 F^ U ] - F» v [SF A i(SF)F] ^ = (SF) 2 F^F^ U = 0, fi < u < a. 

Finally we note that there are no nonlinear spherically symmetric solutions, i.e. if F is a spherically 
symmetric solution then it is a solution of Maxwell's equations. In fact, the most general spherically 
symmetric 2-form in spherical coordinates, originating at the symmetry center, is 

F = f(r, £) dr A d£ + h(r, f ) sin6> dO A dip. 

Now the equation F A *dF = requires h(r,^) = const, and the equation SF A *F = requires 
/(r,£) = const /r 2 . It follows: dF = 0, <5F = 0. 

4 EED: Further properties of the nonlinear solutions 

Let's introduce the notations F14 = u, F24 = p, so, in an F-adapted coordinate system we can write 

F = eudx A dz + ejxiy A <iz + udx A d£ + pdy A d£ 
* F = — pdx A + udy A dz — epdx Ad£ + eudy A 

This form of F shows that every nonlinear solution defines a map 

a F : M -> M 2 , ap(x,y,z,£) = [u(x, y, z, £),p(x, y, z, £)]. 

We endow now IR 2 with the canonical complex structure J and identify (IR 2 , J) with the complex 
numbers C. Hence, our nonlinear solution F defines the complex valued function op = u + ip. This 
correspondence and its inverse will be studied further. 

In the F-adapted coordinate system only 4 of the components are different from zero, namely: 

Qi = -Ql = eQl = -eQ\ =\A 2 \= u 2 +p 2 . 



22 



The complex valued function ctp denned above by the nonlinear solution F has module \ap\ = 



\J 'u 2 + p 2 and phase ip = arctg(j>/u) = arccos (u/y^u 2 + p 2 ) . We shall come now to these two quan- 
tities in a coordinate free way. We show first how the nonlinear solution F defines at every point a 
pseudoorthonormal basis in the corresponding tangent and cotangent spaces. The nonzero 1-forms 
A = A - fC and A* = A* - f*( are normed to A = A/\A\ and A* = A* /\A*\. We note that \ A\ = \ A\ 
and | A* \ = \A*\. Two new unit 1-forms R and S are introduced through the equations: 



R 



-1, A U K V = 0, (A*)^ = 0, ( U R V = e, S = C + eR- 



The only solution of the first 4 equations in the F-adapted coordinate system is R^ = (0,0,-1,0), 
so, (R) 2 = —1. Then for S we obtain S M = (0,0,0, 1), (S) 2 = 1. This pseudoorthonormal co-tangent 
basis (A, A*,R, S) is carried over to a tangent pseudoorthonormal basis by means of rj. 

In this way at every point, where the field is different from zero, we have three frames: the 
pseudoorthonormal (F-adapted) coordinate frame (dx,dy,dz,d^), the pseudoorthonormal frame x° = 
(A, eA*, R, S) and the pseudoorthogonal frame x = (A eA* , R, S). The matrix x^u of x with respect 
to the coordinate frame is 

u —p 
p u 
0-10 
1 



We define now the amplitude 4> > of the solution F by 



\det(x, 



Clearly, in an F-adapted coordinate system 4> = \f u 2 + p 2 
Each of these 3 frames defines its own volume form: 



(40) 



\OtF\ 



\A\. 



lo = dx A dy A dz A d^; uj x ° = A A eA A R A S = —uj; uj x = A A A* A R A S = -{u 2 + p 2 )uJ. 

We proceed further to define the phase of the nonlinear solution F. We shall need the matrix x of 
the frame x° with respect to the coordinate basis. We obtain 



sj v?+p 2 
V 



sj u 2 +p 2 






u 

yju 2 +p 2 











-1 
1 



The trace of this matrix is 



Hxl 



flUJ 



2u 
x/u 2 +p 2 



Obviously, the inequality |^ir(x° )| < 1 is fulfilled. Now, by definition, the quantities ip and ip defined 
by 



<P = \ tr (.X%), tP = arccos(^) = arccos Qtr 



(X 



(41) 



will be called the phase junction and the phase of the solution, respectively. Clearly, in the F-adapted 
coordinate system <p = cos[arg(ap)]- 

Making use of the amplitude function (ft, of the phase function (p and of the phase ip, we can write 



(ficostp, p = <p(±\fl 



>sin^. 



(42) 
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We consider now the pair of 1-forms (A, A*) on a region U C M, where A / 0, A* / 0, and since 
A. A* = 0, we have a 2-dimensional Pfaff system on U. 

Proposition 8. The above 2-dimensional Pfaff system (A, A*) is completely integrable, i.e. the 
following equations hold: 

dA A A A A* = 0, dA* A A A A* = 0. 

Proof. In fact, AAA* = (u 2 + p 2 )dx A dy, and in every term of dA and d^4* at least one of the 
basis vectors dx and dy will participate, so the above exterior products will vanish. The proposition 
is proved. 

Remark. As we shall see further, this property of Frobenius integrability is dually invariant for the 
nonlinear solutions. 

Remark. These considerations stay in force also for those time stable linear solutions, which have 
zero invariants I± = I2 = 0. But Maxwell's equations require u and p to be infinite running plane 
waves in this case, so the corresponding amplitudes will NOT depend on two of the spatial coordinates 
and the phase functions will be also running waves. As we'll see further, the phase functions for the 
nonlinear solutions are arbitrary bounded functions. 

We proceed further to define the new and important concept of scale factor C for a given nonlinear 
solution F with \SF\ / 0. First we recall a theorem from vector bundle theory, which establishes 
some important properties of those vector bundles which admit pseudoriemannian structure [30]. 
The theorem says that if a vector bundle £ with a base manifild B and standard fiber V admits 
pseudoriemannian structure g of signature (p,q),p + q = dimV, then it is always possible to introduce 
in this bundle a riemannian structure h and a linear automorphism ip of the bundle, such that two 
subbundles £ + and £~ may be defined with the following properties: 

1. </(£+, £+) = / l (£+,£+), 

2. «/(£-,£-) = -/»(£-£-), 

3. #(£+,£-) = /*(£+£") = 0. 

The automorphism 92 is defined by 

9x(u x , v x ) = h x (ip(u x ), v x ), u x , v x e T X B, x € B. 

In components we have 

9ij = Kw) <Pi = 9imh mk . 

In the tangent bundle case this theorem allows to separate a subbundle of the tangent bundle if the 
manifold admits pseudoriemannian metric. In the simple case of Minkowski space (M, 77), sign(rj) = 
(—,—,—,+), introducing the standard Euclidean metric h on M we may separate 1-dimensional sub- 
bundle, i.e. a couple of vector fields ±A D , being eigen vectors of 99, and to require rj(X , X Q ) = 1. In 
our canonical coordinates we obtain X Q = d/d£, so, the coordinate free relativistic definitions of the 
electric E and magnetic B 1-forms through F and *F are: 

E = -i(X )F, B = i(X )(*F). 

Now, the scale factor C for a nonlinear solution F with \SF\ ^ is defined by 

\i(X )F\ _ \A\ _ \i(X )(*F)\ _ \A*\ 
\SF\ \SF\ \S*F\ \S*F\ 

24 



(43) 



It is interesting to see how C is denned in non-relativistic terms. Consider the vector fields 



<9B E x B 

^ = rotE + - + ^-^divB, (44) 

-» <9E E x B 

M = rot B - — - — — — div E. (45) 
<9£ |E x B| 

It is obvious that on the solutions of Maxwell's vacuum equations T and M. are equal to zero. 
Note also that under the transformation (6) we get T — > — M and M — ► T. 

We shall consider now the relation between T and E, and between M and B on the nonlinear 
solutions of our equations when J- ^ and M ^ 0. 

Recalling eqns.(ll)-(16) we obtain 

(E x B) x E = -E x (E x B) = — [E(E.B) - B(E.E)] = B(E 2 ), 
and since |E x B| = E 2 = B 2 , we get 

rotE + — J x E + BdivB = 0. 

In the same way we get M x B = 0. In other words, on the nonlinear solutions we obtain that T is 
co-linear to E and M is co- linear to B. Hence, we can write the relations 

T = /i.E, M = f 2 .B, (46) 

where f\ and f 2 are two functions, and of course, the interesting cases are /i / 0,oo; f 2 / 0, oo. 
Note that the physical dimension of fi and f 2 is the reciprocal to the dimension of coordinates, i.e. 
[fi] = [f2] = [length]- 1 . 

We shall prove now that fi = f 2 - In fact, making use of the same formula for the double vector 
product, used above, we easily obtain 

f xB + M xE = 

(<9B \ / <9E \ 

rotE + — J x B + I rotB - — J x E - EdivB - BdivE = 0. 

Therefore, 

f xB + jl4xE = 

= / 1 ExB + / 2 BxE = (h - / 2 )E x B = 0. 

The assertion follows. The relation \F\ = \M\ also holds. It is readily verified now that the scale 
factor C may be defined by 

£<E3) - - if, - § - J§L («) 
l/2| |^| \M\ 

We go back to the 4-dimensional consideration. Obviously, for all nonlinear solutions F with 
\5F\ / we have < C(F) < 00, and for all other cases we readily obtain C — > 00. From the 
definition (43) of C is seen that C depends on the point in general. Note that if the scale factor £, 
defined by the nonlinear solution F, is a finite and constant quantity, we can introduce a characteristic 
finite time-interval T(F) by the relation 

cT(F) = £(F), 
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and also, in some cases, the corresponding characteristic frequency by 

u{F) = 1/T(F). 

In these wave terms the modified scale factor 2irC acquires the meaning of wave length. 

When C = const / we can say that the inverse scale factor C~ 1 is a measure of "how far" a 
nonlinear solution F (with \SF\ / 0) is from the linear Maxwell solutions. 

We modify now the volume form u) x to (—-uj x ). The integral of (— -oj x ) over the 4-region 
M 3 x C, C = const, is well defined for all 3d-finite solutions, and its value is an integral invariant. 
Making use of an F-adapted coordinate system we readily obtain 

-- / u x = - [ (u 2 + p 2 )dx A dy A dz A d£, = -EC = ET, 

where E is the integral energy of the solution in this frame. Note that this integral invariant depends 
only on the 3d- finite nature of the solution and on the finite constant value of the scale factor C. 

5 EED: Explicit non-linear vacuum solutions 

As it was shown with every nonlinear solution F of our nonlinear equations a class of ^-adapted 
coordinate systems is associated, such that F and *F acquire the form (39): 

F = eudx Adz + udx A d£ + epdy Adz + pdy A d£ 
* F = —pdx Adz — epdx A d£ + udy Adz + eudy A d£. 

Since we look for non-linear solutions of (18)- (20), we substitute these F and *F and after some 
elementary calculations we obtain 

Proposition 9. Every F of the above kind satisfies the equation 5F A F — 5 * F A *F = 0. 
Proof. It is directly verified. 
Further we obtain 

5F = - eu z )dx + (p^ - ep z )dy + e(u x +p y )dz + (u x +p y )d^, 

5* F = -e(p 5 - ep z )dx + e{u^ - ep z )dy - (p x - u y )dz - e(p x - u y )d£, 
F^{8F) v dx v = (*F)^(5 * F) v dx v = 
= £ \p(Pt ~ tPz) + u(u£ - eu z )) dz + [p(p£ - ep z ) + u(u^ - eu z )} d£, 
(SF) 2 = (5 * F) 2 = -(u ( - eu z ) 2 - fa - ep z ) 2 = -<f> 2 (^ - e^ z ) 2 . 
We infer that our equations reduce to only one equation, namely 

P(P* " Wz) + <H - £u z) = 2 [(" 2 + P\ ~ £ ( u2 + V 2 )z\ = 0. (48) 

The obvious solution to this equation is 

u 2 +p 2 = 4> 2 {x,y^ + ez), (49) 

where (j) is an arbitrary differentiable function of its arguments. The solution obtained shows that 
the equations impose some limitations only on the amplitude function <p = \olf\ and that the phase 
function ip is arbitrary except that it is bounded: \(p\ < 1. The amplitude is a running wave along 
the specially chosen coordinate z, which is common for all F-adapted coordinate systems. Considered 
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as a function of the spatial coordinates, the amplitude <f> is arbitrary, so it can be chosen spatially 
finite. The time-evolution does not affect the initial form of <p, so it will stay the same in time, but the 
whole solution may change its form due to 92 . Since \ip\ < 1 and the two independent field components 
are given by F14 = u = <p(p, F24 = p = ±<fi \J\ — ip 2 this shows, that among the nonlinear solutions of 
our equations there are (3+1) soliton-like solutions. The spatial structure of <p can be determined by 
initial condition, and the phase function <p can be used to describe additional structure features and 
internal dynamics of the solution. 

We compute SFAF = S*FA*F and obtain 

SF A F = ep(u£ — eu z ) — eu(p^ — ep z ) dx A dy A dz + p(u^ — eu z ) — u(p^ — ep z ) dx Ady A d£. 
In terms of <p and ip = arccos((/?) we obtain 

p(v,£ - eu z ) - u(p^ - ep z ) = -0 2 (V>£ - eipz), 

and so 

SF A F = —e<p 2 {ip£ - eip z )dx Ady Adz - (p 2 (ip£ - eip z )dx Ady A d£. 
Applying * from the left we get 

*(SF A F) = -<P 2 (^ - eip z )dz - e<p 2 {^ - eip z )dt = -e<P 2 {^ - eip z )(. 



Corollary. 5F A F is isotropic, and it is equal to zero iff ip is a running wave along z. 

Corollary. F is a running wave along z iff ip is a running wave along z, i.e. iff SF A F = 0. 

Computing the 4-forms d A A A A ( and d^L* A A* A £ we obtain 

dA A A A C = d^L* A A* A ( = e [u(p£ - ep z ) - p(u^ - eu z )] dx A dy A dz A d£ 

= e<j) 2 (ip^ — eip z )dx Ady Adz A d£. 

Hence, since (SF) 2 = —(p 2 (ip^ — eip z ) 2 , and <p 7^ is a running wave, the nonlinear solutions F satisfying 
\SF\ = \S * F\ = 0, imply ip£ — eip z = 0, i.e. absence of rotational component of propagation. The 
following relations are equivalent: 

1. SF AF = 0. 

2. \SF\ = \S*F\ = 0. 

3. ip is a running wave along z: L^ip = 0. 

4. dA A A A C = dA* A A* A ( = 0. 

We give now two other relations that are equivalent to the above four. First, consider a nonlinear 
solution F, and the corresponding (1, 1) tensor F^ = n UCT F CT ^. We want to compute the corresponding 
Frblicher-Nijenhuis tensor Sf = [F, F] , which is a 2-form on M with values in the vector fields on M. 
The components of Sf in a coordinate frame are given by 

r)F a dF a r)F a dF a ~ 

pa ur v _ jpa /f _ pa Ur v 1 pa l» 

» dx a v dx a a dx» a dx" ' 

We recall the two unit vector fields A and eA*, given by (in a F-adapted coordinate system) 

- d 1 77 d ~ 1 ~ d d 

A = -ip— - Vl-ip 2 — , sA* = ^l-<p 2 —-ip—, 
ox ay ox oy 

and we compute Sf{A, eA*). 

(S F )°A^eA* u = (Sp^i^eA* 2 - A 2 eA* 1 ). 



{ s f)% = 2 
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For (Sp)i2 we get 

(S F )\ 2 = (Sf)1 2 = 0, (S F f 12 = -e{S F )\ 2 = 2e{ P { H - eu z ) - ufa - e Pz )}. 

It is easily seen that A 1 eA* 2 — A 2 sA* 1 = 1, so, the relation Sf(A,eA*) = is equivalent to the 
above four. 

Second, recall that if [A, (+, •)] is an algebra (may graded), then the (anti) derivations V : A — > A 
satisfy: V(a.b) = Va.b + e a a.Vb, where a, b £ A and e a is the parity of a € A- So, the derivations are 
not morphisms of A, and satisfy the generalized Leibniz rule. The difference 

A D (a,b) = Da.b + e a a.Db - D(a.b) ,a,b£A 

is called the Leibniz bracket of the operator D, and D is a (anti) derivation if its Leibniz bracket 
vanishes. If A is the exterior algebra of differential forms on a manifold M and D is the coderivative S 
with respect to a given metric, then the corresponding Leibniz bracket (sometimes called the Schouten 
bracket) is denoted by {, }. So, if F, G are two forms on M then 

{F, G} = 5FAG + (-1) P F ASG — 5(F A G). 

Note that the brackets {F, F} do not vanish in general. 

Now, if the 2-form F on the Minkowski space is a nonlinear solution, then FAF = and 

{F, F} =5F AF+ {-l) 2 F A5F- S(F A F) = 25F A F. 

So, the above relations are equivalent to the requirement that the Leibniz/Schouten bracket {F, F} 
vanishes. 

Finally, all these conditions are equivalent to C = oo. 

We note the very different nature of these seven conditions. The complete integrability of any of 
the two Pfaff 2-dimensional systems (A, Q and (A*, () is equivalent to zero value of \SF\ on the one 
hand, and to the zero value of the quantity Sf(A,eA*) on the other hand, and both are equivalent 
to the vanishing of the Leibniz/Schouten bracket and to the infinite value of C This could hardly 
be occasional, so, a physical interpretation of these quantities in the nonzero case, i.e. when -ip is not 
a running wave, is strongly suggested. In view of the above conclusion that the condition \SF\ = 
implies absence of rotational component of propagation, our interpretation is the following: 

A nonlinear solution will carry rotational component of propagation, i.e. intrinsic an- 
gular (spin) momentum, only if tp is NOT a running wave along the direction of translational 
propagation. 

Natural measures of this spin momentum appear to be 5FAF / 0, or \SF\ / 0, or Sf(A, sA*). The 
most attractive seems to be 5FAF, because it is a 3-form, and imposing the requirement d(5F AF) = 
we obtain both: the equation for tp and the corresponding conserved (through the Stokes' theorem) 
quantity Ti = J i*(5F A F)dx A aly A dz, where i*(5F A F) is the restriction of 5F A F to IR 3 . 

Making use of the relations u = <fi cosip and p = (ftsimp, we get 

A = <p cos ipdx + 4> sin ipdy + f(, A* = —e<j) sin ip dx + ecf) cos ip dy + /*C, 
SF = -(f> sin ip (?/>§ - eip z ) dx + <f> cos ip (V>£ - £ip z ) dy + (u x + p y )( 

= -e(L^)A* + (epL^ + u x +p y )t (50) 
5 * F = —ecp cos tp (ip£ — eip z ) dx — e<p> sin tp (ip£ — eip z ) dy — e{p x — u y )( 

= e(L c i/>)A - e(fL c tp + Px - u y )(, 
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where is the Lie derivative with respect to £. We obtain also 



\SF\ = \S*F\= zn^J^l = _ ei, z \ = <j,\L^\, 
V 1 - r 



(51) 



\bF\ \<pz~eip z \ \^-eif>. 
Now we can write 



SF = ±e^ + {^e f — + (u x + p y )^j C, 5 * F = + (±e^ - e( Px - u y fj (. (52) 



Corollary. Obviously, the following relations hold: 

A^5Fr = 0, (A*) li (5*Fy = 0, A li (6*Fy = ±e^ = -(A*) li (6Fy. 

Finally we note that since the propagating along the given £ nonlinear solutions in canonical 
coordinates are parametrized by one function <p of 3 independent variables and one bounded function 
(p of 4 independent variables, the separation of various subclasses of nonlinear solutions is made by 
imposing additional conditions on these two functions. 



6 Homological Properties of the Nonlinear Solutions 

From pure algebraic point of view we speak about homology (or, cohomology) every time when we 
meet a linear map D in a vector space V over a field (e.g. R, or C), or in a module W over some ring, 
having the property D o D = 0. Then we have two related subspaces, Ker(D) = {x € V : D(x) = 0} 
and Im(D) = D(Y). Since Im(D) is a subspace of Ker(D), we can factorize, and the corresponding 
factor space H(D,Y) = Ker(D) / Im(D) is called the homology space for D. The dual linear map D* 
in the dual space V* has also the property D* o D* = 0, so we obtain the corresponding cohomology 
space H*(D*,Y*). In such a situation the map D (resp. D*) is called boundary operator (resp 
coboundary operator). The elements of Ker(D) (resp. Ker(D*)) are called cycles (resp. cocycles), 
and the elements of Im(D) (resp. Im(D*)) are called boundaries (resp. coboundaries) . 

The basic property of a boundary operator D is that every linear map B : V — > V which commutes 
with D: D o B = B o D, induces a linear map B* : H(D) — > H(D). So, a boundary operator realizes 
the general idea of distinguishing some properties of a class of objects which properties are important 
from a definite point of view, and to find those transformations which keep invariant these properties. 

The basic example for boundary operator used in theoretical physics is the exterior derivative d 
(the so called de Rham cohomology). This operator acts in the space of differential forms over a 
manifold, e.g. the Euclidean space (R 3 ,g), or the Minkowski space-time M = (IR 4 , 77) . The above 
mentioned basic property of every boundary operator D appears here as a commutation of d with the 
smooth maps f of the manifold considered: d o f* = f* o d. A well known physical example for cocycles 
of d (which are called here closed differential forms) comes when we consider a spherically symmetric 
gravitational or electrostatic field generated by a point source. Since the field is defined only outside 
of the point-source, i.e. on the space iV = IR 3 — {0}, then a natural object representing the field is a 
closed differential 2- form u, du; = 0. In view of the spherical symmetry such a (spherically symmetric 
or S'0(3)-invariant) closed 2- form is defined up to a constant coefficient q, and in standard spherical 
coordinates originating at the point-source we obtain uj = qsvnQdO A dip. The Hodge star on iV gives 
*uj = (q/r 2 )dr, which is usually called electric field generated by the point source q. Now, the Stokes' 
theorem establishes the charge q as a topological invariant, characterizing the nontrivial topology of 
the space N. 
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The above mentioned boundary operator d is a differential operator. Pure algebraic boundary 
operators in a linear space V also exist and one way to introduce such operators is as follows. Let V* 
be the dual to V space and (,) denotes the canonical conjugation: (x*,x) — > (x*,x), where x* G V* 
and x G V. Now fix i £ F and let x* be such that (x*,x) = 0. Consider now the decomposable 
element x* ® x £ V* ® V . The isomorphism between V* <8> V and L(V, V) allows to consider x* (g> x 
as a linear map </?( x * )X ) : — ► F as follows: 

¥>(s*,x)(y) = (x*®x)(y) = (x*,y)x, y € V. 

Hence, the linear map v?(x*,x) sends all elements of V to the 1-dimensional space determined by the 
non-zero element x 6 V. Clearly, <£>( X * ;X ) is a boundary operator since 

¥>(x*,x) °¥>(**,x)(y) = VV,*)^*'^) = (»*,J/)y(x',x)(^) = (x*,y){x*,x)x = 0, i/eK 

So, if {ei} and {e- 7 } are two dual bases in the n-dimensional linear spaces V and V* respectively, then 
every couple {e 3 , ej), i 7^ j, determines the linear boundary operators <P( E j >e A = £ J ® e%,i ^= j- 

If g is an inner product in V and the two elements (x, y) are (/-orthogonal: g(x, y) = 0, then 
denoting by g the linear isomorphism g : V — > V* (lowering indices) we can define the linear map 
g(x) ® y, which is obviously a boundary operator in V: 

(g(x) «t/)o (g(x) ®y)(z) = (g(x),z)(g(x) ®y){y) = (g{x), z)g(x,y)(y) =0, z G V. 

In particular, if (V, r/) is the Minkowski space, then every isotropic vector £ : Q = Q 2 = defines 
a boundary operator 77(C) <8> C = C ® C- I n f ac t> 

(r](C) C) o (r](C) C) (x) = (C ® C) o (C ® C)(x) = (C, x)(( C)(C) = <C, x)((, 0(0 = 0, x £ V. 

We note that the duality between V and V* allows to consider the element x* ® x as a linear map 
in as follows (x* <8> x)(y*) = (y*,x)x*, and to build the corresponding boundary operators in V*. 
Also, in the finite dimensional case we always get dimiV) — 2 = dim[Ker(x* <S> x) / Im(x* <8>x)], where 
< x*,x >= 0. 

The above mentioned property that the image space of any such boundary operator x* <g) x, where 
(x*,x) = 0, is 1-dimensional, implies that the natural extensions of these boundary operators to 
derivations in the graded exterior algebras A(V) and A(V*) [28] define boundary operators of degree 
zero in these graded algebras (see further). 

Further in this section we show that the energy-momentum tensors of the nonlinear vacuum so- 
lutions in EED give examples of such boundary operators, and we also present some initial study of 
the corresponding homology. We have the following property of the electromagnetic stress-energy- 
momentum tensor Q(F) 1 ^ established by Rainich [29] 

Q{F)lQ{FZ = \(ll + ll)5^ (53) 

where I\ = \F [IV F IIV and I2 = \F^ U (*FY U are the two invariants. Considering Q(F) 1 ^ as a linear 
map in the module of vector fields, or 1-forms, over the Minkowski space-time, we see that if the two 
invariants are equal to zero, as it is in the case of nonlinear solutions of our equations, we obtain 
a boundary operator at all points of M where F / 0. Our purpose now is to consider how the 
corresponding homology is connected with the structure of the nonlinear solutions F of the vacuum 
EED equations. 

Proposition 10. The image space Im^Qp) coincides with the only isotropic eigen direction of 
Qf- 
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Proof. The linear map Qp is given in an F-adapted coordinate system by 



Qf = ~4> 2 dz ®^-+ s4> 2 dz e<p 2 d£ <g> ^- + 4> 2 d£ <g> (54) 

a£ a£ 

Clearly, Qf = <fi 2 ( <8> C> where £ = s^-z + ^ in the F-adapted coordinate system and £ = — ed z + 
defines the only isotropic eigen direction of Qf- Let in this coordinate system the arbitrary vector 
field X be presented by its components {X^), fi = 1, ... ,4. We obtain 



O r\ Pi Pi 

Q F (X) = -«/, 2 ^(X) — + £0 2 ^(A) — - e</> 2 dt(X) — + <i> 2 dti{X) — 

= ^ 

= e0 2 X 3 C + </> 2 X 4 C = </> 2 (eA 3 + X 4 )(. 
If a = a^dx^ is a 1-form then in the same way we obtain 

Q F (a) = (QF)^a u dx^ = <fi 2 (-ea 3 + a 4 )C- 

The proposition is proved. 

Proposition 11. The kernel space Ker(Qp) coincides with the 3-space spanned by the vectors 
A, A* and £. 

Proof. We know that Qf{A) = Qf{A*) = Qf(C) = 0. Now, if X is an arbitrary vector, then 
from the above Prop.10 follows that Q F {X) = <j) 2 (eX 3 + X 4 )(, so we conclude that Qp(X) will be 
equal to zero only if X is a linear combination of A, A* and The proposition is proved, 

Hence, we may write Ker(Qp) = {A} © {^4*} © {C}- The corresponding factor space 

H(Q F ) = K er(Q f) / Im(Q f) 

is isomorphic to {A} © {A*}. The classes defined by A and A* are given by [A] = A + f( and 
[A*] = A* + /*£, where / and /* are functions. Recall now that Qf(SF) = Qp(S * F) = 0, so 
[SF] = 5F + h( and [S * F] = S * F + h*(, where h and h* are functions, give the corresponding 
homology classes. 

According to the above mentioned property, every symmetry of a boundary operator induces a 
linear map inside the homology space. Therefore, the homology spaces are invariant with respect to the 
linear isomorphisms which commute with the boundary operators. In our case we have to find those 
linear maps $ in the module of vector fields over M, which commute with Qf, i-e. Qp = QpoQ. 
It is readily obtained that in the F-adapted coordinate system every such <3? is given by a matrix of 
the following kind: 



$ = 



a 


b 


c 


—EC 


m 


n 


q 


-eq 


r 


s 


w 





er 


es 





w 



where all nine independent entries of this matrix are functions of the coordinates. It follows that the 
Qp- homology spaces are invariant with respect to all diffeomorphisms (p : M — > M which generate 
isomorphisms dip : TM — > TM of the tangent bundle of M given in the F-adapted coordinate system 
by a nondegenerate matrix of the above kind. 

An important property of the boundary operator Qp is that its image space Im{Qp) is 1- 
dimensional. As it was mentioned earlier, this allows to extend Qf as a boundary operator in the 
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graded exterior algebra of differential forms over M. In fact, recall that a linear map <p in a linear 
space V induces derivation ip A in the exterior algebra A(V) according to the rule 

(p A (x\ A X2 A • • • A x p ) = <f(xi) A a?2 A • • • A x v + x\ A <p(x2) A • • • A x p + • • • + x\ A a?2 A • • • A <p(x p ). 



Remark. If we try to extend Qf to antiderivation with respect to the usual involution 
ui{a) = (— l) p a, a 6 A P (M) we'll find that this is not possible since the necessary condition for 
this, given by Qf{ol) A a + ui{a) A Qf{q) = 0, does not hold for every a G A 1 (M). 

Hence, if Im(<p) = ip(V) is 1-dimensional, then every summond of (p A o (^ A (xi A X2 A • • • A x p ) will 
contain two elements of the kind (p(xi) and <p(xj), and if these two elements are collinear, their exterior 
product is zero and the corresponding summond is zero. In our case ip = Qf and Im(QF) = {(} is 
1-dimensional, so we shall have Q F o Q F (a) = 0, a € A(M). 

Corollary. The extension defines a boundary operator of degree zero in A(M). 

Remark. Further the extension Q F will be denoted just by Qf- 

Corollary. The extension of Qf to derivation in A(M) introduces in A(M) some structure of 
graded differential algebra with corresponding graded homology algebra H(Qf)(A(M)). 

The following relations are readily verified: 

Q F (F) = 0, Q f (*F) = 0, Q F (5FAF) = 0. 

For example, Q F {F) = Q F {Af\Q = Qf(A) A( + AaQ f (() = 0. So, F,*F and 5FAF are Q F -cycles. 
From relations (50) we obtain 

[6F] = -e(L^)[A*], [5 * F] = e(L^)[A], 

i.e. SF and A* define the same QF-homology classes, and 6 * F and A define the same QF-homology 
classes in A 1 (M). 

Proposition 12. The scale factor C = |^4|/|<5.F| depends only on the classes of A and 5F. 
Proof. Since p]| = \A + f(\ = \A\ and \[5F}\ = \5F + h(\ we obtain C = \A\/\8F\ = \[A]\/\[5F}\. 

Clearly, ([A], [A*]) and ([8F], [6 * F]) represent two bases of H(Q F ) in A 1 (M). 

Corollary. The transformation ([A], [A*]) — > ([<5F], [5 * F]) is given by 



eL^ 
-eL c ^ 



(-eL^[A*],eL^[A}) = ([6F],[6*F}). 



(55) 



The above formula (55) shows that the transformation matrix, further denoted by Ai, between these 
two bases is (±eC)~ 1 J, where J is the canonical complex structure in a real 2-dimensional space. This 
fact may give another look on the duality symmetry (Recall Prop. 3), because of the invariance of J 
with respect to the transformation J — > S.J.S^ 1 , where S is given on p. 7 : 



a b 




1 




a 


-b 


1 


1 


—b a 




-1 




b 


a 


a 2 + b 2 


-1 



One could say that the duality symmetry of the nonlinear solutions is a consequence of the null-field 
homology presented. In other words, every initial null-field configuration given by (A, A*, Q, with 
\5F\ / 0, \5 * F\ ^ 0, £ 2 = 0, compulsory has rotational-translational dynamical nature, so, it is 
intrinsically forced to propagate with rotational component of propagation in space-time, because, 
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according to relation (55), the nonzero F^, i.e. the nonzero (A, A*), imply nonzero values of the 
derivatives of F^ v including the nonzero value of L^ip even if ip is time-independent, and the basis 
([A], [A*)] is continuously forced to rotate. In fact, the running-wave character of <fr = \A\ drags the 
solution along the coordinate z and the nonzero L^tp implies cos?/' / 0, sin?/' / 0. The evolution 
obtained is strongly connected with the nonzero finite value of the scale factor C = |L^V| _1 , which, 
in turn, determines rotation in the homology space Hp{Q). This rotation is determined entirely by 
the Lie derivative of the phase i/j with respect to £, and it is intrinsically consistent with the running 
wave translational propagation of the energy-density 4> 2 . It is seen that the field configuration has 
a rotational component of propagation, while the energy density has just translational component of 
propagation. 

It is interesting to see the action of Qf as derivation in A(M). We shall do this in a F-adapted 
coordinate system. Let F define a nonlinear solution and Qf '■ A(M) — > A(M) be the corresponding 
derivation with 4> 2 = u 2 + p 2 the corresponding energy-density. We give first the action of Qf as 
derivation on the bases elements. 

Q F {dx) = Q F {dy) = 0, Q F {dz) = -e<p 2 (, Q F (dO = 4>\ 

QF(dxAdy) = 0, Qp{dx A dz) = -e<j) 2 dx A £, Qp(dy A dz) = -ecj) 2 dy A (, 

QF{dx A d£) = 4> 2 dx A C, Q F {dy Ad£) = 4> 2 dy A(, Q F {dz A d£) = 0. 

QF(dx A dy A dz) = —e<j) 2 dx A dy A (, Qp(dx A dy A d£) = <j) 2 dx A dy A (, 
Q F (dx A dz A d£) = 0, Q F (dy A dz A d£) = 0, 

QF{dx A dy A dz A d£) = 0. 

Let now the arbitrary 2-form G be represented in this coordinate system by G = G^ u dx^Adx u , fi < 
v. Making use of the above given explicit form for the action of Qf on the basis elements as derivation 
we obtain 

Q F (G) = e<t> 2 (-eG 13 + G 14 )dx A + e0 2 (-eG 23 + G 2 i)dy Adz+ 

o o (56) 

<p 2 (-eG 13 + G u )dx A d£ + 2 (-eG 23 + G 2 ±)dy A d^. 

This result makes possible the following conclusions concerning 2-forms: 

1. The space Im(QF) consists of null fields, i.e. every nonlinear solution F determines a 
subspace Im(QF) C A 2 (M) of null-fields . 

2. The space Ker{Qp) consists of 2-forms, which in this coordinate system satisfy: eG\ 3 = 
G14, eG23 = G24, and (G12, G34)-arbitrary. 

3. The eigen spaces of Qf(G) coincide with the eigen spaces of F for every (nonzero) 

G e A 2 (M). 

If G is a 3- form with components G123, G124, G134, G234 in the same F-adapted coordinate system, 
we obtain 

Qf(G) = £0 2 (-eGi23 + G 12 4)dx A dy A dz + 4> 2 {-eG 123 + G l2 A)dx A dy A d£. 

So, Qf(G) is isotropic, and a 3-form G is in Ker(Qp) only if eG\ 23 = G124 in this coordinate system. 
Moreover, since Qp{G) does not depend on Gi^dx Adz Ad^ and G 2 udy Adz Ad^ we conclude that 
the kernel of Qf in this case consists of time-like 3-forms. 
Finally, if G is a 4- form, then Qf{G) = 0. 

Corollary. If G G A(M) lives in Im{QF) where F is a nonlinear solution, then Qf{G) is isotropic. 



33 



This may be extended to the smooth functions / G C°°(M) if we assume QfU) = 0- 
Let's resume. Every space-like (straight-line) direction may be chosen for z-coordinate on M, and 
the 1-form £ = edz + d£ determines an isotropic direction along which a class of null-fields F are 
defined by the formula (39). The corresponding linear map Qf satisfies Qf ° Qf = and defines 
homology in the spaces of 1-forms and of vector fields. 

Since the image space Qf(^M) is 1-dimensional, Qf extends to a boundary operator in the whole 
exterior algebras over the 1-forms and vector fields. The image space of the extended Qf consists of 
isotropic (null) objects. If G is a 2-form then Qp(G) has, in general, the same eigen properties as F. 
Hence, every 2-form F with zero invariants lives in just one such subclass and the whole set of these 
2-forms divides to such nonintersecting subclasses. Moreover, every 2-form G has its (null-field) image 
in every such subclass. 

For every nonlinear solution F(u,p), (5F / 0) the corresponding (f> 2 = u 2 + p 2 propagates trans- 
lationally, i.e. is a running wave, along the space-like direction chosen (considered as the coordinate 
z). The 4-dimensional versions of the corresponding electric and magnetic fields are presented by 
the nonisotropic parts of the homology classes defined by the mutually orthogonal space-like 1-forms 
A and (—A*). The nonlinear solutions with \SF\ = \S * F\ = propagate only translationally, i.e. 
without rotation. Rotational components of propagation, or spin-momentum, may have just those 
nonlinear solutions having nonzero finite scale factor C = \A\/\5F\, or equivalently, satisfying one 
of the conditions given above (pp. 27-28). The isotropic 3-form 5F A F defines a QF-homology class 
since Qf(SF A F) = 0, and it appears as a natural candidate representing locally the spin-momentum 
if we assume the additional equation d(5F A F) = 0, which should reduce to an equation for the 
phase ip. The two mutually orthogonal space- like 1-forms 5F and 5 * F define the same homology 
classes as A* and A respectively. The transformation matrix M. between the two bases (A, A*) and 
(5 * F, SF) defines a complex structure in the 2-dimensional homology space through the scale factor: 
M[A] = M(A + f() =M(A) + f( = 5*F + f(= [5*F] and M = ±e£~ 1 J, where J is the canonical 
complex structure in a 2-dimensional space. The 2-parameter duality symmetry coincides with the 
symmetries of M. and transforms solutions to solutions inside the subclass of solutions propagating 
along the spatial direction chosen. 

7 Structure of the Nonlinear Solutions 

We return now to the correspondence F — > op. 

7.1 Properties of the duality matrices 

We consider the set G of matrices a of the kind 



a = 



a b 
—b a 



where a, b £ 



(57) 



As we mentioned earlier, the nonzero matrices of this kind form a group G with respect to the usual 
matrix multiplication. Together with the zero 2x2 matrix I a they also form a 2-dimensional linear 
space Q = G {J{I } over R with respect to the usual addition of matrices. As is well known Q gives 
the real representation of the field of complex numbers. 

A natural basis of the linear space Q is given by the two matrices 



I 



The group G is commutative, in fact 



1 
1 



J 



1 
-1 



a b 




m n 




m n 




a b 




am — bn (an + bm) 


—b a 




—n m 




—n m 




—b a 




— (an + bm) am — bn 
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Every element a £ Q can be represented as a = al + bJ, a, b G R. Recall the natural representation 
of G in Q given by 

o«(J) + ba(J) . 

From now on we shall consider a(I) and a(J) just as matrix product, so we have 

a(I) = a. I = a = al + bJ, a(J) = a. J = —bl + aJ. 

Since (J" 1 )* = J, J generates a complex structure in Q: Jo J(x) = — x, x E Q. 

The product of two matrices a = al+bJ and (3 = ml+nJ looks like a. (3 = (am—bn)I+(an+bm)J. 
The commutativity of G means symmetry, in particular, every a £ G is a symmetry of J: a.J = J. a. 

Finally we note, that the inner product g e in Q = T e (G), where e is the identity of G, given by 

g e {a,(3) = hr{ao(3*), a,(3eG, 
generates a (left invariant) riemannian metric on G by means of the (left) group multiplication: 

9a = {L a -i)*g e , a € G. 

7.2 The action of G in the space of 2-forms on M 

We consider now the space A 2 (M) - the space of 2-forms on M with its natural basis: 

dx A dy, dx A dz, dy A dz, dx A d£, dy A d£, dz A d£. 

We recall that the Hodge * acts in A 2 (M) as a complex structure * = J and on the above basis its 
action is given by: 

J(dx A dy) = -dz A d£ J(dx A dz) = dy A d£ J(dy A dz) = -dx A d£ 

J(dx A d£) = dy Adz J{dy A d£) = -dx A J^cfe A d^) = dx A dy. 

Hence, in this basis the matrix of J is off-diagonal with entries (1, —1, 1, —1, 1, —1). 

Let now J be the identity map in A 2 (M). We define a representation p of G in A 2 (M) as follows: 

p(a) = p(al + bJ) = al + bj, a G G. (58) 

The map p may be considered as a restriction of a linear map Q — ► L^i M \ to the nonzero elements of 
£7, i.e. to the elements of G. Every p(a) is a linear isomorphism, in fact, its determinant det|| / o(a)|| 
is equal to (a 2 + 6 2 ) 3 . The identity / of G is sent to the identity transformation X of A 2 (M), and the 
complex structure J of the vector space Q is sent to the complex structure J of A 2 (M) . This map is 
surely a representation, because p(a./3) = p(a).p((3), a, /3 G G. In fact, 

p(a.p) = p[(al + bj).(ml + nJ)] = 

p[(om — bn)I + (a£2n + 6m) J] = (am — bn)Z + (an + bm)J. 

On the other hand 

p(a).p((3) = p(al + bJ).p(mI + nJ) 

= (aX + bJ).{mX + n,7) = (am — 6n)X + (an + bm)J. 

We consider now the space A 2 (M, Q) of (/-valued 2-forms on M. Every such 2-form $7 can be 
represented as fi = F\ (g) / + X2 ® J, where F\ and X2 are 2-forms. We have the joint action of G in 
A 2 (M,G) as follows: 

[p(a) x a](U) = p(a).Xi (a" 1 )*^) + p(a).F 2 ® (a" 1 )*^. 



p(a)(a/ + 6J) = a(a^ 1 )*(d) + 6(a- i )*(J) 



1 

a 2 -I- 6 2 
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We obtain 

det(a)\p(a) x a](Q) = 
[(a 2 l + abJ)F x - {b 2 J + abl)F 2 ] ® I + [(6 2 J + a6X)iq + (a 2 J + a&JO-Ffe] J. 
In the special case Q = F <g) I + i/.-F 1 <8> J it readily follows that 

[/?(«) x a] (0) = a (59) 

In this sense the forms fl = F®/ + (8) J are equivariant with respect to the joint action of G. 
Explicitly for a general 2-form F we have 

p(a).F = al.F + bJ.F 

= (aF 12 + bF u )dx A dy + (aiq 3 - bF 2A )dx A dz + (aF 23 + bF u )dy A dz 
+ (aF u - oF 23 )d£ A d£ + {aF 2A + oFi 3 )dy Ad(+ (aF 34 - bF 12 )dz A d£. 

If F £ is of the kind (39) we modify correspondingly the representation as follows: p(a) = dl + ebj, 
and obtain 

p(a).F £ = e(au — bp)dx A dz + e(ap + bu)dy A dz 
+ (era — bp)dx A d£ + (era + bu)dy A d£. 

Recalling now equation (48) the above formula (60) shows that if F £ is a nonlinear solution then 
p{a).F £ will be a nonlinear solution if the quantity 

[(p(a).F £ ) 14 ] 2 + [{ P {a).F £ ) 24 \ 2 = (au - bp) 2 + (ap + bu) 2 

is a running wave along z. But this quantity is equal to (a 2 + b 2 ){u 2 +p 2 ) and since (a 2 + b 2 ) = const, 
we see that p{a).F £ is again a nonlinear solution for any a£G. In other words, the group G acts as 
group of symmetries of our nonlinear equations. Moreover, in view of the conclusions at the end of the 
preceding section, G acts inside every subclass of solutions defined by the chosen space-like direction 
(the coordinate z). Hence, if F is a nonlinear solution, we may write G.F C Qf(A 2 (M)), i.e. any 
orbit G.F lives entirely and always inside the subclass Qf(A 2 (M)). 

As for the phase tp of the solution (60) we obtain (in this coordinate system) 



ip = arccos 



e(au — bp) 



y/(a 2 + b 2 )(u 2 +p 2 ) 



TP(F £ (u,p))+iP(a(a,b)), 



where i[)(F £ (u,p)) is the phase of the solution F £ (u,p) and ip(a(a,b)) = arccos [a/Va 2 + b 2 \ is the 
phase of the complex number a = al + bJ. Now, since ip(a(a, b)) = const we obtain the 

Corollary. The 1-form dtp and the scale factor C = {L^tpl^ 1 are G-invariants: d(tp(F)) = 
d(V(p(a).F)), £{F) = £{p{a).F). 

We are going to show now that the action p(a) : F £ — > p(a).F £ can be generated by an appropriate 
action T:GxM->MofGon the space-time M. The term appropriate here means that the action 
T takes in view some isotropic straight-line direction, which will be considered to lie inside the planes 
(z,£). The action T is defined by 

(a(a, b); x, y, z, £) (a/, y', z, £) = ^ + , ^ ^ ^2 ' z ' ' a ( a > & ) e G - 

Since x = ax' + by', y = —bx' + ay' and writing again (x, y) instead of (x', y')) we obtain: 

T*F £ = e(au — bp)dx Adz + e(ap + bu)dy Adz + (au — op)dx Ac!(+ (era + bu)dy A d£ 
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We obtain also 

T*( = C, T*A = {au-bp)dx + (ap + bu)dy, T*(A*) = -e(ap + bu)dx + e(au - bp)dy, 

5(T*F) = [a(v,£ - eu z ) - b(p^ - ep z )]dx + [b(u^ - eu z ) + a(p 5 - ep z )]dy+ 

+e[a(u x +py) - b(p x - u y )]dz + [a(u x +p y ) - b(p x - u y )]d£, 

\T*A\ 2 = \T*(A*)\ 2 = (a 2 + b 2 )(u 2 + p 2 ), \5(T*F)\ 2 = (a 2 + b 2 )\5F\ 2 . 

The invariance of the scale factor follows directly now from the above relations. 

The corresponding fundamental vector fields Zj and Zj, generated by the basis (I, J) of Q are 

7 ( 9 m 9 \ 7 9 J- 9 

Zi = - [x— + y— , Zj = -y— + x- 



ox Oy J ox oy 

From the above corollary and from the commutation relation d o Lx = Lx o d it follows now that the 
phase ip satisfies the equations 

L Z/ (d^) =L Zj (d4>) = 0. 

The G-invariance of the scale factor C, considered in this coordinate system, implies the relations 
Z[(C) = Zj(C) = 0, which lead to the conclusion that C may depend only on (z,£) in this coordinate 
system. In fact, we obtain the linear system 

xC x + yC y = 0, -yC x + xC y = 0, 

and since x 2 + y 2 / 0, the only solution is C x = C y = 0, i.e. C = C(z, £). 
Now we have the following result: 

Corollary: The scale factor £ is a constant quantity: C = const, only if it is a Lorentz-invariant 
quantity. 

In fact, every constant is Lorentz invariant. On the other hand, the invariance of C with respect 
to rotations in the plane (x, y) is obvious in this coordinate system: 

d£(z,t) x 9C(z,0 _ 
dx dy 

Further, the generators of the proper Lorentz rotations are 

d d d d d d 

So, any of the relations L\(C) = 0, L,2(C) = implies that C does not depend on £, and, finally, the 
relation L^{C) = requires not-dependence on z. 

Corollary. If the scale factor is Lorentz invariant: C = const, then the 3-form F A 5F is closed: 
d(F A SF) = 0. 

In fact, from relation (51) it follows that if C = const then ip£ — eip z = const too. We recall that 

SF AF = -e<p 2 (ip£ - etp z )dx A dy A dz - (j) 2 {ip£ - etp z )dx A dy A d£. 

The assertion now follows from the notice that $ 2 is a running wave along z. Hence, when the scale 
factor C is Lorentz-invariant we obtain another conservative quantity, namely, the integral of the 
restriction of 5F A F on B 3 over the whole 3-space will not depend on time. 
We also note that under the action of a(a, b) we have 

A—>A' = [au- bp, ap + bu, 0,0], A* -> (A*)' = [ap + bu,-(au- bp), 0,0], 
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and this is equivalent to 

A' = aA-bA*, (A*)' = bA + aA*. 

Hence, in the F-adapted coordinate systems the dual transformation, as given above, restricts to 
transformations in the (x, y)-plane, so we have the 

Corollary. The Frobenius integrability of the 2-dimensional Pfaff system (A, A*) is a G-invariant 
property. 

Corollary. The Frobenius nonintegrability of the 2-dimensional Pfaff systems (A,£) and (A*,() 
is a G-invariant property. 

Remark. For a possible connection of SF A F to the Godbillon-Vey closed 3-form r = d6 A 9 see 
further in sec. 11. 

Finally we note that inside an orbit G.F a limited superposition rule holds: if F is a solution and 
if [a(a, b) + (3(m, n)] 6 G then, obviosly, a(a, b).F + (3(m, n).F is also a solution. 

7.3 Point dependent group parameters 

We are going now to see what happens if the group parameters become functions of the coordinates: 
a = a(x, y, z,£)I + b(x, y, z, £)J. Consider the nonlinear solution F given by 

F = eudx A dz + epdy Adz + udx A d£ + pdy A d£, 

and the matrix 

(3 F = 0(u,p) 



u p 
—p u 

Since F is a nonlinear solution then u and p satisfy the equation 



det||/?(u,p)|| =(u 2 +p\-e(u 2 +p 2 ) z = Q, (61) 



where is the Lie derivative with respect to the intrinsically defined by F isotropic vector field Q. 
Consider now the 2-form (in the F-adapted coordinate system) 

F a = dx A C = edx A dz + dx A d£, 

which, obviously, defines a (linear constant) solution, and define a map j3 : M —>■ Q as follows: 

0(x, y, z, = u(x, y, z, £)/ + ep(x, y, z, £)J, 

The action of 

p((3(x, y, z, £)) = u(x, y, z, £)1 + ep(x, y, z, £)J 
on F gives exactly F, i.e. the (linear) solution F Q is transformed to a nonlinear solution: 

p(P(x, y, z, £))(F Q ) = F = eudx Adz + epdy Adz + udx A d^ + pdy A d£. 

This suggests to check if this is true in general, i.e. if we have a solution F defined by the two 
functions u and p, and we consider a map a : M — > G C Q, such that the components a(x,y,z,£) 
and b(x,y, z,£) of a = a(x,y, z,£)I + b(x,y, z,£)J determine a solution, then whether the 2-form 
F = a(x, y, z, £).F = [a(x, y, z,£)% + sb(x, y, z, £)J~\ .F will be a solution? 
For F we obtain 

F = p(a).F = (aX + ebJ).F = 
e(au — ebp)dx Adz + e(ap + ebu)dy Adz + (au — ebp)dx A d£ + (ap + ebu)dy A d£, 
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which looks the same as given by (60), but here a and b are functions of the coordinates. 
Now, F will be a solution iff 

[(au - ebp) 2 + (ap + ebuf] - e [(au - ebp) 2 + (ap + ebuf] z = 0. (62) 

This relation is equivalent to 

[(a 2 + b\ - e(a 2 + b 2 ) z ] (u 2 + p 2 ) + [(u 2 +p\- e(u 2 + p 2 ) z ] (a 2 + b 2 ) = 0. 

This shows that if F(u,p) is a solution, then F(u,p;a,b) = p(a(a,b)).F(u,p) will be a solution iff 
F(a,b) is a solution, i.e. iff p(a(a,b)).F is a solution. So we have the 

Proposition 14. Every nonlinear solution F(a,b) = [aX + ebJ~\F , defines a map <E>(a, b) : 
F(u,p) — ► ($F)(n,p; a, 6), such that if F(u,p) is a solution then (&F)(u,p; a,b) is also a solution. 

Further in this section we shall denote by the same letter F the solution F and the corresponding 
map $(F) : F x -> F.iq. 

The above Prop. 14 says that the set T,(F a ) of nonlinear solutions, defined by the chosen F Q , has 
a commutative group structure with group multiplication 

F(a,b).F(u,p) = p(a(a,b))F .p((3(u,p))F = p(a.(3)F = p((3.a)F = F(u,p).F(a,b). 

Clearly, F Q should be identified with p(I).F a because p(l.I + 0.eJ)(F o ) = T(F ) = F . 

The identification of a nonlinear solution F(u,p) with the corresponding map $(u,p) = p(a(u,p)) = 
ul+epj according to the rule F(u,p) = $>(u,p).F = (uX+epJ).F gives the possibility to introduce in 
an obviously invariant way the concept of amplitude of a nonlinear solution. Together with the solution 
F(u,p) we have its conjugate solution F(u, —p). We identify these two solutions with the correspond- 
ing linear maps $(u,p) and $>(u, —p) in A 2 (M) and compute the quantity \tr [<&(u,p) o <3?(u, —p)\ . 

-tr[<!>(u,p)o<f>(u,-p)] = -tr [(ul + epj) o (ul - epj)] = -tr [{u 2 + p 2 )l] = u 2 + p 2 = (f> 2 (u,p). 

In the general case a = a p .a p -\ . . . a,\ we readily obtain F(a) = F(a p ).F(a p -\) . . . F(a\), and for 
the corresponding amplitude <j)(F p .F p -i . . . F±) and phase ip(F p .F p ^i . . . Fx) we obtain 

v 

0(F p .F p _x . . . Fi)) = <t>(F p )4(F P -i) . . . ^(F P .F P ^ . . . F 1 ) = (63) 

i=l 

The corresponding scale factor is 

C{Fp ' Fp - 1 • • • Fl) = \m = \ 
Denoting £*. = C(Ff t ) we obtain the nonequality 

Ci.C2---C p 1 



C(F p .F p _ 1 ...F 1 ) > 



C 2 .C 3 ...C P + d.Cs-Cp + ... + £i.£ 2 ...£p_i £-! + £~ x + ... + C p 
If a(a, 6) 7^ at some point of M then we have F _1 = ^(q^ 1 ).^ given by 



( f )" 1 = hfi 1 -"^- 7 



Clearly i/)(F ) = 0, which corresponds to C(F a ) = 00. 
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Finally, we note that, in general, there is no superposition inside the subclass of solutions defined 
by a given solution F(u,p). The verification shows that if a(a, b) generates the solution F(a,b) and 
a(m,n) generates the solution F(m,n) inside the same (-subclass, then the sum F(a,b).F(u,p) + 
F(m,n).F(u,p) will be a solution only if the phase difference ipF(a,b) ~ ^F{m,n) ls a running wave 
along the corresponding z-coordinate. In particular, if this phase difference is a constant, then the 
superposition will work, which recalls the interference rule in Maxwell vacuum electrodynamics, that 
two waves may interfere only if their phase difference is a constant quantity. 

The considerations in this section allow the following conclusions and interpretations. The whole 
set of nonlinear solutions divides to subclasses S(F Q ) of the kind (39) (every such subclass is determined 
by the spatial direction along which the solution propagates translationally, it is the coordinate z in our 
consideration, or by the corresponding (). Every solution F of a given subclass is obtained by means of 
the action of a solution a(u,p) of equation (61) on the corresponding F Q through the representation p, 
namely, F = p(a).F Q . We also note that every such subclass divides to subsubclasses G.F, determined 
by a given solution F and the action of G with constant coefficients. Every such subsubclass may be 
considered as one solution represented in different bases of Q. 

Further, the relation p(a(u,p))F = F suggests to consider F , as a "vacuum state along £" of 
the subclass considered, and the action of p(a) on F Q as a "creation operator along provided 
a(x,y, z,£) 7^ 0, b(x, y, z, £) / 0, (x, y, z, £) € D C M, and a = al + bj satisfies (61) in D. Since 
p{ct).p{a~ l ).F = F Q we see that every such "creation operator along £" in D may be interpreted also 
as an "annihilation operator along £" in D, i.e. in the same subclass of solutions. A given a may 
generate a creation operator with respect to some subclass of solutions only if \a\ is a running wave 
along the corresponding isotropic vector (. Finally, the only universal "annihilation" operator is the 
p-image of the zero element of Q. 

Clearly, every solution F = p(a).F of a given subclass may be represented in various ways in terms 
of other solutions to (61) in the same domain D, i.e. we have an example of a nonlinear "superposition" 
inside a given subclass, represented by the above mentioned commutative group structure. 

The interpretation of a given subsubclass determined by a given solution F and the action of G 
with constant coefficients as one solution, suggests also the extension of this interpretation to the point 
dependent case. In fact, at every point m = (x, y, z,£) G M we may consider an isomorphic image Q m 
of Q with basis {ei(m), e2(m)}, where 



ei(m) 



u(m) 
u(m) 



e 2 (ra) 



p(m) 
—p(m) 



Now, the element [e\{m) + e2(m)] G Q m is sent through the representation p, m : Q m — > L A 2 i M \ to the 
linear map $(m) = u(m)X + p(m)J ', and when (u 2 +p 2 ) satisfies equation (61) we obtain the solution 
<&(F ). Hence, all solutions F(u,p) of a given subclass may be considered as one solution represented 
in different point dependant bases {ei(?7z.), 62(771)} of Qmi 

m G M. In fiber bundle terms we could say 
that the nonlinear solutions are determined by sections of the bundle M x Q which satisfy (61) and 
their p-actions on corresponding "vacuums" F . 

Looking back we see that in the frame of linear solutions (Maxwell solutions) we have additive 
commutative group structure (any linear combination with constant coefficients of linear solutions is 
again a solution), while in the frame of the nonlinear solutions we have a multiplicative commutative 
group structure: a product (as given above) of any number of nonlinear solutions along Q is again a 
solution along (. The important difference between these two structures is that the additive structure 
of the linear solutions holds for all solutions, and the multiplicative group structure holds just inside a 
given (-subclass. So, the whole set of nonlinear solutions consists of orbits of the (multiplicative) group 
of those complex valued functions a(x,y, z,£) (the product is point-wise), the module \a\ = V a 2 + b 2 
of which is a running wave along some fixed isotropic direction (. 
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We see also that the amplitude and the phase of a solution in a natural way acquire the inter- 
pretations of the amplitude <j) = \a\ = \/a 2 + b 2 and phase if) = arccos((/?) (see relation (41)) of the 
corresponding complex fields a(x,y, z,£). Clearly, through the well known Moivre formulae, every 
F = p(a).F Q generates (F) n = p(a n ).F and \[F = p(^/a).F , and in this way we obtain new 
solutions of the same subclass. 

8 A Heuristic Approach to the Vacuum Nonlinear Equations 

In this section we give a heuristic, in our view, approach as for how to come to the nonlinear vacuum 
equations of EED. Only two things turn out to be needed: first, accepting that the time dependent 
EM fields propagate along straight lines, second, some experience with differential geometry concepts. 

Our basic assumption reads: The translational component of propagation of time depen- 
dent vacuum EM fields is along isotropic straight lines in Minkowski space-time. 

We shall denote by M the Minkowski space-time, endowed with the canonical coordinates (x, y,z,£ = 
ct). From the above assumption it follows that at any particular case we may find isotropic vector 
field the integral (isotropic) lines of which can be described by the vector field £ = — ed z + or by the 
corresponding through the (flat) Lorentz metric n, sign(rj) = (—,—,—,+), 1-form ( = edz + d£. This 
1-form ( is closed, so, it defines 1-dimensional completely integrable Pfaff system. The corresponding 
3-dimensional differential system (or distribution) on M is spanned by the vector fields 

A = —ud x — pd y , A' = —md x — nd y , ( = —ed z + 

which are ^-images of the space- like 1-forms (A, A'), given by 

A = udx + pdy, A' = mdx + ndy, 

and of C = £dz + d£, where (u,p, m, n) are four independent functions on M. We form now the tensor 
Q = 4> 2 C ® C> where the function <p is still to be determined. This tensor field defines linear map in 
the spaces of vector fields and 1-forms on M. Clearly, ImQ is spanned by £, and KerQ is spanned 
by (A, A',Q, where A and A' are linearly independent. The factor space Ker Q/ImQ is isomorphic 
to the plane spanned by (A, A'), and further we choose n(A,A') = 0, i.e. A' = —pdx + udy, and we 
shall denote this choice of A' by A*. It is easily verified that the 2-dimensional Pfaff system (A, A*) 
is completely integrable, and that the relations A.( = A*.( = hold. 
We define now the 2-forms 

f = aa(, f* = a*a(. 

The following relations hold: 

i(A)F* = 0, i(A*)F = 0, i(QF = i(QF* = 0. 

Moreover, if *F is the Hodge-77 dual to F we have F* = *F. 

Obviously, F A F = and F A *F = 0, so F is a null field. Now, the tensor field 

is equal to (n 2 +p 2 )(n( I/ , so we define <j) 2 = u 2 + p 2 = —A 2 = —(A*) 2 . We compute dF and d * F: 

dF = e{p x — u y )dx A dy A dz + (p x — u y )dx Ady A d^— 
e(u^ — eu z )dx A dz A d£ — e{p^ — ep z )dy Adz A d£, 
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d * F = (p y + u x )dx A dy A dz + e(p y + u x )dx f\dy f\ d£— 

(p^ — ep z )dx A dz A d£ + (u^ — eu z )dy Adz A d£. 

The relation F A *d * F + (*F) A *dF = is now immediately verified. We compute the quantities 
i(F)dF = i(A A C)dF = o i(A)dF and i(*F)d *F = i(A* A ()d * F = oi(A*)d * F, and obtain 

i(A A ()dF = i(A* A C)d * F = \ [(u 2 + p\ - e(u 2 + p 2 ) z ] (. 

The equations 

F A *d * F + (*F) A *dF = 0, i(F)dF = 0, i(*F)d * F = 

now give ^vQ 1 ^ = and that the amplitude function ^ is a running wave along the z-coordinate. From 
these relations we can derive everything that we already know. 

The above considerations in this section made use, more or less, of the metric structure of 
Minkowski space and of the corresponding Hodge *-operator. A more careful look into the details, 
however, shows that we can come to the important relations of EED without using any metric. We 
proceed now to consider this possibility. 

Instead of Minkowski space M we consider the manifold IR 4 together with its tangent tensor and 
exterior algebra, and standard coordinates (x,y, z,£). The main tool to be used is the Poincare 
isomorphism between the p- vectors and (4 — p)-forms on M 4 [28]. For convenience we briefly recall how 
it is introduced. Let A P (V) be the space of p-tensors, i.e. all antisymmetric (contravariant) tensors, 
and A( n ~ p \V*) be the space of (n — p)-forms, i.e. all antisymmetric covariant tensors, over the pair 
of dual n-dimensional real vector spaces (V,V*), and p = 1, ...,n. If {e^} and {e 1 } are two dual 
bases we have the n-tensor to = e\ A • • • A e n and the n-form lo* = e 1 A • • • A e n . The duality requires 
< £ l ,ej >= 8j and < uj*,uj >= 1. If x G V and a £ A P (V*) then by means of the insertion operator 
i(x) we obtain the {p — l)-form i(x)a. If a is decomposable: a = a 1 A a 2 A • • • A a p , where a 1 , . . . , a p 
are 1-forms, then 

i(x)a = < a 1 , x > a 2 A ■ ■ ■ A a p - < a 2 , x > a 1 A a 3 A ■ ■ ■ A a p + . . . 
+ < a p ,x > a 1 A--- Aa p -\ 

Now let x\ A x 2 A ■ ■ ■ A x p £ A P (V). The Poincare isomorphism *p : A P (V) — > A n ~ p (V*) acts as follows: 

*P p (xi A x 2 A • • • A x p ) = i(xp) i(x p ^i) . . . Q i(x\)uj* . 

In the same way 

^(q 1 A • • • A a p ) = i{a p ) . . . Q ^a 1 )^. 
For nondecomposable tensors is extended by linearity. We note also the relations: 

qs n _ p y p = (-i) p ( n ~ p hd, y, n - p y p = {-i) pi - n ~ p \d, 

< <P P (a 4 A • • • A a p ), %$p(x 1 A ■ ■ ■ A x p ) >=< a 1 A • • • A a p , x± A ■ ■ ■ A x p > . 
On an arbitrary basis element *P acts in the following way: 

p 

q3P( e ii A • • • A = (-l)*=i e lp+1 A ■ ■ ■ A e ln , 

where i± < ■ ■ ■ < i p and i p+ \ <■■•<«„ are complementary p— and (n — p)— tuples. 
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On the manifold M 4 we have the dual bases {dx 1 } and {^r} of the two modules of 1-forms and 
vector fields. The corresponding u> and uj* are 

d d d d 
uj = — — A — A — A — , to = dx Ady Adz A <i£. 
ax ay a£ 

Our basic initial objects of interest are the 1-form ( and the vector field £ defined by 

C = e<fa + de, c = -^ + ^- 

Note that no unknown functions are introduced with these two objects. We introduce now the tensor 
field Q := (f) 2 ( © (, where the function <\> is to be determined. The tensor Q defines linear maps in the 
two spaces of vector fields and 1-forms on ]R 4 by the relations 

Q(X) = 4> 2 <(,X>(, Q*(a)=<t> 2 <a,(>(, 

where X is an arbitrary vector field, a is an arbitrary 1-form, and < a, X > is the canonical conjugation 
between TR 4 and T*R 4 . Clearly, the image space ImQ of Q is generated by £, and the image space 
ImQ* of Q* is generated by (. The kernel space KerQ is 3-dimensional and looks like: KerQ = 
{a} © {b} © {£}, where the two vectors (a, b) are linearly independent: aAb / 0. The kernel space 
KerQ* is also 3-dimensional and looks like: KerQ* = {a} © {/?} © {(}, where the two 1-forms (a,/3) 
are linearly independent: a A ft 7^ 0. The corresponding 2-dimensional factor spaces KerQ /ImQ and 
KerQ* /ImQ* naturally arise, and further we shall consider mainly the second one. The elements of 
KerQ* /ImQ* have the form a + f( and [3 + /£. Inside each of these two classes there is only one 
element that has no component along £, we denote these two elements by A and A' correspondingly. 
These two 1-forms live in the planes (x, y). So, we have A = udx + pdy, A' = mdx + ndy, with the 
relation un — pm ^ 0. 

It is now readily verified the following 

Corollary. The 2-dimensional Pfaff system (A, A') is completely integrable, and the two Pfaff 
systems (A, Q and (A, 1 Q are not, in general, completely integrable. 

Hence, the following relations hold: 

A A A A' = 0, dA' A A A A' = 0, dA A A A ( + 0, dA' A A' A C / 0. 

Intuitively, the integrable Pfaff system (A, A') carries information about those properties of the 
physical system under consideration which do not change during the evolution, and the other two 
nonintegrable Pfaff systems carry information about the admissible interaction of the physical system 
with the rest of the world. Therefore, namely these two Pfaff systems, (A, £) and (A f , £), should be 
of primary interest for us, (from the point of view of studying this physical system) because studying 
them, we may get some knowledge about the proper, or identifying, properties of the system (Recall 
Sec.l). The objects A © £ and A' © ( are also expected to be of use. 

Since A and A' are linearly independent we choose them as follows: 

A = udx + pdy, A' = —pdx + udy, so A A A' / 0. 

According to the above remark we turn our attention to the two 2-forms A A ( and A' A (, and 
our first observation is that they are defined on the classes determined by A and A' respectively. Our 
second observation is that the corresponding curvatures, given by 

dA A A A C = dA' A A' A ( = -e[p(v,£ - eu z ) - u(p^ - ep z )] to*, 
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give the same coefficient (Sec. 5) which introduces rotational component of propagation. As we shall 
see further (Sec. 9) this is not occasional. 

We compute now ^ 2 (A A £) and ^ 2 (A' A £) and obtain 



d d 



d d 



d d 



d d 



<#\A A C) = -eu— A — + ep— A — + u— A — p— A — , 



dy <9£ 
d d 



dx dy dz dx dz 



<$ 2 {A' A C) = -u^- A - p^- A + eu^- A ^ + ep^- A ^. 

ax dz ay dz ox ot, dy a£ 

The following relations are readily verified: 

i(A) [¥ 2 (A A 0] = 0, i(A') [y 2 (A' A 0] = 0, 

i(A A C) [V 2 (A A 0] = 0, i(A' A C) [¥ 2 (A A 0] = 0. 
We compute <p 3 [d(A A ()] : 



¥ 3 [d(AA0]=e(p x -u y )( + e 



(n,-eu z )^--( P ,-s Pz )^- 



The two objects A <g> ( and A' ® £ may be considered as 1-forms with values in the tangent bundle, 
so we can use the insertion operator i{X) as follows: i(X)(A (g) £) =< A,X>(. Denoting the two 
vector fields *p 3 (d(^ A C)) and ^ 3 (d(A' A C)) by Z and Z' respectively, we obtain 



i(Z)(A ® C) = ® C) = e[p(« C - eu z ) - u( Pi . - e Pz )](, 



and 



i{Z)(A' ® C) = -«(£') ® C) = " + " ep*)]C- 

Recalling now from Sec. 5 the explicit forms of SF AF = (5*F)A *F, 5F A *F and (5*F) AF, and the 
equations they satisfy, the same equations may be written as follows: the two equations SF A *F = 
and (S * F) A F = are equivalent to 

i{Z)(A'®Q = -i{Z')(A®Q =0, 

and the equation 5F A F = (6 * F) A *F is equivalent to 

i{Z){A®() = i{Z')(A' ®Q. 

The additional equation d(SF A F) = is equivalent to 



«Pi[i(Z)(A®c)] 



o. 



Finally, choosing the two vector fields a and b in the form a = ud x + p5 y and b = — pd x + ud y , the 
function <p in Q = (j) 2 ( ® ( is obtained through the relation 

Q= (<,>,<g>)(4AC,aAC) =<Aa>(®C=(ti 2 +P 2 X®C, 
or through the relation 

Q = (<,>,<g>)(A' AC,b AC) =< A',b > C® C = (« 2 +P 2 )C® C, 

The considerations made so far were limited, more or less, inside a given subclass of (nonlinear) 
solutions, which propagate translationally along the same isotropic 4-direction in M, or along a given 
spatial direction which we choose for z-coordinate. A natural question arises: is it possible to write 
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down equations which would simultaneously describe a set of N such non-interacting solutions, which 
propagate translationally along different spatial directions. The answer to this question is positive, 
and the equations look like: 

N N 

(tF™- 1 A *F 2fc - 1 ) ea.! V e 2k ^ - £ (s * F 2k ~ x A F 2 ^ 1 ) e 2k V e 2fc + 
fc=i fe=i 

N 

{-5F 2k ~ l A F 2 ^ 1 + 6 * F 2 ^ 1 A *F 2fc " 1 ) ® e 2fc _! V e 2fc = 0, 

k=l 

where the index k enumerates the 2-space (M 2 ) fc for the corresponding couple (F k , *F k ). So, for every 
k = 1,2, ...,N we obtain the system (32)-(34), i.e. the corresponding couple (e2fc-i, e2k) defines, 
in fact, the direction of translational propagation of the solution (F k ,*F k ). The case of possible 
interaction of nonlinear solutions (e.g. through "interference") will be considered elsewhere. 



9 Connection-Curvature Interpretations 

The concepts of connection and curvature have proved to be very useful in field theory, for example, 
General Relativity and the later Yang-Mills theory were built entirely in terms of these (geometric in 
nature) concepts. So, it seems interesting to try to find corresponding interpretation of EED. 

The concept of connection manifests itself in three main forms: general connection on a manifold, 
principal connection in a principal bundle, and linear connection in a vector bundle. In the first 
case the general connection gives an appropriate geometric picture of a system of partial differential 
equations, and the corresponding (general) curvature shows when this system of equations is integrable 
and when it is not integrable: the nonzero curvature forbids integrability. The principal connection 
adapts appropriately the concept of general connection to the special case of principal bundles making 
use of the rich geometric structure of Lie groups and free Lie group actions on manifolds. The linear 
connection allows to differentiate sections of vector bundles along vector fields, obtaining new sections 
of the same vector bundle. 

First we are going to introduce the necessary concepts in view of their use for corresponding 
interpretation of EED. 



9.1 Integrability and Connection 

The problem for integration of a system of partial differential equations of the kind 

^ = ft (x k , y b ), i,k = l, ...,p; a, b = 1, q, (64) 

where f?(x k ,y b ) are given functions, obeying some definite smoothness conditions, has contributed to 
the formulation of a number of concepts, which in turn have become generators of ideas and research 
directions, and most of them have shown an wide applicability in many branches of mathematics and 
mathematical physics. A particular case of the above system (nonlinear in general) of equations is 
when there is only one independent variable, i.e. when all x % are reduced to x 1 , which is usually 
denoted by t and the system acquires the form 

^- = f«(y\t), a, 6 = 1, ...,q. (65) 

We recall now some of the concepts used in considering the integrability problems for these equations, 
making use of the geometric language of manifold theory. 
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Let X be a vector field on the g-dimensional manifold M and the map c : I —> M, where I is an 
open interval in R, defines a smooth curve in M. Then if X a , a = 1,2, ... ,q, are the components of X 
with respect to the local coordinates (y 1 , ■■■,y q ) and the equality c'(t) = X(c{t)) holds for every t £ I, 
or in local coordinates, 

c(i) is the integral curve of the vector field X, i.e. c{t) defines a 1-dimensional manifold such that X is 
tangent to it at every point c(t). As it is seen, the difference with (65) is in the additional dependence 
of the right side of (65) on the independent variable t. Mathematics approaches these situations in an 
unified way as follows. The product R x M is considered and the important theorem for uniqueness 
and existence of a solution is proved: For every point p £ M and point r £ R there exist a vicinity 
U of p, a positive number e and a smooth map <I> : (r — e, r + e) x U — > M, <& : (i, y) — > <ft(y), such 
that for every point y £ U the following conditions are met: <p T (y) = y, t — > <^t(y) is an integral curve 
of X, passing through the point y € M; besides, if two such integral curves of X have at least one 
common point, they coincide. Moreover, if (t',y), (t + t',y) and (t,<p(y)) are points of a vicinity U' 
of {0} x R in R x M, we have (p t+t >(y) = tpt(<Pt' (y)) • This last relation gives the local group action: 
for every i G / we have the local diffeomorphism ip t : U — > ft(U). So, through every point of M there 
passes only one trajectory of X and in this way the manifold M is foliated to non-crossing trajectories 
- 1-dimcnsional manifolds, and these 1-dimensional manifolds define all trajectories of the defined by 
the vector field X system of ODE. This fibering of M to nonintersecting submanifolds, the union of 
which gives the whole manifold M, together with considering t and y(t) as 1-d submanifolds of the 
same manifold, is the leading idea in treating the system of partial differential equations (64), where 
the number of the independent variables is more than 1, but finite. For example, if we consider two 
vector fields on M, then through every point of M two trajectories will pass and the question: when 
a 2-dimensional surface, passing through a given point can be built, and such that the representatives 
of the two vector fields at every point of this 2-surface to be tangent to the surface, naturally arises. 
The answer to this problem is given by the following Frobenius theorems. 

For simplicity, further we consider regions of the space WP xR', but this is not essentially important 
since the Frobenius theorems are local statements, so the results will hold for any (p + g)-dimensional 
manifold. 

Let U be a region in MP xM q , and (x 1 , ...,x p ,y 1 = x p+l , ...,y q = x p+q ) are the canonical coordinates. 
We set the question: for which points (xo,yo) of U the system of equations (64) has a solution 
y a = (p a (x l ), defined for points x, sufficiently close to xq and satisfying the initial condition tp(xo) = yo? 
The answer to this question is: for this to happen it is necessary and sufficient the functions ff on 
the right hand side of (64) to satisfy the following conditions: 

Bf a 8f a df a df a 

+ *)■■/?(*> V) = -£{*,V) + ^(x,y).fKx,y). (66) 

This relation is obtained as a consequence of two basic steps: first, equalizing the mixed partial 
derivatives of y a with respect to x l and x 3 , second, replacing the obtained first derivatives of y a with 
respect to x % on the right hand side of (64) again from the system (64). If the functions /" satisfy 
the equations (66), the system (64) is called completely integrable. In order to give a coordinate free 
formulation of (66) and to introduce the curvature, as a measure for non-integrability of (64), we shall 
first sketch the necessary terminology [23]. 

Let M be an arbitrary n = p + q dimensional manifold. At every point x £ M the tangent space 
T X {M) is defined. The union of all these spaces with respect to the points of M defines the tangent 
bundle. On the other hand, the union of the co-tangent spaces T*(M) defines the co-tangent bundle. 
At every point now of M we separate a p dimensional subspace A X (M) of T X (M) in a smooth way, 
i.e. the map x — > A x ix smooth. If this is done we say that a p-dimensional distribution A on M is 
defined. From the elementary linear algebra we know that every p-dimensional subspace of T X {M) 
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defines unique (n — p) = q dimensional subspace A* of the dual to T x (M) space T* (M) , such that all 
elements of A* annihilate (i.e. send to zero) all elements of A x . In this way we get a g-dimensional 
co- distribution A* on M. We consider those vector fields, the representatives of which at every point 
are elements of the distribution A, and those 1-forms, the representatives of which at every point are 
elements of the co-distribution A*. Clearly, every system of p independent vector fields, belonging to 
A, defines A equally well, and in this case we call such a system a differential p-system V on M. The 
corresponding system V* of q independent 1-forms is called ^-dimensional Pfaff system. Clearly, if 
a G V* and X G V, then a(X) = 0. 

Similarly to the integral curves of vector fields, the concept of integral manifold of a p-dimensional 
differential system is introduced. Namely, a p-dimensional submanifold V p of M is called integral 
manifold for the p-dimensional differential system V, or for the p-dimensional distribution A, to which 
V belongs, if the tangent spaces of V p at every point coincide with the subspaces of the distribution 
A at this point. In this case V p is called also integral manifold for the g-dimensional Pfaff system 
V* ■ If through every point of M there passes an integral manifold for V , then V and V* are called 
completely integrable. 

Now we shall formulate the Frobenius theorems for integrability [25,26]. 

A differential system V is completely integrable if and only if the Lie bracket of any two vector 
fields, belonging to V , also belongs to V. 

So, if (Xi, ...,X p ) generate the completely integrable differential system V , then 

[X l ,X J ]=C%X k , (67) 

where the coefficients Cfj depend on the point. If we form the p-vector X\ A X2 A • • • A X p then the 
above relatin (67) is equivalent to 

[X i} Xj] A Xi A X 2 A • • • A X p = 0, i,j = l,2,...,p. 

We turn now to the integrability criterion for the corresponding Pfaff systems. If the 1-forms 
(a 1 ,..., a 9 ) define the g-dimensional Pfaff system V* the following integrability criterion holds (the 
dual Frobenius theorem): 

The Pfaff system V* is completely integrable if and only if 

da a = K% c a b A a c , b < c. (68) 

Now, it is easily shown, that the above equations (68) are equivalent to the following equations: 

(da a ) A a 1 A . . . A a a A ... A a 9 = 0, a = 1, q. (69) 

We note that these criteria do not depend on the special choice of the vector fields Xi, i = 1, . . . ,p, or 1- 
forms ct a , a = 1, . . . , q, where p + q = dimM , which represent the distribution V and the codistribution 
V*. 

When a given Pfaff system V* , or the corresponding differential system V , are not integrable, 
then the relations (67)-(69) are not fulfilled. From formal point of view this means that there is at 
least one couple of vector fields (X, Y), belonging to V , such that the Lie bracket [X, Y] does not 
belong to V . Therefore, if at the corresponding point x £ M we choose a basis of T X (M) such, that 
the first p basis vectors form a basis of A x , then [X, Y] will have some nonzero components with 
respect to those basis vectors, which belong to some complimentary to A^ subspace T x : T X (M) = 
AajGT^. As a rule, the distribution V is chosen to be integrable, so the choice of a projection operator 
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V x : T X (M) — > T x (M),x £ M with constant rank defines a new distribution A = Ker(V), such 
that KerV x (M) © Y X [M) = T X (M). Then the distribution V is called vertical and the distribution 
A = KeriV) is called horizontal. 

Now, the curvature K, of the horizontal distribution A is defined by 

)C(X,Y) = V([X,Y}), X,Y € A. (70) 

It is clear, that the curvature /C is a 2-form on M with values in the vertical distribution T. In fact, 
if / is a smooth function and X, Y are two horizontal vector fields, then 

K(X, fY) = V([X, fY}) = V(f[X, Y] + X(f)Y) = 

= fV([X, Y\) + X{f)V(Y) = fK(X, Y) 

because V(Y) = 0, Y — horizontal. 

The corresponding (/-dimensional co-distribution A* (or Pfaff system V*) is defined locally by the 
1-forms (9 1 , ...,9 q ), such that 9 a (X) = for all horizontal vector fields X. In this case the 1-forms 9 a 
are called vertical, and clearly, they depend on the choice of the horizontal distribution A. In these 
terms the non-integrability of A means 

d6 a / K£ c 9 b A 9 C , b< c. 

This non-equality means that at least one of the 2-forms d9 a is not vertical, i.e. it has a nonzero 
horizontal projection H*d9 a , which means that it does not annihilate all horizontal vectors. In these 
terms it is naturally to define the curvature by 

H*d9 a = d9 a - K£ c 9 b A 9 C , b< c. 

Now we see how this picture is defined by the equations (64). 

Consider now the manifold M = W x R q with coordinates (x 1 , . . . , x p , y 1 , . . . , y q ). The vertical 
distribution is defined by 

9 

— , a = p+l,...,p + q = n. 

Now, making use the equations of the system (64), i.e. the functions ff , we have to define the 
horizontal spaces at every point {x,y), i.e. local linearly independent vector fields Xi,i = 1, ...,p. The 
definition is: 

x '-t< +f 'w (71) 

The corresponding Pfaff system shall consist of 1-forms 9 a ,a = p + 1, ...,p + q and is defined by 

a = aly a - f?dx\ (72) 

In fact, 

»■<*) = *• (It) + «' - it) - ^ (/^ 

= o + /»«-/?«?-o = o. 



Remark. The coordinate 1-forms dy a are not vertical and the local vector fields are not 
horizontal with respect to the so defined horizontal distribution (X\, . . . ,X p ). 

In this way the system (64) defines unique horizontal distribution. On the other hand, if a hori- 
zontal distribution is given and the corresponding vertical Pfaff system admits at least 1 basis, then 
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this basis may be chosen of the kind (72) always. Moreover, in this kind it is unique. In fact, let 
(9 fl , ...,9' q ) be any local basis of A*. Then in the adapted coordinates we'll have 

9 ,a = A%dy b + Bfdx\ 

where A^,Bf are functions on M. We shall show that the matrix A% has non-zero determinant, i.e. 
it is non-degenerate. Assuming the opposite, we could find scalars A a , not all of which are equal to 
zero, and such that the equality XaA^ = holds. We multiply now the above equality by A a and sum 
up with respect to a. We get 

\ a e' a = \ a B a % dx\ 

Note now that on the right hand side of this last relation we have a horizontal 1-form, while on the 
left hand side we have a vertical 1-form. This is impossible by construction, so our assumption is 
not true, i.e. the inverse matrix (A^)~ l exists, so multiplying on the left 9 la by (A^)" 1 and putting 
{Al)~ 1 6' b = 9 a we obtain 

6 a = dy a + {AD^B^dx 1 . 

We denote now (A^Bf = -/? and get what we need. The uniqueness part of the assertion is proved 
as follows. Assume there is another basis {0 a )" of the same kind. So, there must be a non-degenerate 
matrix C fc a , such that {6 a )" = C%9 h . We get 

dy a - (ftY'dx* = C a h dy b - Clj\dx\ 

and from this relation it follows that the matrix C% is the identity. 

Let's see now the explicit relation between the integrability condition (66) of the system (64) and 
the curvature H*d9 a of the defined by this system horizontal distribution. We obtain 



d9 a 



dfl 
dxi 



dx j A dx l 



dfl 

8y b 



dy b A dx\ 



In order to define the horizontal projection of d9 a we note that H*{dx l ) = dx l , so, H*{dy a 
H*(9 a + ffdx 1 ) = ffdx\ because 9 a are vertical. That's why for the curvature Q, = H*d9 a <g> ^ 
get 

'df* 



we 



n = (H*d9 a ) 



d 

QyCl 



dfl 



df a 

dx i 3x3 ' dy b 1 dy bJj 



+ 



f?-^Kfi)dx l A dx> 



d 

QyCl ■ 



I < J. 



(73) 



It is clearly seen that the integrability condition (66) coincides with the requirement for zero curvature. 
The replacement of dy a , making use of the system (64), in order to obtain (66), acquires now the status 
of "horizontal projection". 

We verify now that the curvature, defined by V([Xi, Xj]) gives the same result. 



V([X i ,X j ]) = V 



dx i Jl dy a ' dxi T J J Qy b _ 



= V 







dJl_d_ 

dx % dy b dxi dy a ' J % dy a dy b 



fb dfi_d_ 

J i Q y b Q y a 



+ 



-r 



d JlAJL 

,b Jj I f)„,a ' 



dx i dxi dy b 1 dy b 3 ) dy a 

The basic moment in these considerations is that the system (64) defines a vertical projection 
operator V and a horizontal projection operator H = id — V in the tangent bundle T(M). Explicitly, 



H is defined in the basis ^) by 



H 



d _d_ 

dx l ' dy c 
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I 




\dx l 


dy aJ 


ft o 


-( 



— + r— o o 

dx l ^J* '•••' 
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where I denotes the identity (p x p)-matrix and denotes the zero (p x p) and the zero (q x q) matrices. 
So, H x projects the basis (^», ^3-) at x £ M on the horizontal subspace basis (X±(x), . . . ,X p (x)). 
The projection V x = (idx x M ~ H x ) projects every tangent space T X M to the corresponding vertical 
subspace VT X , so, the transposed matrix (V x )* = (id,T x M — H x )* projects every cotangent space to the 
vertical cotangent subspace. Explicitly we have (we omit to write the point x € M) 



(id- H)*(dx\dy a ) = (dx\dy a ) 



-f c 
I 



(0,...,0,-/<W + dy a ). 



Hence, a general connection on a manifold M is defined by a projection operator H (or V = (idxiM) ~ 
H)) in T(M) with constant rank. As a rule, if H denotes the horizontal projection as it is in our case, 
the image space of the corresponding vertical projection V = (idx(M) ~ H) is considered to define 
an integrable distribution, so the kernel subspace Ker(id T ( M ^ — H) = KerV defines the horizontal 
distribution, and the corresponding curvature carries information about the integrability of KerV. 

9.2 General Connection Interpretation of EED 

We go back now to EED. According to Proposition 8. the 1-forms (A, A*), which live in the planes 
(x,y), define a completely integrable 2-dimensional Pfaff system. We choose the vertical subspaces 
to live inside also in the (x, y)-planes, so, accordingly, the coordinates are given in a different order, 
namely: (z,£,x, y). The horizontal spaces are defined by 

„ d d d v d d d 

oz ox oy o£, ox Oy 



which corresponds to the partial differential system 

dx dx 



dy dy 



Recalling from Sec. 4 the unit 1-forms R and S, and denoting by prime (') the corresponding (through 
the metric 77) vector fields, we have the 

Corollary: The following relations hold: 

X 1 = R' + (i(R')F)', X 2 = S' + (i(S')F)'. 



The horizontal projection in TM is defined as follows (in the F-adapted coordinate system) 



-ff* 



d d d d 



dz <9£ dx dy J \dz <9£ dx dy 



d d d d 



1 














1 








eu 


u 








ep 


P 









= (X 1 ,X 2 ,0,0). 



In coordinate free form, i.e. as a section of T*(M) (g) T(M), this projection is given by 

H* = -R ® [R' + (i(R')F)'] + S®[S' + (i(S')FY] . 
The corresponding vertical 1-forms are 

6 1 = dx — eudz — ud^, 6 2 = dy — epdz — pd£, 
which corresponds to the vertical projection V* = (id — H*)*: 



(id — H*)* (dz,d£,dx,dy) = (dz,d^,dx,dy) 



—eu —ep 

— u —p 

1 

1 



(0,0,^,0 
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Respectively, for the horizontal projection of the coordinate 1-forms we obtain 



H*(dz, d£, dx, dy) 



1 














1 








eu 


u 








ep 


P 









(dz\ 

<% 

dx 

\dyj 



(dz, d£, eudz + ud£, epdz + pd£). 



Recalling our isotropic 1-form ( = edz + d£, we may write 

H*(dz)=dz, H*(d£)=d£, H*(dx) = dx -9 1 = u(, H*{dy) = dy -6 2 =p(. 
Making use of these relations we obtain 

= e(n e - eu z )dz A df, H*d9 2 = e{ Pi - ep z )dz A d£. 
Hence, for the curvature Q, determined by the nonlinear solution F(u,p) we obtain 

d d 
ft{u,p) = e(u£ — eu z )dz A d£ ® — + e(p^ - ep z )dz A d£ <g> —. 



Corollary. The nonlinear solution F(u,p) is a running wave only if the curvature Q(u,p) is equal 
to zero. 

Denoting by tj{X{) and ^(-X^) the ry-corresponding 1-forms to X\ and _X" 2 , we compute the expres- 
sions < rj(Xi), Q(u,p) > and < 77(^2 ), £l(u,p) >: 

d d 
< r?(Xi),n(u,p) >= [e(u s -eu z ) < ri(X±),— > +e(p^ - ep z ) < ^{X^,— > ]dzAd£ = 

= -\[{u 2 + p\ - e(u 2 + p 2 ) z ]dz A d£ = e < r/(X 2 ), ft(u,p) > . 

So, relations < rf(Xi),fl(u,p) >=< rj(X 2 ), >= give equation (48). 

For the square of O we obtain 

ft^f = (n 5 - e^) 2 + - e^) 2 = 2 (^ - e^ z f, /x < v. 

Recalling now the formulas in Sec. 5 we have 

Corollary: The square of the tensor field O is equal to —(5F) 2 , so \Q\ = \5F\. 

Corollary: A nonlinear solution may have rotational component of propagation only if the cur- 
vature form Q is time-like, i.e. if Q 2 > 0. 

In order to define the scale factor C in these terms we observe that the vertical components of the 
two coordinate vector fields d z and d^, given by 



V* ipz) = d z — H* (d z ) = d z — Xi = -eu d x - ep d y , 



and 



F* (<%) = - H* = <% - X 2 = -u d x - p d y , 
have the same modules equal to y/u 2 + p 2 . So, the scale factor C is given by 

r _\V.(d z )\ _ \V.(dt)\ 



\n\ 
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Another way to introduce the scale factor only in terms of and is the following. Recall that 
the space Lw of linear maps in a linear space W is isomorphic to the tensor product W* <8> W, so, the 
space Lyy* of linear maps in W* is isomorphic to W® W*. We shall define a map J~ : L\\? x L\y* — > L^w ■ 
Let <& £ and ^ G Ly^* and {e^} and {e J } be two dual bases in W and W* respectively. Then we 
have the representations 

$ = (gi ej , * = *™e m (g) e n . 
Now we define T($>, as follows: 

<jr) = jF($j e * ^ ej , ^™e m ® e n ) = < e*, e m > e n ® ej = $j*^ n ® ^ 

For the composition ir o J 7 (<1>, \E f ) we obtain 

trojF($,*) = 

We compute now the matrix A of Till*, (id — #*)*) making use of the corresponding matrix repre- 
sentations of these two linear maps given above. 



A{F{H*,(id-H*)*)) 



— eu —ep 

— u —p 

-2u 2 -2up 

-2up -2p 2 



, tr(A) = -2{u 2 +p 2 ). 



Hence, the relation L^(tr(A)) = gives equation (48): (u 2 +p 2 )% — e(u 2 +p 2 ) z = 0, and we can define 
the scale factor C as 



C = 



\\tr o F(H*, (id - 



On the other hand since C = \L^ip\ 1 (cf. relation (51)), where ip is the phase of the solution, we 
obtain that tp satisfies the differential equation 



1^1 



±\troF(H*,(id-H*)*) 



9.3 Principal Connection Interpretation 

The principal connections are general connections on principal bundles [30] appropriately adapted 
to the available bundle structure of the manifold. A principal bundle (V, tt, B, G) consists of two 
manifolds: the total space V and the base space B, a Lie group G, and a surjective smooth map 
7r : V — ► B. There is a free smooth right action R : V x G — > V of G on V, such, that the fibers 
tt^ 1 (x),x G B, are orbits G x of G through x G B, and every fiber G x is diffeomorphic to G. Since 
locally V is diffeomorphic to U a x G, where U a C B is an open subset of B, the right action of 
G on "P is made appropriately consistent with the group product G x G — ► G, namely, there is a 
coordinate representation (U a ,ip a ), where ip a : U a x G — > 7r _1 (L 7 Q ) is a diffeomorphism such, that 
R(tp a (x, a), b) = ip a (x, a).b = ip a (x, ab); a,b £ G. 

If is the Lie algebra of G then the right action R carries the elements h £ Q to fundamental 
(vertical) vector fields on V, and Z^, h G G, define integrable distribution Vg on V. So, choosing 
additional to Vg distribution H, we obtain a general connection on V . To make this general connection 
a principal connection we make it consistent with the right action R of G on V in a natural way: the 
subspace H za at the point z.a G "P coincides with the image of H z through (dR a ) z , i.e. H za = 
(dR a ) z (H z ). The projection on the horizontal distribution is denoted by H*. 

The horizontal distribution can be defined by a connection form co on V, which is (/-valued 1-form 
and satisfies the conditions: u)(Zh) = h, h G Q; oo o H* = and R*uj = Ad(a" 1 ) ou, a G G. Then the 
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curvature of the connection lu is given by f2 = da; o H m = du> + ^[uj, u>\. If the group is abelian, as 
it is our case, then = duj. 

If the bundle is trivial, i.e. V = M x G then the projection ir is the projection on the first member: 
ir(x, a) = x,x G M,a G G. In this case with every connection form uj can be associated a ^-valued 
1-form 9 on the base space such, that ui(x, e; X, Z^) = h + 6(x; X), where x G M, X G T X M, h G Q and 
e is the identity of G [30, p. 290]. For the curvature in the abelian case we obtain 17 = Ti*d9. 

We go back now to EED. The base manifold of our trivial principal bundle V = M x G is, of 
course, the Minkowski space-time M, and the abelian Lie group G is given by the (2 x 2)-real matrices 

u 2 + p 2 / 0. 

The corresponding Lie algebra Q, considered as a vecttor space, has the natural basis (/, J) and, as 
a set, it differs from G just by adding the zero (2 x 2)-matrix. According to the above, in order to 
introduce a connection on V = M x G it is sufficient to have a ^-valued 1-form on M. We recall now 
that a nonlinear solution F in the F-adapted coordinate system defines the 1-form C = edz + and 
the two functions (u,p). So, we define 9 as 

9 = u( <g> / + p( <8> J. 

For the curvature we obtain 

n = 9, 1 ®I + {l 2 ®J = d6= [u x dx A C + Uydy A £ — e(u^ - eu z )dz Ad(\ ® I 

+ [p x dx A C + Pydy AC — £{P( — sp z )dz A d£] J. 

Since is horizontal and, in our case, it is also G-invariant, there is just one (/-valued 2-form $ on M 
such, that 7r*<I> = fl Further we identify $ and Q and write just since in our coordinates they look 
the same. So, we can find *S7 with respect to the Minkowski metric in M: 

*S7 = (*&) ® I + (*^ 2 ) ® J = e [u x dy A ( — u y dx A £ - s(v,£ - eu z )dx A dy] ^ I 

+e [p x dy A£— p y dx A C — e(j>£ — ep z )dx A ciy] ® J. 

We recall that the canonical conjugation a — > a* in £/, given by (I, J) — > (I, — J), defines the inner 
product in Q by < a,/3 >= \tr(ao(3*). We have < 1,1 >= 1, < J, J >= 1, < /, J >= 0. We compute 
the expressions * < £1, *f2 >, < 9, *£l >, A(9, and obtain respectively: 

* < ft, *tt >= A < I, I > + * (Q 2 A *S1 2 ) < J, J > 

= - eu z f + - eu z f = -{5Ff. 
< 9, *Q >= — [u(ii(: — eu z ) + p(p(. — ep z )] dx A dy A dz 
-e[u(u£ - eu z ) +p(p£ - ep z )]dx A dy A d£, 

A(9, = | [p(ii£ — eu z ) — u(p£ — ep z )] dx Ady A dz 
+e [p(ii£ — eu z ) — u(p^ — ep z )] dx A dy A <iC j ® I A J = 5F A F ® / A J. 

Corollary. Equations (18)-(20) are equivalent to the equation < 9, *Q >= 0; a non-linear solution 
may have rotational component of propagation only if < $7, *Q 0; the equation d(SF A F) = is 
equivalent to d[ A (9, = 0. 

Since our principal bundle is trivial, the fundamental vector fields Zi and Zj on M x G, generated 
by the basis vectors (J, J) in Q coincide with the corresponding left invariant vector fields on G. Since 



a(u,p) 



u p 
—p u 
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we have a linear group action, a left invariant vector field Z a on G, at the point a 6 G, generated by 
a G is given by Z a = (a, a o a). So, in our case, at the point (x,[3 = (u,p)) G M x G we obtain 
Zf = ul + pj; Zj = —pl + uj. Hence, the metric on G gives = |Zj| = yju 2 +p 2 . Therefore, for 
the scale factor C we readily obtain 

\Zi\ \Zj\ 



c = 



Corollary: The 1-form 9, the curvature 2-form O and the scale factor C are invariant with respect 
to the group action (M x G, G) -> (M x G). 

Finally we note that we could write the 1-form 9 in the form 

9 = uC,®I + pC®J = C,®uI + C®pJ = t®Z I . 
If we start with the new 1-form 

9' = C <8> Zj = ( <g> (-p/) + C <8> (uJ) = -pC ® I + u( <8> J 
then, denoting d9' = Q', we obtain 

* < n',*fi' >= * < n,*n >, < 9',*n' >= - < o,*n >, a(6»',*o / ) = - a(o,*q). 

Hence the last corollary and the definition of C stay in force with respect to 9' too. 

Remark. The group G acts on the right on the basis (I, J) of Q. So, the transformed basis with 
f3(u,ep) G G is Rp{I, J) = (ul — epJ, epl + uJ). Now let u> be 1-form on M such that uj 2 < and 
A *C = 0. Then the corresponding generalized field is given by 

F u I + *F W J == {uj A C) <8> (ul - epJ) + *(w AC)® (ep/ + nJ). 

In particular, u> = dx defines the field in an F-adapted coordinate system. 

9.4 Linear connection interpretation 

Linear connections V are lst-order differential operators in vector bundles. If such a connection V is 
given and a is a section of the bundle, then Va is 1-form on the base space valued in the space of 
sections of the vector bundle, so if X is a vector field on the base space then i(X)Va = Vja is a new 
section of the same bundle. If / is a smooth function on the base space then V(/cr) = df (g) a + /Vcr, 
which justifies the differential operator nature of V: the components of a are differentiated and the 
basis vectors are linearly transformed. 

Let e a and e b ,a,b = 1,2, ... ,r be two dual local bases of the corresponding spaces of sections: 
< e b ,e a >= 5 b a , then we can write 

a = a a e a , V = d <g> id + T^dx* (g> (e a <g> e b ), V(e n ) = T^Jx^ <g> e b , 

so 



V(a m e m ) = da m ®e m + o m T b dx» < e a , e m >®e b 



and r^ a are the components of V with respect to the coordinates {x^} on the base space and with 
respect to the bases {e a } and {e b }. Since the elements (e a (g> e b ) define a basis of the space of (local) 
linear maps of the local sections, it becomes clear that in order to define locally a linear connection it 
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is sufficient to have some 1-form 9 on the base space and a linear map 4> = <P b a £a ® e f> m the space of 
sections. Then 

V(cr) = dcj a © e a + (9 © </>(ct) 

defines a linear connection with components = # M </>k in these bases. So, locally, a linear connection 
V may be written as 

V = d © (e a © e a ) + ^dx^ © (e a © e 6 ). 

If *i and *2 are two such 1-forms then a map (^i, ^2) — ► (A, ©)(*i, ^2) is defined by (we shall write 
just © for (A, ©) and the usual o will mean just composition) 

©(*!, * 2 ) = ^i)U^mdx» A dx v © [ o (e a © e „, e m © e n )] 

= (^i)U^ m dx^ A dx v © [ < e a , e n > (e m © e b )] = ^ l ) b ^ 2 ) a um dx^ A dx" © (e m © e b ), \i < v. 
Now, in the case of trivial vector bundles, the curvature of V is given by [30, p. 328] 

d(^ a dx^) ©(e a ©e a )+®(^,^). 

We go back now to EED. The vector bundle under consideration is the (trivial) bundle A 2 (M) of 
2-forms on the Minkowski space-time M. From Section 7 we know that a representation p' of Q is 
defined by relation (58). So, if a(u,p) £ Q then p'(a) is a linear map in A 2 (M), and recalling our 
1-form ( = edz + d£ we define a linear connection V in A 2 (M) by 

V = d © id A 2( M ) + C © p'(a(u,p)) = d © id A 2( M ^ + C © (ul + pj), a e G- 

Two other connections V and V* are defined by 

p'(a(u,p)) = p'(a(u,-p)) =ul-pj, and (p')*(a(u,p)) = p'(a.J) = p'(a(-p,u)) = -pl + uj, 

and we shall denote: 

X = ul + pj, x = ul - pj, x* = -pi + uj. 

So, we may say that x represents F and x* represents *F. 

Denoting * = ( © x, * = C ® X> ^* = C © X*> w e note that (since C A ( = 0) 

©(#, *) = ©(*, *) = ©(*, **) = 0. 

Now, since \& = u(®T+p(®J', ^> = u£© J — p(,®J, ^* = — pC^I+uC^J for the corresponding 
curvatures we obtain 

1Z = d(u() © J + d(pQ © J, K = d«) © J - d(pC) © J, =d(-K)0l + d(<) 

Remark. We have omitted here £ in front of j\7, but this is not essential since, putting p — > in 
the expressions obtained, we easily restore the desired generality. 

By direct calculation we obtain: 

*\ Tr [®(*,*d*)] = -e[u(u£ -eu z ) + p(p£ - £p z )\dz - [u(u^-eu z ) + p(p^ - £p z )]d£; 

\tt [©(**, *d\I/)] = £ p(ii£ — eu z ) — u(p£ — ep z ) dx A dy A dz+ 
p(v,£ — eu z ) — u(p£ — ep^ dx Ady A d£ = 5F A F; 
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Denoting by \TZ\ 2 the quantity \ \ * TV [®{TZ A *1Z)] | we obtain 

\1Z\ 2 = i| * Tr [®(d% *d*)] | = K - en z ) 2 + - ep z ) 2 = |,5F| 2 . 
Finally, since in our coordinates 

^tr(x) = ^tr(ul + pj) = u, and ^tr[(x o x)] = ^tr[(ul + pj) o (ul - pj)] = u 2 + p 2 , 
for the phase i/j and for the scale factor C we obtain respectively 



\tr X r y/h tr (x°x) y/tr( X °X) 
ip = arccos — ; =, L = -= = — . 

hr( X ox) J\\n\ \ n \ 



61 

These results allow to say that choosing such a linear connection in A 2 (M) then our nonlinear equations 
are given by Tr *d\£)] = 0, and that the non-zero value of the squared curvature invariant \TZ\ 2 

guarantees availability of rotational component of propagation. 



As a brief comment to the three connection-curvature interpretations of the basic relations of EED 
we would like to especially note the basic role of the isotropic 1-form £: in the general connection 
approach ( defines the horizontal projection of the coordinate 1-forms dx,dy: H*(dx) = u( and 
H*(dy) = p(; in the principal connection approach ( defines the 1-form 6 = ( ® Zi, or 9' = Q <g) Zj\ 
in the linear connection approach it defines the components of V through £ ® p'(a(u,p)). It also 
participates in defining the 2-form F Q = dx (g> £, which gives the possibility to identify a nonlinear 
solution F(u,p) with an appropriately defined linear map p'(a(u,p)) =ul + pj in A 2 (M). 

This special importance of £ is based on the fact that it defines the direction of translational 
propagation of the solution, and its uniqueness is determined by the EED equations (18)-(20): zero 
invariants F llv F tlv = F^ U {*F)^ U = are required for all nonlinear solutions. 

For all nonlinear solutions we have SF ^ 0, and all finite nonlinear solutions have finite amplitude 
(f>: < 4> 2 = \tr{F o F) = \tr(x ° x) < °°- The scale factor C separates the finite nonlinear solutions 
to two subclasses: if C — ► oo, i.e. \SF\ = = \TZ\ = 0, the solution has no spin properties; if C < oo, 
i.e. \5F\ = |f2| = \1Z\ ^ 0, the solution carries spin momentum. 

Hence, in terms of curvature we can say that the nonzero curvature invariant or \1Z\, is 
responsible for availability of rotational component of propagation, in other words, the spin properties 
of a nonlinear solution require non-zero curvature. 



10 Nonlinear Solutions with Intrinsic Rotation 

Before to start with spin-carrying solutions we briefly comment the nonlinear solutions with running 
wave character. These solutions require \5F\ = \S * F\ = 0, so, for the two spatially finite functions u 
and p we get in an F-adapted coordinate system u = u(x, y, £ + ez) and p = p(x, y, £ + ez). Whatever 
the spatial shape and spatial structure of these two finite functions could be the whole solution will 
propagate only translationally along the coordinate z with the velocity of light c without changing its 
shape and structure. In this sense this class of nonlinear solutions show soliton-like behavior: finite 
3d spatial formations propagate translationally in vacuum. If we forget about the spin properties of 
electromagnetic radiation, we can consider such solutions as mathematical models of classical finite 
electromagnetic macro-formations of any shape and structure, radiated by ideal parabolic antennas. 
Maxwell equations can NOT give such solutions. 

Now we turn to spin-carrying solutions. The crucial moment here is to find reasonable additional 
conditions for the phase function tp, or for the phase tp = arccos ip. 
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10.1 The Basic Example 

The reasoning here follows the idea that these additional conditions have to express some internal 
consistency among the various characteristics of the solution. A suggestion what kind of internal 
consistency to use comes from the observation that the amplitude function (p is a first integral of the 
vector field £, i.e. 

In order to extend this consistency between £ and </> we require the phase function p> to be first integral 
of some of the available F-generated vector fields. Explicitly, we require the following: 

The phase function p is a first integral of the three vector fields A, A* and S: 

Afa) = A» = S(p) = 0, 
and the scale factor £ is a first integral o/R : R(£) = 0. 

The first two requirements A(ip) = A*(p) = define the following system of differential equations 
for p: 



di P A 2 di P n R 2 di P di P n 

-fjr- - VI - <P = 0, a/1 - - Vjt = °- 

ox ay ox oy 



Noticing that the matrix 



-If -y/i - ip 



ip 2 -ip 

has non-zero determinant, we conclude that the only solution of the above system is the zero-solution: 

dp = dp = Q 
dx dy 

We conclude that in the coordinates used the phase function p may depend only on (z,£). The third 
equation S(p) = requires <p not to depend on £ in this coordinate system, so, p = p(z). For C we 
get 

l^zl 

Now, the last requirement, which in these coordinates reads 

\ip z \ 

means that the scale factor C is a pure constant: C = const. In this way the defining relation for C 
turns into a differential equation for p: 

(74) 

\p z \ oz £ 

The obvious solution to this equation is 

<p(z) = cos ^k— + const^j , (75) 

where k = ±1. We note that the naturally arising in this case spatial periodicity 27r£ and characteristic 
frequency v = c/2ir£ have nothing to do with the corresponding concepts in CED. In fact, the quantity 
C can not be defined in Maxwell's theory. 
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The above considerations may be slightly extended and put in terms of the phase ip, and in these 
terms they look simpler. In fact, we have the equation 

^-e^z = «^) k = ±1, 

where C = const. So, we get the two basic solutions 

z £ 
ipi = k— + const, tp2 = k— + const. 

We get two kinds of periodicity: spatial periodicity along the coordinate z and time-periodicity 
along the time coordinate £. The two values of k = ±1 determine the two possible rotational structures: 
left-handed (left polarized), and right-handed (right polarized). Further we are going to concentrate 
on the spatial periodicity because it is strongly connected with the spatial shape of the solution. In 
particular, it suggests to localize the amplitude function cj) inside a helical cylinder of height 2nC, so, 
the solution will propagate along the prolongation of this finite initial helical cylinder in such a way 
that all points of the spatial support shall follow their own helical trajectories without crossings. For 
such solutions with \SF\ ^ we are going to consider various ways for quantitative description of 
the available intrinsic rotational momentum, or the spin momentum, of these solutions. We call it 
spin-momentum by obvious reasons: it is of intrinsic nature and does not depend on any external point 
or axis as it is the case of angular momentum. And it seems quite natural to follow the suggestions 
coming from the corresponding considerations in Sec. 5. 



10.2 The ^-Approach 

In this first approach we make use of the corresponding scale factor C = const, of the isotropic 1-form 
C and of the two objects Zj = ul + pJ and Zj = —pi + uJ, considered as (/-valued functions on M. 
By these quantities we build the following Q A (/-valued 1-form H: 

H = 2itk^(®(Z i AZj). (76) 



In components we have 
In our system of coordinates we get 



H? = 27r K ^UZfZ b J -Z b I Z a J ). 



£ ,2 



H = 2iTK-(f) z (edz + d£) <g> I A J, 



hence, the only non-zero components are 



H\ 2 = 2iTKe-<j) 2 , H} 2 = 2irK-(f) 2 . 

c c 

It is easily seen that *H is closed: d * H = 0. In fact, 

r 

d * H = 2ttk- [(4> 2 )^ - £(4> 2 ) z ] (dxAdyAdA d£) ® (I A J) = 
because 4> 2 is a running wave along the coordinate z. We restrict now *H to M 3 and obtain 

(*H) R 3 = 2iTK—(f) 2 (dx Ady A dz) ® (I A J). 
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According to Stokes theorem, for finite solutions, we obtain the finite (conserved) quantity 

H = I *H R3 = 2ttk-E = kET I A J, 

JR3 C 

which is a volume form in Q, T = 2tt£/c, and E is the integral energy of the solution. The module 
|H| of H is |H| = ET. 

We see the basic role of the two features of the solutions: their spatially finite/concentrated 
nature, giving finite value of all spatial integrals, and their translational-rotational dynamical nature 
with \SF\ / 0, allowing finite value of the scale factor C. 



10.3 The FN-Bracket Approach 

We proceed to the second approach to introduce spin-momentum. We recall the Froliher-Nijenhuis 
bracket Sf of the finite nonlinear solution F, evaluated on the two unit vector fields A and eA*: 

(S F )%A' i eA* v = (S F r i2 (^eA* 2 - A^A* 1 ). 

For (Sf)i2 we g e t 

(S F )\ 2 = (S F ) 2 12 = 0, (S F f 12 = -e(S F )i 2 = 2e{ P { H - eu z ) - u( Pi - e Pz )}. 

It is easily seen that the following relation holds: A 1 eA* 2 — A 2 eA* 1 = 1. Now, for the above obtained 
solution for ip we have 

u = 4>{x, y, £ + ez) cos ^k— + constj , p = (p(x, y, £ + ez) sin ^k— + const^j . 

We obtain 

(S F f 12 = -e(S F )i 2 = "2^0 2 , 



(S F )°A»eA* 



0,0,-2e-^,2-( 



Since (f> 2 is a running wave along the z-coordinate, the vector field S F (A,eA*) has zero divergence: 
V„ [S F {A, eA*)]^ = 0. Now, defining the helicity vector of the solution F by 

X F = ^S F (A,eA*), 



then T,p has also zero divergence, and the integral quantity 

{T l p) i dxdydz 
does not depend on time and is equal to kET. 



10.4 The d(F ASF) = Approach 

Here we make use of the equation d(F A SF) = and see what restrictions this equation imposes on 
ip. In our system of coordinates this equation is reduced to 

d(F A SF) = s$ 2 + Vzz - 2eV z f) dx A dy A dz A d£ = 0, (77) 

i.e. 

% + *l>zz - tel>zt = ~ e ^) 5 - £ - e ^) z = 0. (78) 
Equation (78) has the following solutions: 
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1°. Running wave solutions ip = ip(x, y, £ + ez), 
2°. ip = £.g(x, y,i + ez) + b(x, y), 
3°. V = z.g(x, y,£ + ez) + b(x.y), 

4°. Any linear combination of the above solutions with coefficients which are allowed to depend 
on (x,y). 

The functions g{x, y, £ + ez) and b{x, y) are arbitrary in the above expressions. 

The running wave solutions tp±, defined by 1°, lead to F A 5F = and to \SF\ = 0, and by this 
reason they have to be ignored. The solutions %p2 and 1/13, defined respectively by 2° and 3°, give the 
scale factors C = l/\g\, and since C is dually invariant it should not depend on (x, y) in this coordinate 
system. Hence, we obtain g = g(£ + ez), so, the most natural choice seems g = const, which implies 
also C = const. A possible dependence of ip on (x,y) may come only through b(x,y). Note that the 
physical dimension of C is length and b(x, y) is dimensionless. 

We turn now to the integral spin-momentum computation. In this approach its density is given 
by the correspondingly normalized 3-form F A 8F. So we normalize it as follows: 

(3 = 2vr— F A SF = 2vr— [-e(f) 2 {^ - eip z )dx Ady Adz- 4> 2 {4>^ - eip z )dx A dy A d£] . (79) 

The physical dimension of (3 is "energy-density x time". We see that (3 is closed: d/3 = 0, so we may 
use the Stokes' theorem. The restriction of (3 to IR 3 is: 

£2 

(3^3 = 2tt — [— £0 2 (V ; g — eip z )dx Ady A dz] . 

Let's consider first the solutions 3° of equation (78) with C = const. The corresponding phase 
ip = K-j* + b(x,y),K = ±1, requires spatial periodicity of 2ir£ along the coordinate z. So, if we 
restrict the spatial extension of the solution along z to one such period l = 2irC, our solution 
will occupy at every moment a one-step part of a helical (screw) cylinder. Its time evolution will 
be a translational-rotational propagation along this helical cylinder. So, we have an example of an 
object with helical spatial structure and with intrinsical rotational component of propagation, and 
this rotational component of propagation does NOT come from a rotation of the object as a whole 
around some axis. 

On the contrary, the solutions defined by 2°, are NOT obliged to have spatial periodicity. Their 
evolution includes z-translation and rotation around the z-axis as a whole together with some rotation 
around themselves. 

For the case 3° with C = const we can integrate 

2ir£ 2 

Ar 3 = K( P dx Ady A dz 

c 

over the 3-space and obtain 

f 9tt C 

I 13 = kE = kET = ±ET, (80) 

where E is the integral energy of the solution, T = 2nC/c is the intrinsically defined time-period, 
and k = ±1 accounts for the two polarizations. According to our interpretation this is the integral 
spin-momentum of the solution for one period T. 

10.5 The Nonintegrability Approach 

Here we make use of the observation that the two Pfaff systems (A,() and (A*,Q are nonintegrable 
when C / 0. From Sec. 5 we have 

dA A A A C = dA* A A* A ( = e [u(p^ - ep z ) - p(u^ - eu z )] dx Ady Adz A d£ 

= e(t> 2 (ip^ — eip z )dx A dy A dz A d£. 
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Integrating the 4-form 
on the 4-volume M 3 x C we obtain kET. 



We recall also that for finite solutions the electromagnetic volume form lu x = —-A A sA* A R A S 
gives the same quantity ET when integrated over the 4-volume M 3 x £, C = const. 

10.6 The Godbillon-Vey 3-form as a Conservative Quantity 

According to the Frobenius integrability theorems having a completely integrable (n — l)-dimensional 
differential system on a n-manifold M is equivalent to having a suitable completely integrable 1- 
dimensional Pfaff system on the same manifold. This Pfaff system is determined by a suitable 1-form 
u, defined up to a nonvanishing function: fco, fix) /0,x£ M, and u) satisfies the equation duo Auj = 
(so obviously, fco also satisfies d(fuj) A fu = 0). From this last equation it follows that there is 1-form 
9 such, that duo = 9 A oj. Now, the Godbillon-Vey theorem [31] says that the 3-form (3 = d9 A 9 is 
closed: d(3 = d(d9 A 9) = 0, and, varying 9 and to in an admissible way: 9 — > (9 + gcu); to — > /u;, leads 
to adding an exact 3-form to /?, so we have a cohomological class T defined entirely by the integrable 1- 
dimensional Pfaff system. From physical point of view the conclusion is that each completely integrable 
1-dimensional Pfaff system on Minkowski space generates a conservation law through the restriction 
of (3 on R 3 . 

Recall now the following objects on our Minkowski space-time: A, A* and £. These are 1-forms. 
We form the corresponding vector fields through the Lorentz-metric and denote them by A, A*,(,. 
Let's consider the 1-form u = f( = efdz + fd£, where / is a nonvanishing function on M. We have 
the relations: 

u(A) = 0, lo{A*) = 0, lu{() = 0. 
Moreover, since ( is closed, u = f( satisfies the Frobenius integrability condition: 

dco A u = fdf A C A C = 0. 

Therefore, the corresponding 1-dimensional Pfaff system, defined by u = /£, is completely integrable, 
and there exists a new 1-form 9, such that duo = 9 A u, and d(d9 A 9) = 0. 

From the point of view of generating a conservative quantity through integrating the restriction 
i* (5 of P = d9 A 9 to M 3 , through the imbedding i : (x, y, z) — > (x, y, z, 0) it is not so important whether 
r is trivial or nontrivial. The important point is the 3-form (5 to have appropriate component /3i23 
in front of the basis element dx A dy A dz, because only this component survives after the restriction 
considered is performed, which formally means that we put d£ = in (3. The value of the corresponding 
conservative quantity will be found provided the integration can be carried out successfully, i.e. when 
(«*/3)i23 has no singularities and (i*/3)i23 is concentrated in a finite 3d subregion of ]R 3 . 

In order to find appropriate 9 in our case we are going to take advantage of the freedom we have 
when choosing 9: the 1-form 9 is defined up to adding to it an 1-form 7 = gto, where g is an arbitrary 
function on M, because 9 is defined by the relation du = 9 A u>, and [9 + guj) Aw = 9 Au always. The 
freedom in choosing to consists in choosing the function /, and we shall show that / may be chosen in 
such a way: u = f £, that the corresponding integral of i*f3 to present a finite conservative quantity. 

Recalling that d( = 0, we have 

du = d(/C) = df A C + fd( = df A C- 
Since du must be equal to 9 A lo we obtain 

duJ = dfAC = g AuJ = A ( /C ) = /0 A C- 
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It follows 



9 A C = jdf A C = d(ln /) A C = [d(ln f) + K 



AC, 



where h is an arbitrary function. Hence, in general, we obtain 9 = d(ln f) + hQ. Therefore, since 
d9 = dh A £ for d9 A 9 we obtain 

d# A0 = d(lnf) A d/i A C- 

Denoting for convenience (Zn f) = (p for the restriction i*/3 we obtain 

z*/3 = e(ip x h y — ip y h x )dx A dy A dz. 

In order to find appropriate interpretation of i* (5 we recall that *(5F A F) = — £(j) 2 (ifj^ — £ip z )£, so, 
*(<5-F A -F) is of the kind / £, and it defines the same 1-dimensional Pfaff system as / £ does. We recall 
also that if the scale factor C = 1/|^ — eip z \ is a nonzero constant then 4> 2 {^^ — £ip z ) is a running wave 
and the 3-form 5F AF is closed. Hence, the interpretation of i*(3 as i*(5F AF) = — jr£K(j) 2 dx Ady Adz 
requires appropriate definition of the two functions / and h. So we must have 



ip x hy — h x (p y 



K 



If we choose 



/ = exp(ip) = exp 



(ft 2 dx 



C 



y + const 



b 2 ) y and h xy 



h 



l)X 



0. 



all requirements will be fulfilled, in particular, ip xy = (p yx 

Hence, the above choice of / and h allows the spatial restriction of the Godbillon-Vey 3-form (3 to 
be interpreted as the spatial restriction of F A 5F. So, the curvature expressions found in the previous 
sections, as well as the corresponding spin-properties of the nonlinear solutions being available when 
SF A F ^ 0, are being connected with the integrability of the Pfaff system u = f (. 

On the two figures below are given two theoretical examples with k = — 1 and k = 1 respectively, 
amplitude function <fi located inside a one-step helical cylinder with height of 2nC, and phase function 
cp = cos(kz/£). The solutions propagate left-to-right along the coordinate z. 




Figure 1: Theoretical example with k = — 1. The Poynting vector is directed left-to-right. 
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In the case n = — 1 the magnetic vector is always directed to the rotation axis, i.e. it is normal 
to the rotation, and the electric vector is always tangent to the rotation, so, looking from behind (i.e. 
along the Poynting vector) we find clock-wise rotation. In the case k = 1 the two vectors exchange 
their roles and, looking from behind again, we find anti-clock-wise rotation. From structural point of 
view the case k = — 1 is obtained from the case k = 1 through rotating the couple (E, B) anti-clock- 
wise to the angle of tt/2, hence we get the (dual) transformation (E, B) — > (— B,E). The dynamical 
roles of the two vectors are exchanged: now the electric vector drags the points of the object towards 
rotation axis, the magnetic vector generates rotation. In both cases the Poynting vector "pushes" the 
object along the rotation axis. 

The above pictures suggest the interpretation that the rotation axis directed vectors keep the object 
from falling apart, and such kind of dynamical structure seems to be the only possible time-stable 
one. 

If we project the object on the plane orthogonal to the Poynting vector we shall obtain a sector 
between two circles with the same center, and this sector has nontrivial topology. One of the two 
vectors is always directed to the center of the circles and this stabalizes the solution, and the other is 
tangent to the circles and correspondingly oriented. The stability of the construction is in accordance 
also with the fact that, when z runs from zero to 2ir£ any of the two vectors performs just one full 
rotation. Any other rotaional evolution would make the center-directed vector leave its directional be- 
havior and this would bring to structural changes and, most probably, to falling apart of the structure. 
Choosing orientation and computing the corresponding rotation numbers of E and B we shall obtain, 
depending on the orientation chosen, (+1), or (—1), and in definite sense these values guarantee from 
mathematical viewpoint the dynamical stability of the solution-object. 

11 Explicit Non- vacuum Solutions 

Before to start searching solutions with nonzero a l ,i = 1,2,3,4, we must keep in mind that the field 
F now will certainly be of quite different nature compare to the vacuum case, and its interpretation 
as electromagnetic field is much conditional: no more propagation with the Lorentz invariant speed, 
no more photon-likeness, etc. The solution is meant to represent a field interacting continuously 
with other physical system, so the situation is quite different and, correspondingly, the properties of 
the solution may differ drastically from the vacuum solutions properties. For example, contrary to 
the nonlinear vacuum case where the solutions propagate translationally with the velocity of light, 
here it is not excluded to find solutions which do not propagate with respect to an appropriate 
Lorentz frame. In integrability-nonintegrability terms this would mean that at least some of the 
proper integrability properties of the vacuum solutions have been lost. On the other hand, the proper 
integrability properties of the external field are, at least partially, guaranteed to hold through the 
requirement for the Frobenius integrability of the 2-dimensional Pfaff systems {a 1 ,a 3 ). 

As a first step we expect possible physical interpretation just in the frame of macro phenomena, 
where the continuity of the energy-momentum exchange is easily understood, but we do not exclude 
interpretations in the frame of microworld in view of the large quantity of variables we have (in fact, 22 
functions), and in view of the universality of the physical process described, namely, energy- momentum 
exchange. 

From purely formal point of view finding a solution, whatever it is, legitimizes the equations con- 
sidered as a consistent system. Our purpose in looking for solutions in the nonvacuum case, however, 
is not purely formal, we'd like to consider the corresponding solutions as physically meaningful, in 
other words, we are interested in solutions, which can be, more or less, physically interpretable, i.e. 
presenting more or less approximate models of real objects and processes. That's why we'll try to meet 
the following. First, the solutions must be somehow physically clear, which means that the anzatz 
assumed should be comparatively simple and its choice should be made on the base of a preliminary 
analysis of the physical situation in view of the mathematical model used. Second, it is absolutely 
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obligatory the solutions to have well defined local and integral energy and momentum. Third, the 
solutions are desirable to be in the spirit of the soliton-like consideration of the real natural objects 
when a corresponding physical interpretation is meant to be done. Finally, it would be nice, the solu- 
tions found to be comparatively simple and to have something to do with "popular" and well known 
solutions of "well liked" equations. 

11.1 Choice of the anzatz and finding the solutions 

Let's now get started. We shall be interested in time-dependent solutions, therefore, the "electric" 
and the "magnetic" components should present. The simplest (F, *F), meeting this requirement, look 
as follows (we use the above assumed notations): 

F = -udy Adz - vdy A <2£, *F = vdx A dz + udx A d£, (81) 

where u and v are two functions on M. Second, when choosing the 1-forms a 1 we shall follow the 
requirement, that the "medium" does not influence the F <-> *F exchange, so we put 

a 2 = a 3 = 0. (82) 

Third, since the rest two 1-forms , a 1 and a 4 must be linearly independent (they must define a 2- 
dimensional Pfaff system), we choose them to be mutually orthogonal and let a 4 be time-like, and the 
other to be space-like. The simplest time-like 1-form looks like a = A(x, y, z, £,)d£,, moreover, any such 
1-form defines integrable 1-dimensional Pfaff system: 

da A a = {A x dx Ad£ + A y dy + A z dz A d£) A Ad£ = 0. 

So, we choose a 4 = A{x,y,z,^)d^. Finally, for a 1 , denoted further by (3, we obtain in general, 
(3 = (3\dx + fody + fodz. Clearly, (3 2 < 0, and a. [3 = 0, so they are linearly independent. Now, j3 
participates in the equations through the expression f3 A *F, and from the explicit form of *F is seen, 
that the coefficient f3\ in front of dx does not take part in /3A*F, therefore, we put (3\ = 0. So, putting 
02 = b and = — B for convenience, we get 

a 1 = (3 = bdy- Bdz, a 4 = a = Ad£, a 2 = a 3 = 0. (83) 

We note, that the so chosen (3 does not define in general an integrable 1-dimensional Pfaff system, so, 
accordingly: da A a = 0, 6.(3 A (3 ^ 0. 
At these conditions our equations 

5*FAF = aAF, 5F A*F = (3 A*F, 5 * F A *F - 5F A F = 0, 

daAaA/3 = 0, d(3 A a A (3 = 
take the form: 5 * F A *F — 5F A F = is reduced to 

— vu y + uvy = 0, —uv x + V u x = 0, 

the Frobenius equations da A a A (3 = 0, d(3 A a A (3 = reduce to 

{-b x B + B x b) .A = 0, 

5*FAF = aAFis reduced to 

u (u^ — v z ) = 0, v (u^ — v z ) = 0, uu x — vv x = Au, 
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Bu. 



finally, SF A *F = (3 A *F reduces to 

v (v£ — u z ) = —bv, u (v£ — u z ) = —bu, uu y — vv y 

In this way we obtain 7 equations for 5 unknown functions u, v, A, B, b. 

The equations — vu y + uv y = 0, — uv x + vu x = have the following solution: 

u{x,V,z,0 = f(x,y)U(z,£), v{x,y,z,C) = f(x,y)V(z,£). 

That's why 

AU = f x {U 2 -V 2 ), BU = f y {U 2 -V 2 ), f{Vz-U g ) = -b, Us-V z = 0. 

It follows that b should be of the kind b(x,y, z,£) = f (x , y)b° '{z , £) , so the equation B x b — Bb x = 
takes the form 

ffxy = fxfy 

The general solution of this last equation is f(x,y) = g(x)h(y). The equation gh(V^ — U z ) = —b 
reduces to 

Vz-U x = -b°. 

The relations obtained show how to build a solution of this class. Namely, first, we choose the 
function V(z,£), then we determine the function U(z,£) by 



U(z,0 = J V z d£ + l(z), 



where l(z) is an arbitrary function, which may be assumed equal to 0. After that we define b° = U z — V^. 
The functions g(x) and h(y) are arbitrary, and for A and B we find 

, U 2 - V 2 U 2 -V 2 
A{x, y, z,£)=g {x)h(y) , B(x, y, z, f ) = g(x)h (y) . 

In this way we obtain a family of solutions, which is parametrized by one function V of the two 
variables (z,£) and two functions g(x), h(y), each depending on one variable. 

In order to find corresponding conserved quantities we sum up the nonzero right-hand sides of the 
equations and obtain *(a A F + (3 A *F). We represent this expression in a divergence form as follows: 



*(a A F + A *F) = Audx - Budy - budz - bvd£ 
= \{U 2 - V 2 ) [{gh) 2 ] x dx - \{U 2 - V 2 ) [{gh)\dy- 



-(9hf(J Ub°dz^j dz-(gh) 2 ^J Vb d^di==-^P^dx\ 

where the interaction energy-momentum tensor is defined by the matrix 

-\{gh) 2 Z 





o Ughfz 



o 






^2 T/2 ; c 



(gh) 2 JUb°dz 

(gh) 2 f Vb°d^ 



and the notation Z = U — V is used. For the components of the full energy tensor = + 
we obtain 

Ub°dz-\{U 2 + V 2 ) 



rp3 

1 3 



(ghY 



/■ 



rp4 

J 3 



-Tf 



2 

(gh) 2 UV, 



J Vb°d£ + ^(U 2 + V 2 



and all other components are zero. 
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11.2 Examples 



In this subsection we consider some of the well known and well studied (l+l)-dimensional soliton 
equations as generating procedures for choosing the function V(z,£), and only the 1-soliton solutions 
will be explicitly elaborated. Of course, there is no anything standing in our way to consider other 
(e.g. multisoliton) solutions. We do not give the corresponding formulas just for the sake of simplicity. 

We turn to the soliton equations mainly because of two reasons. First, most of the solutions have 
a clear physical sense in a definite part of physics and, according to our opinion, they are sufficiently 
attractive for models of real physical objects with internal structure. Second, the soliton solutions 
describe free and interacting objects with no dissipation of energy and momentum, which corresponds 
to our understanding of the Frobenius integrability equations as criteria for availability of integral 
properties of the medium considered. 



1. Nonlinear (1+1) Klein-Gordon Equation. In this example we define our functions U and V through 



the derivatives of the function k(z, £) in the following way: U = k z , V = k^. 
U i -V z 
kzz ~ kf 



Then the equation 

k z £ — k^ z = is satisfied automatically, and the equation U z — V£ = b° takes the form 
= b°. Since b° is unknown, we may assume b° = b°(k), which reduces the whole problem to 
solving the general nonlinear Klein-Gordon equation when b° depends nonlinearly on k. Since in this 
case V = k^ we have 



Vb°{k)d£ = J k^b°{k)di 
For the full energy density we get 



d_ 



kl + kj + 2 



b°{k)dk 



d£ 



J b°(k)dk. 



b°(k)dk 



Choosing b°(k) = m 2 sin(k), m = const, we get the well known and widely used in physics Sine-Gordon 
equation, and accordingly, we can use all solutions of this (l+l)-dimensional nonlinear equation to 
generate (3+l)-dimensional solutions of our equations following the above described procedure. When 
we consider the (3+1) extension of the soliton solutions of this equation, the functions g(x) and h(y) 
should be localized too. The determination of the all 5 functions in our approach is straightforward, 
so we obtain a (3+l)-dimensional version of the soliton solution chosen. As it is seen from the above 
given formulas, the integral energy of the solution differs from the energy of the corresponding (1+1)- 
dimensional solution just by the (x, y)-localizing factor [g(x)h(y)] 2 . 
For the 1-soliton solution (kink) we have: 



k(z, £) = Aarctg <^exp 
±2m 



±-(*--0 

7 c 



7 



U(z,0 = k z = 

A = g'(x)h(y) 



m 
' 7 



X) 



±2777-7 



ch 



, B = g(x)h'(y) 



w 2 



±2mw 



C1 ch [±2 ( z - H£) 
±27777 



ch 



b° = U z - V e 



and for the 2-form F we get 



-27T7 2 


sh 




ch 







(gh) 2 4m 2 



l 2 ch 2 ±^(z-H C ) 



±2mq(x)h(y) , , w ±2mq(x)h(y) 
F = r yy ' yy ' n dy Adz + yy ' yyj dy A d£. 



jch ±f(z-fC) 



7C/7 
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In its own frame of reference this soliton looks like 

= ±2mg(x)h{y) 
ch(±mz) 



dy A dz. 



From this last expression and from symmetry considerations, i.e. at homogeneous and isotropic 
medium, we come to the most natural (but not necessary) choice of the functions g(x) and h(y): 

1 i 
9( x ) = -T7 — \) Kv) = 



ch(mx) ' ch(my) 

2.Korteweg-de Vries equation. This nonlinear equation has the following general form: 

ft + aiff z + a 2 fzzz = 0, 
where a\ and a>2 are 2 constants. The well known 1-soliton solution is 



t , a a T „ / 3a 2 ca a ai 

fM = c* [i-zzy i = 2 fc w = — 

where a Q is a constant. We choose V(z,£) = f(z,£) and get 

v>ch*[i-%ty Uz ^ ^Lw c)eh*[i-%ty 

Tt = {ghf a ^ l + L ) 



2w 2 Lch 4 [, _ ^] ■ 

3. Nonlinear Schrddinger equation. In this case we have an equation for a complex- valued function, 
i.e. for two real valued functions. The equation reads 

ift + fzz + \f\ 2 f = 0, 

and its 1-soliton solution, having oscillatory character, is 

2 exp [-i (2az + 4(a 2 - /3 2 )£ - 9)] 
" 2(3 ch(2(3z + 8afc-5) ' 

where a, (3, 5 and are constants. The natural substitution f(z, £) = ^/p.exp(iip) brings this equation 
to the following two equations 

Pt + {2p<p z ) z = 0, 4p+ 2pPzz 2 ~ Pz =4(yg + ^). 



For the 1-soliton solution we get 

4/3 4 



r ch 2 (20z + 8af3£ - 5) 
We put U = p, V = —2pip z and obtain 



, <p = - [2az + 4(a 2 - /3 2 )£ - 6] . 



rr 4/3 4 16a/3 4 



ch 2 (2(3z + 8a/?£ - S) ' ^* c/i 2 (2/?z + 8a/?£ - ,5) ' 

„ , 4/3 4 (l-16a 2 ) „ , .,„ , 4/? 4 (l-16a 2 ) 

,4=g'(x)%) , 2 /V> — -/ jr , B = g{x)h'(y)- 



ch 2 (2[3z + 8a(3Z- 5)' y ' y vw/ c/i 2 (2/3 z + 8a(3£ - 5) ' 
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bQ _ 16/^(16a 2 /3 2 - l)sh{2f3z + 8a(3£ - 5) 



T% = (gh) 



ch 3 (2pz + 8al3£-5) 
16/3 8 



c/i 4 (2/3z + 8a/?£-5)' 

We note that the solution of our equations obtained has no the oscillatory character of the original 
Schrodinger 1-soliton solution, it is a (3+l)-localized running wave and moves as a whole with the 
velocity 4ca, so in the relativistic frame we must require 4ca < c. 

A.Boomerons. The system of differential equations, having soliton solutions, known as boomerons, 
is defined by the following functions: K : TZ 2 — > 7£ 3 , H : TZ 2 — > 7£, and besides, two constant 
3-dimensional vectors r and s, where s is a unit vector: |s| = 1. The equations have the form 

- s.K z = 0, K 2 £ = H zz s + r x K z - 2K Z x (K x s). 

Now we have to define our functions U(z,£), V(z,£) and b°(z,£). The defining relations are: 

U = H z = \s\ 2 H z , V = s.K z = (s.K) 2 , b° = -s. [r x K z - 2K Z x (K x s)] . 

Under these definitions our equations — V z = 0, V^ — U z = —b° look as follows: 

[H z - s.K z ] z = 0, s. [Kx - H zz s - r x K z + 2K Z x (K x s)] = 0. 

It is clear that every solution of the "boomeron" system determines a solution of our system of 
equations according to the above given rules and with the multiplicative factors g(x) and h(y). The 
two our functions A{z,£) and B(z,£) are easily then computed. 

Following this procedure we can generate a spatial soliton solution to our system of equations by 
means of every solution, in particular every soliton solution, to any (l+l)-soliton equation, as well as to 
compute the corresponding conserved quantities. It seems senseless to give here these easily obtainable 
results. The richness of this comparatively simple family of solutions, as well as the availability of 
corresponding correctly defined integral conserved quantities, are obvious and should not be neglected 
as a generating procedure for obtaining (3+1) one-soliton, as well as, (3+1) many-soliton solutions. 
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12 Retrospect and Outlook 



In trying to understand our observational knowledge of the real world we must be able to separate the 
important structural and behavioral properties of the real objects from those, the changes of which 
during time-evolution do not lead to annihilation of the objects under consideration. One of the basic 
in our view lessons that we more or less have been taught is that the physical objects are spatially finite 
entities, and that for their detection and further study some energy-momentum exchange is necessarily 
required. So, every physical object necessarily carries energy-momentum and every interaction between 
two physical objects has such an energy-momentum exchange aspect. The second lesson concerning 
any interaction is that, beyond its universality, energy-momentum is a conserved quantity, so NO loss 
of it is allowed: it may only pass from one object to another. This means that every annihilation 
process causes creation process (es), and the full energy- momentum that has been carried by the 
annihilated objects, is carried away by the created ones. Energy-momentum always needs carriers, 
as well as every physical object always carries energy-momentum. Hence, the energy-momentum 
exchange abilities of any physical object realize its protection against external influence on one side, 
and reveal its intrinsic nature, on the other side. Therefore, our knowledge about the entire complex 
of properties of a physical object relies on getting information about its abilities in this respect and 
finding corresponding quantities describing quantitatively these abilities. 

The spatially finite nature of a physical object implies spatial structure and finite quantity of 
energy-momentum needed for its creation, so NO structureless and infinite objects may exist. The 
approximations for "point object" and "infinite field", although useful in some respects, seem the- 
oretically inadequate and should not be considered as basic ones. More reliable appears to be the 
approximation "finite continuous object", which we tried to follow throughout our exposition. This 
last approximation suggests that nonlinear partial differential equations should be the basic tool for 
building mathematical models of local nature of such objects. The natural physical sense of these 
equations is naturally supposed to be local energy-momentum exchange. 

Another useful observation is that physical objects are many-aspect entities, they have complicated 
structure and their very existence is connected with internal energy- momentum exchange/redistribution 
among the various structural components. So, the mathematical model objects should be many- 
component ones, and with appropriate mathematical structure. Of basic help in finding appropriate 
mathematical objects is having knowledge of the internal symmetry properties of the physical object 
under consideration. This "step by step" process of getting and accumulating important information 
about the physical properties of natural objects reflects in the "step by step" process of refining the 
corresponding mathematical models. 

The greatest discovery at the very beginning of the last century was that the notion of electro- 
magnetic field as suggested by Maxwell equations is inadequate: the time dependent electromagnetic 
field is not an infinite smooth perturbation of the aether, on the contrary, it consists of many individ- 
ual time-stable objects, called later photons, which are created/destroyed mainly during intra-atomic 
energy-transition processes. Photons are finite objects, they carry energy-momentum and after they 
have been radiated outside their atom-creator, they propagate as a whole translationally by the speed 
of light. Moreover, their propagation is not just translational, it includes rotational component, which 
is of intrinsic and periodical nature. The corresponding intrinsic action for one period T is h = ET, 
where E is the full energy of the photon, and all photons carry the same intrinsic action h. During 
the entire 20th century physicists have tried to understand the dynamical structure/nature of photons 
from various points of view, and this process is still going on today. 

Extended Electrodynamics, as presented in the preceding sections, is an attempt in this direc- 
tion. The basic starting observation for approaching the problem is that the energy-momentum local 
quantities and relations of Maxwell theory do agree with the experiment, but the free field equations 
5F = 0, 5 * F = give non-realistic free time-dependent solutions: they are either strongly time- 
unstable, or infinite. Hence, these solutions can not be used as mathematical models of photons, since 
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the latter are time-stable and finite objects. The basic idea of writing down new nonlinear equa- 
tions was to pass to local energy-momentum relations, describing how the internal energy- momentum 
exchanges are carried out. The mathematical structure of Maxwell equations and local conserva- 
tion relations allowed the extension procedure to be used, so the new nonlinear equations contain all 
Maxwell solutions as exact solutions. This feature is important from the point of view of applications. 

In order to come to the new equations we made use of the dual symmetry of the linear solutions, 
which suggested to consider the field as having two vector components (F,*F), and introduced the 
concept of G-covariant Lorentz force as a natural generalization of the known Lorentz force. The G- 
covariant Lorentz force has three vector components, and each component has a clear physical sense 
of corresponding local energy-momentum change, as it should be. The two components F^ V 5F V and 
(*F)fj, u (5 * F) v determine how much energy- momentum the field is potentially able to give out to 
some other physical object through any of the two components F and *F, and the third component 
F^ u {5 * F) u + (*F) tlu SF u characterizes the internal energy-momentum exchange between the two 
components F and *F. We concluded that in case of energy-momentum isolation, i.e. in vacuum, 
these three components must be zero, so, the dynamics obtained is of intrinsic nature. 

In studying the vacuum nonlinear solutions, i.e. those with SF ^ and 5 * F / 0, we found 
their basic property: every nonlinear solution has zero-invariants: F^F^ = [*F) llv F^ u = 0. Then 
the eigen properties of F, *F and of the energy-momentum tensor Q u ^ determine unique isotropic 
direction ( along which the solution propagates translationally as a whole, which fits well with the 
photons' way of propagation . The obtained simple form of F = A A ( and *F = A* A ( allowed 
a complete analysis of the nonlinear solutions to be made. The whole set of nonlinear solutions 
consists of nonoverlapping subsets, and each subset is characterized by the corresponding isotropic 
eigen direction. Every solution of a given subclass is uniquely determined by two functions: the 
amplitude function 4> which is arbitrary with respect to the spatial variables and is a running wave 
along (; the phase function <p, \<p\ < 1, where (p depends arbitrarily on all space-time variables (we 
made use of the function tp = arccos{(p)). The components of F have the form 4>cosip and (psinip, so, 
finite solutions with photon-like behavior are allowed. 

A basic characteristic of the nonlinear solutions with \SF\ / turned out to be the scale factor C 
defined by the relation £ = \A\/\8F\ and having physical dimension of [length]. Hence, every nonlinear 
solution with \b~F\ ^ defines its own scale. As it was shown further, the case C = const is allowed, 
so corresponding intrinsically defined time-period T = C/c and frequency v = T~ l can be introduced, 
and the corresponding finite solutions with integral energy E admit the characteristic intrinsic action 
given by ET. 

The natural question " do there exist nonlinear finite solutions with rotational component of prop- 
agation" was answered positively. Seven equivalent conditions of quite different nature, determining 
when this is possible, were found to exist. It is remarkable that the condition C ^ is one of them. 
So, for finite solutions with C = const, we have a natural and intrinsically defined measure of this 
rotational component of propagation, namely, the elementary action ET, which recalls the correspond- 
ing invariant characteristic of photons h, the Planck constant. Anyway, these solutions deserve to be 
called " photon- like" . It is also remarkable that completely integrable and nonintegrable 2-dimensional 
Pfaff systems can be associated with a nonlinear solution, and a natural integral characteristic of the 
nonintegrable Pfaff systems is also the quantity ET. This point shows that both, integrable and non- 
integrable features of a solution are important, the nonintegrability characteristics give information 
about the contacts between the system and the physical environment. 

It was very interesting to find that the energy-momentum tensor Q of a nonlinear solution defines 
a boundary operator in the tangent and cotangent bundles of the Minkowski space-time. The corre- 
sponding homology /cohomology spaces are 2-dimensional, the classes represented by the electric and 
magnetic components of the field F form a basis of the corresponding homology/ cohomology space. 
Moreover, since for each nonlinear solution the corresponding image spaces of Q are 1-dimensional, Q 
is extended to boundary operator in the whole algebra of exterior forms and antisymmetric tensors. 
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The Q-image of any differential 2-form is collinear to the nonlinear solution that generates Q, and the 
Q-image Q(a) of any p-form a is isotropic: [(3(a)] 2 = 0. 

The natural representation of G in the space of 2- forms leaves the scale factor C invariant. As 
a consequence we obtained that C is Lorentz- invariant only if it is a constant, and in this case the 
3-form SF A F is closed, which generates a conserved quantity through Stokes theorem, and this 
conserved quantity is proportional to the elementary action ET. If the group parameters depend on 
the space-time points, then the commutative group structure of G generates a group structure inside 
the subset of solutions with the same £. In such a case a "vacuum state" F Q can be defined such, that 
every solution of the subclass is defined by an action upon F Q of a point-dependent group element 
with determinant (a 2 + b 2 ) having running wave character along £. 

The connection-curvature interpretations showed that rotational component of propagation may 
have only those nonlinear solutions which generate non-zero curvature in the three cases of general, 
principal amd linear connection, so, roughly speaking, nonzero curvature means nonzero spin. 

We gave five ways to compute the integral spin ET of a nonlinear solution. It worths noting that 
an appropriate representative of the Godbillon-Vey class determined by the completely integrable 1- 
dimensional Pfaff system £, or f(, can be made equal to 5F A F, and so it can be used to compute 
the elementary action ET. 

Finally, in presence of external fields, we showed that our general system of nonlinear equations 
together with the additional Pfaff equations is consistent, and we found a large family of solutions. 
This family is parametrized by one function of one-space (say z) and one-time independent variables: 
V(z,£,), and two other functions g(x) and h{y). So, choosing V(z, £) to be any (one, or many) 
soliton solution of any soliton equation, and g(x)&nd h{y) to be finite, we obtain its (3+1) image as a 
finite/concentrated solution of our nonlinear equations with well defined energy-momentum quantities. 
We illustrated our approach with examples from the 1-soliton solutions of the well known Sine-Gordon, 
KdV and NLS equations. 

We would like to note also that the nonlinear vacuum equations obtained follow the idea that 
was formulated in Sec. 1.1: the admissible changes (in our case 5 * F and 5F, or correspondingly, 
dF and d * F) are projected upon the field components (F,*F), and these projections are assumed 
to be zero, or intrinsically connected. This leads directly to nonlinear equations with corresponding 
physical sense of energy-momentum balance relations. From mathematical point of view this suggests 
a natural generalization of the geometrical concept of parallelism, and this idea is developed and 
illustrated with many examples from differential geometry and theoretical physics in the Appendix. 
This generalization, in turn, suggests natural ways to nonlinearization of important physical linear 
and nonlinear equations. 

The development of theoretical physics during the last century, and especially during the past 
30 years, shows growing interest to nonlinearization of the widely used linear equations. At the 
beginning of the century Einstein declared [32] that the fundamental equation of Optics = must 
be replaced by nonlinear equation(s), and he had worked in trying to find appropriate one(s). General 
Relativity made a decisive step along the road of nonlinearization. Mie's [44] and Born-Infeld's [9] 
nonlinearizations of CED are also steps in this direction. One of the greatest achievements was the 
Yang-Mills approach, which dominates nowadays the various models of field theory (the "standard 
model" ) . Unfortunately, spatially finite solutions of the vacuum Yang- Mills equations with soliton 
behavior were not found. 

The later studies with appropriate interaction term in the Yang-Mills lagrangian brought the 
monopole solutions [33]. A somewhat curious example is presented by the instanton solutions [34] 
because of their mathematical significance on one side, and because of the interpretational difficulties 
coming from the necessary zero value of their energy-momentum tensor, on the other side. Some 
development should be noted of the theory of the so called " nontopological solitons", or Q-balls, 
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[35,36,37,38], dyons [39], vortices [40], sphalerons [41], skyrmeons [42], knots [43], which in most cases 
give spherically symmetric static solutions of equations obtained from correspondingly nonlinearized 
lagrangians. We have to mention that all these examples follow the rule: "add nonlinear interaction 
term to the " free" lagrangian giving linear equations and see what happens" . We do not share the view 
that this is the right way to pass to (3+l)-dimensional nonlinear equations. We explained our view 
in the preceding sections and showed that our approach works in the important case of the vacuum 
photon- like finite solutions and in presence of external fields. All free and not-free time-stable objects 
in Nature are finite, so, the widely used "free" terms in any lagrangian could play some misleading 
role, especially those "free" terms that generate, or imply, D'Alembert-like wave equations. 

Finally, we just mention (without any comments) the (super)string and brane approaches, which 
also try to incorporate in theoretical physics spatially finite solutions. 
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A Appendix: Extended Parallelism. Applications in Physics 



Abstract 

This appendix considers a generalization of the existing concept of parallel (with respect to a 
given connection) geometric objects and its possible usage as a suggesting rule in searching for 
adequate field equations and local conservation laws in theoretical physics. The generalization tries 
to represent mathematically the two-sided (or dual) nature of the physical objects, the change 
and the conservation. The physical objects are presented mathematically by sections IP of vector 
bundles, the admissible changes are described as a result of the action of appropriate differential 
operators D on these sections, and the conservation proprieties are accounted for by the requirement 
that suitable projections of on *ff and on other appropriate sections must be zero. It is shown 
that the most important equations of theoretical physics obey this rule. Extended forms of Maxwell 
and Yang-Mills equations are also considered. 

A.l The general Rule 

We are going to consider here a more general view on the geometrical concept of parallelism, more 
or less already implicitly used in some physical theories. The concept of parallelism appropriately 
unifies two features: change and conservation, through suitable differential operators and suitable 
"projections". In some cases this concept may be applied to produce dynamical equations, and in other 
cases the corresponding physical interpretation is conservation (balance) equations (mainly energy- 
momentum balance). The examples presented further show how it has been used so far, and how it 
could be used as a generating tool for producing new useful relations and for finding natural ways to 
extension of the existing equations in physics. An important feature, that deserves to be noted even 
at this moment, is that the corresponding equations may become nonlinear in a natural way, so we 
might be fortunately surprised by appearing of spatially finite (or soliton-like) solutions. 

We begin with the algebraic structure to be used further in the bundle picture. The basic concepts 
used are the tensor product <g) of two linear spaces (we shall use the same term linear space for a vector 
space over a field, and for a module over a ring, and from the context it will be clear which case is 
considered) and bilinear maps. The bilinear maps used are as a rule assumed to be nondegenerate [1] . 
Also, all linear spaces, manifolds and bundles are assumed to be finite dimensional, and the dimensions 
will be pointed out only if needed. 

Let (Ui,Vi), (U 2 , V 2 ) and (U3, V3) be three couples of linear spaces. Let $ : U\ x U 2 — > U3 and 
tp : V\ x V 2 — > V3 be two bilinear maps. We consider the decomposable elements 

(ui <g> v±) G Ui <g> Vi and (u 2 <8> v 2 ) G U 2 <8> V 2 . 
Now we make use of the given bilinear maps as follows: 

($, <p)(ui <g> vi,u 2 <8> v 2 ) = u 2 ) <8> ip(yi,v 2 ), 

and on nondecomposable elements we extend this map by linearity: 

($, ¥>)((«!)* ® («l)i, {U 2 ) j ® {V 2 )j) = $((«!)*, (U 2 ) j ) ® ¥>((ui)i, {V 2 )j), 

where summation over the repeated indecies is assumed. The obtained elements are in U3 ® V3. Note 
that if {ej} is a basis of V then for every element u®v£U<£>Vwe have u® v = u® v l e.{ = v l u ® ei = 
K % <g> a, where K l = v l u G U, so further we may use such a representation. 

We give now the corresponding bundle picture. Let M be a smooth real manifold. We assume 
that the following vector bundles over M are constructed: £i,r/j, with standard fibers Ui,Vi and sets 
of sections Secfa), Sec(rji),i = 1,2,3. 



74 



Assume the two bundle maps are given: ($,zc/m) : x — > ?3 an d (v^^m) : Vi x V2 - > %• Let 
{e^} and {Tj} be local bases of Sec{r]i) and Sec{r]2) respectively. Then if A* (g) and -B- 7 <8> Tj are 
sections of £1 ® 771 and £2 ® ^2 respectively, we can form an element of Sec{^ ® 773) by the rule : 

(^^(^(gei.B^Tj) = *(A i ,B J ')® ¥ >(e i ,r j ). (84) 

Note that this rule assumes that the summation commutes with the bilinear maps $ and ip. 

Let now £ be a new vector bundle on M and u> G Sec{^). Assume that we have a differential 
operator D : Sec(£) — > Sec{^2 ®ffi) s ° we can form the section Dm 6 >Sec(£2 <8) 772) an d represent it as 
(Dw) 3 ® Tj G Sec(^2 <8> ^2)- Applying now the two bilinear maps we obtain 

${A\ (Dw) j ) ® ^(ei,rj) G <Sec(£ 3 ® 773). 

We give now the following 

Definition: The section w G Sec(£), called further the section of interest, will be called ($>,tp;D)- 
parallel/ conservative with respect to A 1 ® e$ G Sec(£;i ® 771) if 

(S.^DXA*® <*,«;) = ($,^(,4'® ®7}) = $(A\ (£w;) J ') ® r,) = 0. (85) 

The term parallel seems more appropriate when, making use of (85), the relation obtained is interpreted 
as equation of motion (field equation) , and the term conservative is more appropriate when the relation 
obtained is interpreted as local conservation law. 

The above relation (85) we call the GENERAL RULE (GR). The bilinear maps ip "project" 
the "changes" Dw of the section w on the section A 1 ® e% in the bases chosen. If the relation (85) 
holds we understand it in the following sense: those characteristics of w which do feel the operator D 
may be NOT constant with respect to D, but the "projections" of the possible nonzero changes Dw 
on an appropriate section A 1 ® e, do vanish. In particular, if A 1 ® a = w, then we may say that w is 
(<3? , ip; D)- autoparallel. In this sense the above definition (85) goes along with our concept of physical 
object as given in Sec.l of this paper. 

As an example of a differential operator we note the exterior derivative d : h p (M) — > A P+1 (M). 
In the case of the physically important example of Lie algebra 0-valued differential forms, with "$ = 
exterior product" and "<p = Lie bracket [,]", £1 = A P (M), £ 2 = A P+1 (M), f = A P (M) ® g, 771 = 772 = 
M x g, the GR (85) looks as follows: 

(A, [, ]; d x id fl )(a* ® ® £7) = (A, [, ])(a* ® d/3 J ' ® £,■) = a* A d/3 J ' ® Ej\ = 0, 

where {.Ej} is a basis of g, and a summation over the repeated indexes is understood. Further we are 
going to consider particular cases of the (GR) (85) with explicitly defined operators D and additional 
operators whenever they participate in the definition of the corresponding section of interest. 

A. 2 The General Rule in Action 

A. 2.1 Examples from Geometry and Mechanics 

We begin with some simple examples from differential geometry and classical mechanics. 

1. Integral invariance relations 

These relations have been introduced and studied from the point of view of applications in me- 
chanics by Lichnerowicz [2]. 

We specify the bundles over the real finite dimensional manifold M: 

£1 = TM; £ 2 = T*(M); 771 = 772 = £3 = m = M x R, denote Sec(M xR) = C°°(M) 
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3>= substitution operator, denoted by i(X),X G Sec(TM); 
y=point-wise product of functions. 

We denote by 1 the function f(x) = l,x G M. Consider the sections 
J0l€ Sec(TM (g) (M x R)); a®l 6 Sec(T*M (8 (M x R)). Then the GR leads to 

<p)(X <8 1, a <8 1) = i(X)a <8 1 = i(X)a = 0. (86) 

We introduce now the differential operator d: if a is an exact 1-form, a = df, so that £ = M x R, 
the GR gives 

($, <p; d)(X ® 1, /) = ® 1, d/ ® 1) = i(X)d/ (1 1) = i(X)d/ = = 0. 

1. e. the derivative of / along the vector field X is equal to zero. So, we obtain the well known relation, 
defining the first integrals / of the dynamical system determined by the vector field X. In this sense 
/ may be called ($, ip, d)- conservative with respect to X, where <3? and <p are defined above. 

2. Absolute and relative integral invariants 

These quantities have been introduced and studied in mechanics by Cartan [3]. By definition, a 
p-form a is called an absolute integral invariant of the vector field X if i(X)a = and i(X)da = 0. 
And a is called a relative integral invariant of the field X if i(X)da = 0. So, in our terminology (the 
same bundle picture as above), we can call the relative integral invariants of X (<E>, ip; d)- conservative 
with respect to X, and the absolute integral invariants of X are additionally (<1>, ip)- conservative with 
respect to X, with (<&, ip) as defined above. A special case is when p = n, and oj G A n (M) is a volume 
form on M. 

3. Symplectic geometry and mechanics 

Symplectic manifolds are even dimensional and have a distinguished nondegenerate closed 2-form 
lo, duj = 0. This structure may be defined in terms of the GR in the following way. Choose £i = 
tji = t)2 = MxR,^i|2 = A 3 (T*M) ®R, and d as a differential operator. Consider now the section 
1 G Sec(M x R) and uj G Sec(A 2 (T*M), with u - nondegenerate. The map <& is the product f.uj and 
the map ip is the product of functions. So, we have 

($, (p; d)(l ®l,w) = $(1, da;) (g) 99(1, 1) = do; = 0. 

Hence, the relation do; = is equivalent to the requirement a; to be (<3?, 99; d)-parallel with respect to 
the section 1(8)1. 

The hamiltonian vector fields X are defined by the condition L^w = dz(X)o; = 0. If 5> = </? is the 
point-wise product of functions we have 

</?; d)(l (8 1, i(X)w) = ($, <p)(l (8 1, di(X)(o;) 1) = L x o; 1 = = 0. 

In terms of the GR we can say that X is hamiltonian if i(X)u is (3>, d)-parallel with respect to 
1 ® 1. 

The induced Poisson structure {/, g}, is given in terms of the GR by setting <I> = a; -1 , where 
a; -1 . a; = idxM, v?=point-wise product of functions, and 1 G Sec(M x R). We get 

($,p)(d/<g>l,d 5 <g>l) =uj- 1 (df,dg)®l. 

A closed 1-form a, da = 0, is a first integral of the hamiltonian system Z, di{Z)u) = 0, if i(Z)a = 
0. In terms of the GR we can say that the first integrals a are (i, (^-parallel with respect to Z: 
(i,ip)(Z (8 l,a <8 1) = i(Z)a (8 1 = 0. From Lzoj = it follows Lz^ 1 = 0. The Poisson bracket 
(a, (3) of two first integrals a and (3 is equal to [— d(w~ 1 (a, /?))] [6]. The well known property that 
the Poisson bracket of two first integrals of Z is again a first integral of Z may be formulated as: the 
function u;" 1 (a,/3) is (i, <p; d)-parallel with respect to Z, 

(i,ip;d)(Z (gsl,^- 1 ^,^)) = i(Z)duj- 1 (a,P)®l = 0. 
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A. 2. 2 Probenius integrability theorems and linear connections 

l.Frobenius integrability theorems [4] 

Let A = (X±, . . . , X r ) be a differential system on M, i.e. the vector fields Xi, i = 1, . . . , r define 
a locally stable submodule of Sec(TM) and at every point p £ M the subspace A£ C T p (M) has 
dimension r. Then A r is called involutive (or integrable) if G A r ,i,j = 1, . . . , r. Denote by 

A£- r C T P (M) a complimentary subspace: A r p © A r p l ~ r = T P (M), and let vr : T P (M) -► A"~ r be the 
corresponding projection. So, the corresponding Frobenius integrability condition means ir([Xi, Xj]) = 
= 1, • • • ,r. 

In terms of the GR we set D{X,j) = it o Lx tl 3>=" product of functions and vector fields", and ip 
again the product of functions. The integrability condition now is 

= <p)(l © 1, ir([Xi, Xj] © 1)) = l.ir{[Xi, Xj}) ©1 = 0, i,j = l,...,r. 

In the dual formulation we have the Pfaff system A*_ r , generated by the linearly independent 
1-forms (ai, . . . , a n _ r ), such that a m (Xi) = 0, i = 1, . . . r; m = 1, . . . n — r. Then A*_ r is integrable if 
da A a.\ A • • • A a n _ r = 0, a G In terms of GR we set ip the same as above, $ = A and d as 

differential operator and obtain 

(<E>, <p; d)(ai A ■ ■ ■ A a n _ r © 1, a) = (A, </?) (ai A ■ ■ ■ A a n _ r © 1, da © 1) 

= da A «i A • • • A a n _ r © y>(l, 1) = da A ol\ A ■ ■ ■ A a n _ r = 0, 

1. e. a is (<&, d)-parallel with respect to ai A • • • A a ra _ r © 1. 

2. Linear connections 

The concept of a linear connection in a vector bundle has proved to be of great importance in 
geometry and physics. In fact, it allows to differentiate sections of vector bundles along vector fields, 
which is a basic operation in differential geometry, and in theoretical physics the physical fields are 
represented mainly by sections of vector bundles. We recall now how one comes to it. 

Let / : W l —> R be a differentiable function. Then we can find its differential df. The map / — > df 
is R-linear: d(n.f) = K.df, k G R, and it has the derivative property d(f.g) = fdg + gdf. These two 
properties are characteristic ones, and they are carried to the bundle situation as follows. 

Let £ be a vector bundle over M. We always have the trivial bundle („ = M x R. Consider now 
/ G C°°(M) as a section of £ G . We note that Sec(£ G ) = C°°(M) is a module over itself, so we can 
form df with the above two characteristic properties. The new object df lives in the space A 1 (M) 
of 1-forms on M, so it defines a linear map df : Sec(TM) — > Sec(£ ),df(X) = X(f). Hence, we 
have a map V from Sec(^ ) to the 1-forms with values in Sec(£ ), and this map has the above two 
characteristic properties. We say that V defines a linear connection in the vector bundle £ . 

In the general case the sections 5ec(£) of the vector bundle £ form a module over C°°(M). So, a 
linear connection V in £ is a R-linear map V : Sec(£) — > A 1 (M, £). In other words, V sends a section 
a G 5'ec(£) to a 1-form Va valued in Sec(^) in such a way, that 

V{ka) = kV{a), V(/ a) = df © a + f V(a), (87) 

where fc G R and / G C°°(M). If X G Sec(TM) then we have the composition i(J)oV,so that 

i(X)oV(fa) = X(f)a + fV x (a), 

where Vx(ff) G 5ec(£). 
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In terms of the GR we put ^ = TM = f and £ 2 = A 1 (M) ® £, and 771 = r? 2 = Co- Also, 
$(A", Vcr) = Vx<7 and ip(f,g) = f.g. Hence, we obtain 

($, ¥>; V)(X 1, a) = ($, </>)(X ® 1, (Vff) 1) = V x a ® 1 = V x ct, (88) 

and the section a is called \ '-parallel with respect to X if V^cr = 0. 

3. Covariant exterior derivative 

The space of ^-valued p-forms A P (M, £) on M is isomorphic to A p (M)(g)Sec(£). So, if (<7i, . . . , cr r ) is 
a local basis of Sec(£), every ^ G A P (M, £) is represented by ip l ® Cj, i = 1, . . . , r, where ■0* £ A P (M). 
Clearly the space A(M,£) = £™ =0 A P (M,£), where A°(M,£) = Sec(0, is a A(M) = EJ =0 AP(M)- 
module: q.^ = q A ^ = (a A tp 1 ) <8> Cj- 

A linear connection V in £ generates covariant exterior derivative D : A P (M, £) — ► A P+1 (M, £) in 
A(M, £) according to the rule 

D* = D(V>* ® <7i) = <7i + (-l)V A V(<Ti) 

= (d^ + (-l)V A r^-d^) ® ^ = (Dtf )' ® <7i. 

We may call now a ^-valued p-form ^ V '-parallel if = 0, and (X, V)-parallel if i(X)D\l/ = 0. This 
definition extends in a natural way to q- vectors with q < p. Actually, the substitution operator i(X) 
extends to (decomposable) (/-vectors X\ A X2 A • • • A X q as follows: 

A X 2 A • • • A = o o • • • o 

and extends to nondecomposable (7- vectors by linearity. Hence, if is a section of A q (TM) we may 
call * (6, V)-parallel if i(0)D* = 0. 

Denote now by the vector bundle of (linear) homomorphisms (n, id) : £ — > £, and let n G 
Sec(Lt). Let % G Sec(A q (TM)®L £ ) be represented as 0®n. The map $ will act as: $(0, *) = z(0)*, 
and the map 93 will act as: <p(H, <7j) = n(<jj). So, if V(crfc) = T 3 ^ k dx^ ® 0j, we may call \& (V)-parallel 
with respect to % = (g) n if 

($, 99; D)(© <g> n, * = f ® a-i) = ($, v?)(© n, (d*)* ® oi) = i{e)(D^y ® n(^) = 0. (89) 

If we have isomorphisms ® P TM ^ ® p T*M,p = 1,2,..., defined in some natural way (e.g. through 
a metric tensor field), then to any p-form a corresponds unique p- vector a. In this case we may talk 
about autopaparallel objects with respect a (point-wise) bilinear map <p : (£ x £) — ► 77, where 77 

is also a vector bundle over M. So, ^ = a k Ofc G A P (M, £) may be called (i, ip; V)-autoparallel with 
respect to the isomorphism if 

(i,y?; V)(d fc 0d t ,« m ®ffm) 

= i(a fc )da m <p(a k ,a m ) + (-l) p i(a k )(^ A I^ds") <p(a k ,a m ) (90) 
= [i(a fc )da m + (-l) p i(a fc )(^ A T™^)] ® ^(a k ,a m ) = 0. 

Although the above examples do not, of course, give a complete list of the possible applications of the 
GR (85), they will serve as a good basis for the physical applications we are going to consider further. 

A. 3 Physical applications of GR 

1. Autoparallel vector fields and 1-forms 

In nonrelativistic and relativistic mechanics the vector fields X on a manifold M are the local 
representatives (velocity vectors) of the evolution trajectories for point-like objects. The condition 



78 



that a particle is free is mathematically represented by the requirement that the corresponding vector 
field X is autoparallel with respect to a given connection V (covariant derivative) in TM: 

i(X)VX = 0, or in components, X a V a X^ + T» V X° X v = 0. (91) 

In view of the physical interpretation of X as velocity vector field the usual latter used instead of 
X is u. The above equation (91) presents a system of nonlinear partial differential equations for the 
components X^, or u^. When reduced to 1-dimensional submanifold which is parametrized locally by 
the appropriately chosen parameter s, (91) gives a system of ordinary differential equations: 

j2 h j v j v 

< 92 » 

and (92) are known as ODE defining the geodesic (with respect to T) lines in M. When M is 
riemannian with metric tensor g and T the corresponding Levi-Civita connection, i.e. Vg = and 
Ty a = Ta U , then the solutions of (92) give the extreme (shortest or longest) distance ds between 
the two points a, b € M, so (92) are equivalent to 




dx» dx v \ () 
ds ds J 

A system of particles that move along the solutions to (92) with g-the Minkowski metric and diiu^r^- > 
0, is said to form an inertial frame of reference. 

It is interesting to note that the system (91) has (3+l)-soliton-like (even spatially finite) solutions 
on Minkowski space-time [5]. In fact, in canonical coordinates (x 1 , x 2 , x 3 , x 4 ) = (x,y,z,£ = ct) let 

= (0, 0, ±-/, /) be the components of u, where < v = const < c, and c is the velocity of light, so 

^ < 1 and u a u a =^1 — f 2 > 0. Then every function / of the kind 

f(x, y, z,£) = f ^x, y, a{z =F > a = const, for example a = 



defines a solution to (91). If u a u°' = then equations (91) (X = u), are equivalent to u^{du)^ v = 
0, where d is the exterior derivative. In fact, since the connection used is riemannian, we have 
= V M | [u v u v ) = u v \i ' pUy, so the relation u u X/ ' yU^ — u u X7 ^Uy = holds and is obviously equal to 
u^{du)^ v = 0. The soliton-like solution is defined by u = (0, 0, ±/, /) where the function / is of the 
form 

f{x,y,z,C) = f(x,y,zT&- 

Clearly, for every autoparallel vector field u (or one-form u) there exists a canonical coordinate system 
on the Minkowski space-time, in which u takes such a simple form: = (0, 0, af, f),a = const. The 
dependence of / on the three spatial coordinates (x, y, z) is arbitrary , so it is allowed to be chosen 
soliton-like and, even, finite. Let now p be the mass-energy density function, so that V CT (pu a ) = 
gives the mass-energy conservation, i.e. the function p defines those properties of our physical system 
which identify the system during its evolution. In this way the tensor conservation law 

V a (pu a u>*) = {V a pu a )u» + pu a V a u^ = 

describes the two aspects of the physical system: its dynamics through equations (91) and its mass- 
energy conservation properties. 

The properties described give a connection between free point-like objects and (3+1) soliton-like 
autoparallel vector fields on Minkowski space-time. Moreover, they suggest that extended free ob- 
jects with more complicated space-time dynamical structure may be described by some appropriately 
generalized concept of autoparallel mathematical objects. 
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2. Electrodynamics 

2.1 Maxwell equations 

The Maxwell equations dF = 0, d * F = in their 4-dimensional formulation on Minkowski space- 
time (M, 77), sign(rf) = (— , — , — , +) and the Hodge * is defined by rj, make use of the exterior derivative 
as a differential operator. The field has, in general, 2 components (F,*F), so the interesting bundle 
is A 2 (M) (g V, where V is a real 2-dimensional vector space. Hence the adequate mathematical field 
will look like SI = F g) e± + *F g> where (ei, 62) is a basis of V. The exterior derivative acts on SI 
as: dO = dF (g ei + d * F g) e2 , and the equation dfl = gives the vacuum Maxwell equations. 

In order to interpret in terms of the above given general view (GR) on parallel objects with respect 
to given sections of vector bundles and differential operators we consider the sections (see the above 
introduced notation) (1 <g 1, SI <g 1) and the differential operator d. Hence, the GR acts as follows: 

<p; d)(l <g> 1, SI <g> 1) = <p)(l g> 1, dft <g 1) = (l.dfi g> 1.1) = dSl. 

The corresponding ($,</?; d)-parallelism leads to dSl = 0. In presence of electric j and magnetic m 
currents, considered as 3-forms, the parallelism condition does not hold and on the right-hand side 
we'll have non-zero term, so the full condition is 

<p)[l <g 1, (dF <g> ei + d * F <g e 2 ) 1] = <£>)[1 <g 1, (m <g ei +j <g e 2 ) g> 1] (93) 

The case m = 0, F = dA is, obviously a special case. 

2.2 Extended Maxwell equations 

The extended Maxwell equations (on Minkowski space-time) in vacuum read [6]- [7]: 

F A *dF = 0, (*F) A (*d * F) = 0, F A (*d * F) + A (*dF) = (94) 

They may be expressed through the GR in the following way. On (M, 77) we have the bijection between 
A 2 (TM) and A 2 (T*M) defined by 7?, which we denote by F <-> F. So, equations (94) are equivalent 
to _ _ 

i(F)dF = 0, i(*F)d*F = 0, i(F)d * F + i(*F)dF = 0. 

We consider the sections = F (g) ei + *F (g e 2 and f2 = F (g ei + *F (g e 2 with the differential operator 
d. The maps $ and if are defined as: $ is the substitution operator i, and 93 = V is the symmetrized 
tensor product in V. So we obtain 

($, v?; d)(F (g ei + *F (g e 2 , F (g ei + *F (g e 2 ) 

— ~ (95) 

= i(F)dF ei V ei + i(*F)d * F e 2 V e 2 + (i(F)d * F + z(*F)dF) ei V e 2 = 0. 

Equations (95) may be written down also as ((i, V)fi)df2 = 0. 

Equations (95) are physically interpreted as describing locally the intrinsic energy-momentum 
exchange between the two components F and *F of SI: the first two equations i(F)dF = and 
z(*F)d * F = say that every component keeps locally its energy-momentum, and the third equation 
7(F)d*F + i(*F)dF = says (in accordance with the first two) that if F transfers energy-momentum 
to *F, then *F transfers the same quantity of energy-momentum to F. 

If the field exchanges (loses or gains) energy-momentum with some external systems, Extended 
Electrodynamics describes the potential abilities of the external systems to gain or lose energy- 
momentum from the field by means of 4 one-forms (currents) J a ,a = 1,2,3,4, and explicitly the 
exchange is given by [7] 

i(F)dF = i(Ji)F, i(*F)d * F = i(J 2 )F, i(F)d * F + i(*F)dF = i(J 3 )F + i(J 4 ) * F. (96) 
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It is additionally assumed that every couple (J a , </{,) defines a completely integrable Pfaff system, i.e. 
the following equations hold: 

dJ a A J„ A J fe = 0, a, 6 = 1,..., 4. (97) 

The system (95) has (3+l)-localized photon-like (massless) solutions, and the system (96)-(97) admits 
a large family of (3+l)-soliton solutions [7]. 

3. Yang-Mills theory 

3.1 Yang-Mills equations 

In this case the field is a connection, represented locally by its connection form uj G A 1 (M) ® g, 
where is the Lie algebra of the corresponding Lie group G. If D is the corresponding covariant 
derivative, and £1 = Dtu is the curvature , then Yang-Mills equations read D * f2 = 0. The formal 
difference with the Maxwell case is that G may NOT be commutative, and may have, in general, 
arbitrary finite dimension. So, the two sections are 1 <8> 1 and *fi <8> 1, the maps <3? and (p are product 
of functions and the differential operator is D. So, we may write 

($, <p; D)(l ® 1, *fi ® 1) = D * ® 1 = 0. (98) 

Of course, equations (98) are always coupled to the Bianchi identity DS1 = 0. 

3.2 Extended Yang-Mills equations 

The extended Yang-Mills equations are written down in analogy with the extended Maxwell equa- 
tions. The field of interest is an arbitrary 2-form \& on (M, rf) with values in a Lie algebra g, dim(g) = r. 
If {Ei}, i = 1, 2, . . . , r is a basis of g we have ^ = tp l <8> Ei and * = ip % ® E{. The map $ is the substi- 
tution operator, the map <p is the corresponding Lie product [, ], and the differential operator is the 
exterior covariant derivative with respect to a given connection uj: = d^ + [oj, We obtain 

($, tp; T>W ® Ei, & ® Ej) = i(^)(dV> m + c^' A V fc C$) ® [S ro , £i] = 0, (99) 

where CJl are the corresponding structure constants. If the connection is the trivial one, then uj = 
and D — > d, so, (99) gives 

i(^)d^ Cg ® E k = 0. (100) 

If, in addition, instead of [, ] we assume for 92 some bilinear map / : g x g — ► g, such that in the basis 
/ is given by f(Ei, Ei) = Ei, and f(Ei, Ej) = for i / j the last relation reads 

itVfyd^* ® -Ej = 0, t = l,2,...,r. (101) 

The last equations (101) define the components as independent 2-forms (of course ip l may be 
arbitrary p- forms). If the bilinear map (f is chosen to be the symmetrized tensor product V:gxg^ 
g V g, we obtain 

i(V^)dV> j ®EiVEj = 0, i S j = 1, • • • , r. (102) 

Equations (100) and (102) may be used to model bilinear interaction among the components of ^. If 
the components of i(^*)d^ ® Ei V Ej have the physical sense of energy-momentum exchange we may 
say that every component ip % gets locally as much energy- momentum from ip 3 as it gives to it. Since 
C\j = —C^, equations (100) consider only the case i < j, while equations (102) consider i < j, in fact, 
for every i,j = 1, 2, . . . , r we obtain from (102) 

iffidip' = 0, and i{^ l )d^ + i(^')d^ = 0. 

Clearly, these last equations may be considered as a natural generalization of equations (94), so spatial 
soliton-like solutions are expectable. 



81 



4. General Relativity 

In General Relativity the field function of interest is in a definite sense identified with a Lorentz 
pseudometric g on a 4-dimensional manifold, and only those g are considered as appropriate to describe 
the real gravitational fields in vacuum which satisfy the equations R^ u = 0, where are the 
components of the Ricci tensor. The main mathematical object which detects possible gravity is the 
Riemann curvature tensor R ail $ v , which is a second order nonlinear differential operator R : g — > R(g). 
We define the map $ to be the contraction, or taking a trace: 

$ : (g a /3, Ran,/3u) = 9 al3 Rafi,f3v = R^u, 

so it is obviously bilinear. The map <p is a product of functions, so the GR gives 

{*,<p;R)(jg®l,g) = ® 1] = Q\g,R(g)]®<p(l,l) = Ric[R(g)} = 0. (103) 

In presence of matter fields \l/ a , a = 1, 2, . . . , r, the system of equations is 

«(v-5«-r)-o. 

It is easily obtained through the GR if we modify the differential operator R aph ^ v to 

Raixfiv — 2 (Tapg^y + T^gap - T av g^p — T^pg av ) + — (g a /3g/iu - gaug^p) T, 

where k is the gravitational constant, T txv {^ a ) = T Vfl {^ a ) is the corresponding stress energy momen- 
tum tensor, and T = g^T^. 

5. Schrodinger equation 

The object of interest in this case is a map * : E 4 — ► C, and M 4 = M 3 x R is parametrized by the 
canonical coordinates (x,y,z;t), where t is the (absolute) time "coordinate". The operator D used 
here is 

d 

D = ih— - H, 

dt 

where H is the corresponding hamiltonian. The maps and ip are products of functions, so the GR 
gives 

($,¥>;D)(1®1,¥) = (l® (ih^-WSfyj 01 = 0. (104) 

6. Dirac equation 

The original free Dirac equation on the Minkowski space-time (M, rj) makes use of the following 
objects: C 4 - the canonical 4-dimensional complex vector space, L C 4-the space of C-linear maps 
C 4 -► C 4 , * e Sec{M x C 4 ), 7 G Sec(T*M <g> L&), and the usual differential d : ^ <g> a -► <ty* <g> e*, 
where {e«}, i = 1,2, 3, 4, is a basis of C 4 . We identify further L C 4 with (C 4 )* C 4 and if {e 1 } is a basis 
of (C 4 )*, dual to {ei}, we have the basis e % ® ej of L C 4. Hence, we may write 

7 = ^T.dxf ® (e* ® ej ), 

and 

7(*) = 7^A M ® ® e i)(^ fe ® e fc ) 

= -yji^ ® V fc < e\ e fc > ej = ^dx" ® ^*<S£ ej = -fapdaf ® ej. 
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The 4 matrices 7^ satisfy ^i^v+lvlii = 2?7 (tu/ zc/c4, so they are nondegenerate: det^^) 7^ 0, \i = 1, 2, 3, 4, 
and we can find (7^) _1 and introduce 7 _1 by 

7" 1 = ((7^)- 1 )l^®(e J ®e J ) 

We introduce now the differential operators V ± : Sec(M x C 4 ) -> S , ec(T*M(g)C 4 ) through the formula: 
P ± = id ± ^mj -1 , i = \/— T, m € M. The corresponding maps are: $ = 77, </? : L C 4 x C 4 — > C 4 given 
by <^(a* <8> /3, p) =< a*,p> (3. We obtain 



(<&, P ± )( 7 , *) = (cD, ^(7^ ® (e* ® e,), ® e k ± -m( 7i 7 1 )^^ (e r ® e s )^ m e m ) 

• <9v; fc 1 
= ^UTTUVidx^dx") < e\e k > ej ± -m^^y^^dx^ ,dx u ) < e\e s > ej 

In terms of extended parallelism we can say that the Dirac equation is equivalent to the requirement 
the section \& G Sec(M x C 4 ) to be (77, y?; D ± )-parallel with respect to the given 7 6 Sec(M x L C 4). 
Finally, in presence of external electromagnetic field A = A^dx^ 1 the differential operators modify 
to 2) ± = (id — eA) ± 1 jm7~ 1 , where e is the electron charge. 



A. 4 Conclusion 

It was shown that the GR, defined by relation (85), naturally generalizes the geometrical concept of 
parallelism, and that it may be successfully used as a unified tool to represent formally important 
equations in geometry and theoretical physics. If * is the object of interest then the GR specifies 
mainly the following things: the change of the object with respect to which we consider the 
change, the "projection" of the change on ^1 through the bilinear map <£, and the bilinear map <p 
determines the space where the final object takes values. When f = $1 we may speak about autopar- 
allel objects in the corresponding sense, and in this case, as well as when the differential operator D 
depends on ^ and its derivatives, we obtain nonlinear equation(s). In most of the examples considered 
the main differential operator used was the exterior derivative d and its covariant generalization. 

In the case of vector fields and one-forms on the Minkowski space-time we recalled our previous 
result that among the corresponding autoparallel vector fields there are finite (3+1) soliton-like ones, 
time-like, as well as isotropic. This is due to the fact that the trajectories of these fields define straight 
lines, so their "transverse" components should be zero if the spatially finite (or localized) configuration 
must move along these trajectories. Moreover, the corresponding equation can be easily modified to 
be interpreted as local energy-momentum conservation relation. 

It was further shown that Maxwell vacuum equations appear as d-parallel with somewhat trivial 
projection procedure. This determines their linear nature and leads to the lack of spatial soliton-like 
solutions. The vacuum extended Maxwell equations are naturally cast through the GR in the form 
of autoparallel (nonlinear) equations, and, as it was shown in our former works and in this review, 
they admit photon-like (3+1) spatially finite solutions, and some of them admit naturally defined spin 
properties [8]. The general extended Maxwell equations (96) may also be given such a form if we 
replace the operator d with (d — i(J a )), and (3+1) soliton solutions of this system were also found [7]. 

The Yang-Mills equations were also described in this way. The corresponding extended Yang-Mills 
equations are expected to give spatial soliton solutions. The Einstein equations of General Relativity 
also admit such a formulation. 
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In quantum physics the Schrodinger equation admits the "parallel" formulation with trivial pro- 
jection. A bit more complicated was to put the Dirac equation in this formulation, and this is due to 
the bit more complicated mathematical structure of 7 = ^^dx^ ® (e l <8> ej). 

These important examples make us think that the introduced here extended concept for ($, tp; D)- 
parallel/ 'conservative objects as a natural generalization of the existing geometrical concept for V- 
parallel objects, may be successfully used in various directions, in particular, in searching for appro- 
priate nonlinearizations of the existing linear equations in theoretical and mathematical physics. It 
may also turn out to find it useful in looking for appropriate lagrangians in some cases. In our view, 
this is due to the fact that it expresses in a unified manner the dual change-conservation nature of 
the physical objects. 



References to the Appendix 



[1]. Greub, W. Multilinear Algebra, Springer, 1967 

[2]. Lichnerowicz, A., Les Relations Integrales d'Invariance et leurs Applications a la Dinamique, 

Bull. Sc. Math., 70 (1946), 82-95 

[3]. Cartan, E., Legons sur les Invariants integraux, Hermann & Fils, Paris, 1922 

[4]. Godbillon, C, Geometrie Differentielle et Mechanique Analytique, Hermann, Paris, 1969 

[5]. Donev, S., Autoclosed Differential Forms and (3+l)-Solitary Waves, Bulg. J. Phys., 15 

(1988), 419-426 

[6]. Donev, S., Tashkova, M., Energy-momentum Directed Nonlinearization of Maxwell's Pure 
Field Equations, Proc. R. Soc. Lond. A, 443 (1993), 301-312 

[7]. Donev, S., Tashkova, M., Energy-momentum Directed Nonlinearization of Maxwell's Equa- 
tions in the Case of a Continuous Madia, Proc. R. Soc. Lond. A, 443 (1995), 281-291 

[8]. Donev, S., Screw Photon-like (3+l)-Solitons in Extended Electrodynamics, The European 
Physical Journal "B", vol.29, No.2, 233-237 (2002). 



84 



